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ABSTRACT 

The  linearized  angular  motion  of  a  symmetric  missile  is  developed  in  some 
detail  with  some  consideration  of  the  rolling  motion,  drag,  aerodynamic  jump 
and  the  effect  of  varying  coefficients.  The  tricyclic  motion  of  a  missile 
with  misaligned  control  surfaces  is  briefly  considered.  This  linear  theory  is, 
then,  applied  to  the  analysis  of  ballistic  range  data. 

Next,  simple  cubic  nonlinearities  in  static  moment  and  Magnus  moment  are 
treated  by  a  quasi-linear  analysis  and  these  cubic  coefficients  obtained  from 
ballistic  range  flight  data.  More  generalized  relations  for  arbitrary  symmetric 


nonlinear  terras  are  derived  and  their  use  in  the  construction  of  amplitude" 
planes  indicated.  These  amplitude  planes  have  proven  to  be  quite  useful  for  the 
prediction  of  missile  flight  performance.  Finally,  the  influence  of  strongly 
nonlinear  static  moments  is  determined  by  a  perturbation  method  which  makes  use 
of  two  quasi-constants  of  the  motion  --  total  energy  and  angular  momentum. 
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preface 


iiaterial  of1  this  report  was  prepared  in  the  form  of  five  chapters  of 
a  book,  The  Theory  and  Application  of  Ballistic  Ranges,  which  is  being  written 
under  the  joint  authorship  of  G.  V.  Bull  and  C.  H.  Murphy  for  publication  by 
Prentice-Hall.  These  chapters  were  written  by  C*  K.  Murphy  and  form  a  complete 
treatment  of  the  angular  motion  of  symmetric  missiles  acted  on  by  both  linear 
and  nonlinear  forces  and  moments  with  an  emphasis  on  the  analysis  of  ballistic 
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publications,  this  unified  form  is  felt  to  be  sufficiently  useful  for  separate 
publication  as  a  BRL  Report. 


Each  section  is  numbered  by  a  decimal  which  indicates  its  chapter  and 
position  in  that  chapter.  Thus  section  6.2  is  the  second  section  of  Chapter 
VI.  Each  equation  could  then  be  identified  by  its  section  and  position  in 
that  section.  Equation  (6. 2.k)  is  the  fourth  equation  in  section  6.2.  For 
simplicity,  however,  the  number  identifying  the  chapter  has  been  omitted  from 
all  equation  numbers  except  when  a  reference  is  made  to  an  equation  in  a 
different  chapter. 
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5-1  Introduction 

Historically,  the  first  interest  of  the  ballistician  was  the  drag  force 
exerted  on  his  projectile  by  the  air.  Thus,  the  first  measurements  made  on 
ballistic  ranges  were  drag  measurements.  These  measurements  also  have  the 
attractive  characteristic  of  being  direct  and  conceptually  very  simple.  Since 
the  rolling  motion  has  a  similar  data  analysis,  it  will  also  be  considered  in 
this  chapter. 


5.2  The  Drag  Equation 

The  drag  force,  D.  is  usually  expressed  in  terms  of  a  drag  coefficient  by 


~  ^  4-U-v  fnl  InTHnrr 

means  U1  one  Iuj.±uw4.ii^  ucx  iiunuii. 


2 

D  =  (±/d)pV  £><J 


(2.1) 


where  p  is  air  density 
V  is  velocity 
S  is  reference  area  and 
CD  is  the  drag  coefficient 

The  velocity  of  the  missile  must  then  satisfy  the  equation 


/ 1  /r\ 


t  O 


2 


\ 
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In  most  ballistic  range  work,  distance  is  a  more  convenient  independent 
variable  than  time.  For  pitching  and  yawing  motion,  it  also  introduces  an 
important  simplification  in  the  equations  of  motion  (see  Chapter  VI) .  If  l 
is  used  to  denote  a  characteristic  length,  a  dimensionless  arclength  may  be 
defined  by  the  relation 


=  I  V  dt 

*  J 

0 


(2.3) 


Tn  most  cases,  the  reference  area  is  taken  to  be  the  maximum  bodv  cross 


sectional  area  and  the  reference  length  the  maximum  body  diameter.  F Or  this 


_ _  —XV  XI _ _ _ V  . 


new  variable,  Equation  (2.2)  transforms  to 
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t*y\'  -  L  - 

' >  v  " 


/pSJ\ 

vary 


CD  5  -  CD 


I  \ 


where  prime  denotes  derivative  with  respect  to  s. 


_>C  B 

pkJ 


The  factor,  ,  will  appear  a  number  of  times  in  later  chapters.  Since 

the  volume  of  the  missile  usually  varies  from  (2/j)Si  for  a  sphere  to  lOSi  for 

a  missile  with  high  fineness  ratio,  the  inverse  of  this  factor  is  a  measure  of 

the  model  mean  density  relative  to  that  of  air.  The  relative  density  of  plastic 

3  5 

is  10  and  that  of  steel  is  8  x  10  .  Thus  this  relative  density  factor  has  a 
maximum  of  10  ^  for  a  plastic  sphere  and  is  usually  of  the  order  of  10-1+. 

For  constant  C_,  Equation  (2.4)  has  a  solution  of 


V  =  V  e 
o 


* 

-CDS 


(2.5) 


Integrating  Equation  (2.5),  a  distance-time  relation  results 


t  =  t  + 
o 


* 

CDS 

Ufa  D 


5T 

VoCD 


- 


(2.6) 


=  t  + 
o 


(t) 


* 


2 

s  + 


It  is  important  to  note  that  iu  »ny  data  analysis  the  origin  for  s  must  be 
selected  as  the  midpoint  of  the  timing  stations  so  that  the  linear  and  the 
quadratic  coefficients  are  related  to  well  determined  velocities  and  retardations 


There  are  two  approaches  to  the  data  analysis  problem.  One  is  that  of  the 

5-3  5-1 

engineer  and  the  other  is  that  of  the  mathematician  .  Fortunately,  the 

results  of  the  two  methods  are  practically  identical.  In  the  first  method,  the 

engineer  makes  use  of  his  love  for  plotting  data  and  exercising  his  personal  o- 

pinion  as  to  the  best  fit.  This  is  done  by  use  of  Equation  (2.4).  He  plots  the 

logarithm  of  the  average  velocity  between  a  pair  of  timing  points  versus  the 

location  of  the  midpoint  of  the  interval.  He  then  draws  a  straight  line  through 

the  points  and  computes  the  drag  coefficient  from  its  slope. 
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The  mathematician  displays  his  characteristic  fear  of  plotted  data  and 
human  decision  making.  He  fits  his  time-position  data  by  a  cubic  least  squares 
process.  For  the  usual  flat  trajectories,  s  can  be  well  approximated  by  x^/ l 
where  x.,  is  distance  measured  downrange. 


t  =  io  +  a1x1  +  +  a;x3 


(2-7) 


where 


a  =  V 
1  o 


-1 


pS 

a2  TuriT 


PS  ) 


(?) 


(2.6) 


The  cubic  coefficient  absorbs  any  Mach  number  variation  in  as  well  as  errors 
in  the  expansion  of  the  exponential  in  Equation  (2.6).  In  those  cases  where  the 
cubic  fits  well 
first  approach. 


cubic  fits  well,  the  resulting  Cp  differs  very  little  from  that  obtained  from  the 


The  accuracy  of  a  measured  drag  coefficient  may  be  estimated  statistically 
when  a  least  squares  reduction  has  been  performed  or  determined  by  the  physical 
size  of  the  drag  "disturbance".  Since  the  velocity  is  so  much  better  determined 


oiic  errw  in  is  directly  proportional  to  the  error  in 


ag  alone.  This  error  may  be  estimated  from  the  corresponding  diagonal  term  of 
the  inverse  matrix  formed  from  the  least  squares  normal  equations. 


(2.9) 


where 


e 


At  =  t  -  t 

calc.  measured 

n  =  number  of  measurements  and 

is  the  diagonal  element  corresponding  to 

of  the  inverse  matrix  of  the  normal  equations. 


The  use  of  the  size  of  the  measurable  effect  of  drag  introduces  physical 
insight  into  this  error  analysis  which  is  sadly  lacking  in  the  above  cold 
statistical  formula.  If  drag  were  not  present,  the  missile  would  travel  at 
constant  velocity  indefinitely.  The  drag  changes  the  time  of  passage  of  the 
model  at  the  various  observation  points.  The  maximum  drag-induced  time 
decrement  with  respect  to  constant  velocity  flight  at  the  average  velocity, 

V  ,  occurs  at  the  ends  of  the  range.  If  VQ  is  measured  at  the  midpoint  of  the 
range*  and  the  range  length  is  IL,  the  endpoints  are  located  at  s  =  —  (L/2). 
Maximum  time  decrement  may  then  be  computed  for  Equation  (2.6) 


* 

CDS 

~(e  -  1) 
If 


V  c_ 

o  D 


is 

V. 


*  2 
v 


psi2 

1  Art]y 

O 


(2.10) 


The  engineer  can  now  make  the  reasonable  assumption  that  the  percentage 
error  in  is  given  by  the  ratio  of  measurement  error  to  the  maximum  time 
decrement,  e/d  .  Usually,  the  time  measurements  are  less  accurate  than 
distance  measurements  and  maximum  estimated  error  in  time  can  be  used.  Should 
the  range  survey  be  the  limiting  factor,  the  distance  error  can  be  converted 
to  an  equivalent  time  error  by  division  by  VQ  and  the  result  compared  with  d^_. 

The  presence  of  length  squared  in  the  time  decrement  shows  the  value  of 
long  flight  paths.  Exploitation  of  this  property  is  limited,  however,  by  the 
presence  of  Mach  number  dependence  in  the  drag  coefficient.  It  is  also 


The  midpoint  is  not  the  point  at  which  the  actual  velocity  is  exactly  equal 
to  the  average  velocity.  A  small  correction  can  be  employed  to  eliminate  this 

minor  difficulty"' 
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interesting  to  note  that  higher  velocity  tends  to  degrade  drag  accuracy.  In 

any  event,  it  is  the  common  experience  that  ballistic  ranges  with  lengths  of 

one  thousand  calibers  can  attain  better  than  l/2  nercent  accuracy  in  C  at 

'  ^  D 

high  supersonic  velocities. 


The  question  of  an  optimum  distribution  of  timing  stations  for 


Lver 


length  can  be  easily  answered.  Since  the  drag  appears  in  the  quadratic  term 
of  a  parabola,  it  is  determined  by  the  curvature  of  the  time-distance  curve 
with  respect  to  the  midrange  point.  Thus,  one  half  of  the  observations  should 
be  at  the  center  and  one  half  at  the  ends.  It  can  be  shown  theoretically^ 
that  this  distribution  of  l/h  -  l/2  -  l/k  will  yield  15  -  20  percent  less  error 
in  C-p,  than  a  symmetric  distribution.  Such  a  small  imnrovement  is  of  little 

U 
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5. 5  Effect  of  Pitching  and  Yawing  Motion 

If  the  angular  motion  is  not  small,  variations  in  the  drag  coefficient 
are  observed.  Usually  this  functional  dependence  on  angle  is  well  represented 
by  a  quadratic  function. 


C 


D 


(5-1) 


where 


5C 


a  =  angle  of  attack 
3  =  angle  of  sideslip 

2 

As  will  be  shown  in  the  next  chapter,  the  variation  of  5  for  a  symmetric 
missile  is  described  by  the  equation 

62  =  +  SKjK^os  %  (3.2) 

where  K.  =  K,  e^js 
J  Jo 

0 s  =  0.  +  0  .s 

'  J  •  J 

0  =  P-l  -  4 

f 

and  Xj  and  0^  are  constants. 


For  a  quadratic  drag  coefficient,  Equation  (2.4)  can  be  integrated  once 
to  yield 


P  O  \ 


m 


r  -  -  (sr)  (%  s  +cDr2  J  •*) 


(3-3) 


This  can  be  expressed  in  exponential  form,  expanded  in  a  power  series,  and 
integrated  again 

/PSi2C  psi2c 

1 t  +f  s  -  1-inr)  s  +  a  »  r<s>  + 

O  '  o  /  o 


s  q 

where  l(s)  =  F  F  8  2 


(3.4) 


J 

o  o 


5  dr  dq 


The  distance  over  which  drag  deceleration  is  measured  is  at  least  several  times 

2 

longer  than  the  period  of  8  amd,  hence,  the  influence  of  the  periodic  part  of 

2 

8  on  l(s)  can  be  neglected.  If  the  effects  of  damping  are  also  neglected, 

I(s)  assumes  the  simple  form 


Tfq'i  r  £_  (v2  +  v?  1 

o  '"10  ^20 ; 


(3.5) 


and  ( —  \  =  c  *  C  0  (K2  +  ) 

V*/  \aiJ  Do  \z  10  *>' 


(3.6) 


At  first  glance,  it  seems  that  Equation  (3*6)  could  have  been  directly 
derived  by  replacing  the  constant  in  Equation  (2,8)  by  the  average  value  of 
C^,  i.e. , 


"D  average  D  D~2 


L/2 

i  r  2 

where  8  =  —  /  8  ds . 

L  Ji 
-L/2 


(3.7) 
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It  Gc'cms  Intuitively  obvious  that  fitting  a  cubic  to  the  time-distance  data 
should  have  this  effect  of  "averaging"  the  drag  coefficient  in  this  way.  It 
is  clear  that  for  the  case  of  no  damping  this  is  correct.  As  we  shall  see, 
it  is  not  correct  for  large  damping. 

The  least  squares  process  includes  in  a„  the  "best"  quadratic  approximation 
of  I(s).  Thus  an  average  8  is  defined  by  the  requirement  that 


L/2  r  2  2 

b  f  s  ~2 

j  l(s)  -  —  &  ds  is  a  minimum. 

■L/2 


A  relation  for  5  may  be  obtained  by  differentiating  the  above  with  respect  to 
8  ,  setting  the  result  equal  to  zero,  and  solving  for  8  . 


(L/2)' 


U 


s  l(s)  ds 


(5.8) 


-  4  +  4 44  +  44> +  ••• 


The  mean  squared  angle  defined  after  Equation  (5*7)  has  the  expansion 


^  ■  4  +  4  44  *44> +  ••• 


(  ^  Q  ^ 

w  y ) 


2  ”5 

Thus,  the  proper  average  value  of  8  is  not  the  same  as  8  and  the  more  precise 
form  of  Equation  (3.6)  is 


ia  11\  -  n 

pS  \B1J  -  -Dr 


ran^e 


„  2 
V 


(3.10) 


According  to  Equation  (3.10),  the  drag  coefficient  obtained  from  the  range 


data  reduction  is  a  linear  function  of  8  .  Thus,  if  CT 


is  plotted  versus 


range 


82  =  K^\  +  for  different  model  tests  at  the  same  Mach  number,  the  points 
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should  fall  on  a  straight  line.  This  technique  has  been  used  for  a  number  of 

years  to  determine  and  C„  0  for  a  large  variety  of  symmetric  configurations 

u  u~  2  ~ 

o  o 

with  complete  success. 


In  Figure  5.1,  the  data  for  a  large  angle  of  attack  program' 


5-5 


is  plotted 


as  an  illustration  of  this  technique.  In  this  program,  a  20mm  bullet  was 

o  ^2 

tested  at  angles  up  to  25  and  a  good  linear  correlation  of  and  5  was 

obtained.  C  „  was  determined  from  the  slope  of  the  line  to  be  4.63.  This 

value  agreed  very  well  with  a  wind  tunnel  measurement  of  4.71. 

As  can  be  seen  from  a  comparison  of  liquations  (3*8)  and  (3*9) >  the 

? 

distinction  in  the  proper  average  5  is  rather  minor.  The  determination  of 
the  proper  average  6^  is  much  more  important  for  other  angle  dependent 
aerodynamic  coefficients.  It  is  for  this  reason  that  the  concept  of  6  was 
introduced  here. 


5-4  Experimental  Results 

The  first  systematic  ballistic  range  drag  measurements  were  of  spheres. 

Charters  and  Thomas^  tested  9/l6"  spheres  at  Mach  numbers  varying  from  0.3 

to  4.0  in  the  BRL  Aerodynamics  Range.  By  tests  in  the  NOL  Pressurized 

5-7 

Ballistics  Range,  May  and  Witt  were  able  to  vary  the  Reynolds  number  by 
more  than  a  hundred  for  a  fixed  Mach  number.  The  results  of  this  program  are 
shown  in  Figure  5*2,  which  is  a  contour  plot  of  on  the  Mach  number  -  Reynolds 
number  plane. 

The  presence  of  a  saddle  point  In  the  vicinity  of  M  =  2.0  and  R  =  30.000 

e  '  ' 

is  clearly  shown  in  this  figure.  In  passing  along  the  line  Rg  =  30,000  a 
maximum  is  found  near  this  point,  while  a  minimum  is  encountered  on  the  line 
M  =  2.0. 


The  influence  of  wall  temperature  on  the  skin  friction  has  been  studied  in 
a  quite  elegant  manner  by  Sommer  and  Short'’”0.  The  test  vehicles  were  spin- 
stabilized  hollow  cylinders  with  fineness  ratios  of  either  1,4  or  1.8.  The 
nose  contours  were  either  a  double  wedge  or  a  double  circular  arc.  (see 
Figures  5-3  and  5.4)  Similar  cylinders  with  a  fineness  ratio  of  0.4  were  used 
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as  taro  models.  The  difference  in  drag  coefficients  of  the  Lest  models  and 
the  taro  models  is,  except  for  small  corrections,  a  measure  of  the  average 
skin  friction  drag  of  the  added  length  of  the  test  model.  Boundary  layer 
trips  were  used  to  ensure  the  presence  of  turbulent  boundary  layers. 


The  tests  were  made  in  the  Ames  Free  Flight  Wind  T'mnei  at  Mach  number 
of  2.8  and  3-9  in  still  air  and  7*2  in  a  Mach  number  2  air  stream.  Since 
the  model  flight  time  was  less  than  .02  sec,  its  wall-temperature  heat  rise 
was  less  than  45°F.  Thus,  the  ratio  of  wall-temperature  to  boundary- layer-edge 
temperature,  Tu/Tn ,  was  determined  by  the  local  temperature  in  the  tunnel, 
i.e.,  1.05  for  otill  air  and  1.3  or  1.7  for  air  flow.  Since  calculated  recovery 
temperatures  were  two  and  a  half  to  six  times  greater  than  the  estimated  wall 
temperatures,  considerable  heat  transfer  was  in  process.  In  Figure  5*5j  the 
results  of  the  program  are  plotted  and  compared  with  a  curve  based  on  zero  heat 


i-  -V-.  —1  M  M  «r  *  d  v-»  S3  4"  '  ^  V“»  /—V  "l  d  n  4"  o 

Ul'tUitji-  til  WlliU  ULUUICX  ULClUCi. 


T  4-  een  V\/->  r?  O  4~  Vi  o  4"  4*  V>  a  v\v»cimcinno  V,  not  +  voncf  nr 
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can  increase  the  skin  friction  by  over  30  percent. 

By  a  somewhat  similar  approach,  base  pressures  have  been  extracted  from 

S-Q 

ballistic  range  drag  data.  Charters  and  Turetsky^  "  determined  the  base  pressure 

for  20°  cones  and  cone  cylinders  by  computing  head  drag  and  skin  friction  drag 

and  subtracting  the  results  from  the  total  drag.  The  base  pressure  may  also 

be  determined  by  a  measurement  of  the  wake  angle  from  a  spark  photograph.  The 

pressure  outside  the  wake  is  computed  by  characteristics  and  assumed  to  be  the 

same  as  at  the  base  of  the  model.  A  comparison  of  ratio  of  base  pressure  to 

free  stream  pressure,  F_/p, ,  as  obtained  from  these  two  methods,  is  given  in 

13  -L 

Table  I. 


TABLE  I 

P_/P,  for  M  =  1.84 
B'  ± 


Model 

Fineness 

Ratio 

Total  Drag 
Method 

Wake  Angle 
Method 

H  r\v\  <2* 

O  Q 

C\  ^ C\ 

n  l.n 

U  Ol 

y 

v  •  yw 

•  ~*y 

Cone-cylinder 

3.2 

0.53 

0.50 

Cone- cylinder 

5.0 

0.6l 

0.59 
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A  final  example  of  ballistic  range  drag  data  is  shown  in  Figures  and 

5*7.  It  can  be  seen  that  the  separated  flow  over  the  nose  spike  of  this  fin- 

stabilized  round  can  assume  two  different  configurations.  The  model  with  the 

larger  separated  flow  has  one-third  less  drag  than  the  identical  shape  with  the 

smaller  separated  region;  This  low  drag  flow  occurs  when  the  initial  launch 

angle  is  less  than  two  or  three  degrees.  For  larger  angles,  the  high  drag  flow 

5-10 

is  established  and  persists 
5.5  The  Roll  Equation 

The  aerodynamic ist  usually  writes  the  roll  moment  due  to  control  deflection 
and  aerodynamic  damping  in  the  form 


“x  = 


(l/2)pV2b 


>C/R6f 


St°i 


(5.1) 


where  b  is  span 

Sfi  is  profile  area  of  deflected  fins 

St  is  total  profile  area  of  fins 

6f  fin  deflection  angle  for  pair  of 
differentially  canted  fins 

n  =  -roll  T'n+'.e 

P  dt - - 


Equation 


f  -i  \ 

W  ■  - 


as  two  reference  areas  and  two  reference  lengths,  b  and  b/2. 
If  the  test  missile  possesses  surfaces  with  two  different  spans,  the  decision 
as  to  which  is  the  reference  length  becomes  quite  difficult.  In  ballistic 
range  work,  it  is  quite  convenient  to  use  one  reference  area,  S,  and  one 
reference  length,  i.  With  this  in  mind,  we  can  redefine  the  roll  moment  as 


lr  /  .  / r-\  \  ..2  4aa 

=  il/z;pv  is 


-  -  pi 

6fui  +  V 


(5.2) 


Conversion  relations  for  the  aerodynamic  coefficients  may  be  obtained  by  a 
comparison  of  Equation  (5.1)  and  (5-2).  All  ballistic  range  data  should  contain 
a  specification  of  the  reference  area  and  reference  length  which  were  used. 


The  equation  for  rolling  motion  is 


•  *  p 

T  fk  —  (  1  /O  \  f*r>  o 

pv  ^  o 


<j«  ~  1 


+  r  c*  j 

o 


(3.  >) 


If  drag  can  not  be  neglected.  Equation  (5*3)  is  linear  with  variable 
coefficients.  For  a  constant  drag  coefficient,  transforming  the  independent 
variable  from  time  to  distance  will  provide  us  with  a  linear  equation  with 
constant  coefficients.  Sin~e  distance  is  the  natural  variable  for  ballistic 
ranges,  this  transformation  is  doubly  attractive,  we  therefore  introduce 
the  variable  s,  which  was  defined  by  Equation  (2.5). 


f  =  22  ,  22  31  .  i0 

v  -  ds  dt  ds  "  v  50 


(5*4) 


'•$■{>)"  [*•»'  ?  v] 


(5*5) 


Substituting  Equation  (5*4  -  5*5)  in  Equation  (5*3),  we  have  the  roll  equation 
in  the  form 


1»  | 

0  +  K  0  -  Kg  =  0 


(5.6) 


nQ  0  1-0 

where  K  =  -  — — ^  k  C .  +  CL 

p  2m  a  £  D 

*—  n  -1 


K  -  £^1 

6  21  bfS 

x  6 


V# 


axial  radius  of  gyration 


The  solution  to  Equation  (5.6)  is 


0  =  0Q  +  0ss  +  A(e  P  -  1) 


(5*7) 
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where  0_ 


K- 

"5 


is  the  steady  state  roll. 


A  =  (01  -  01)  =  (K  -  0 '  K_)  Kl 

t>  up  O  0  p  p 


0_,  01  are  the  conditions  at  the  middle  of 
u  u  the  observed  flight. 

If  K^s  is  small  enough,  the  exponential  in  Equation  (5.7)  can  be  expanded 


a.s  a  cubic 


0  =  0„  +  0‘ s  +  (AKT) 


/  2  K  s  \ 

IS---V 


(5.8) 


Equation  (5.8)  can  be  fitted  by  the  same  least  squares  reduction  that  was  used 


Unfor  innately }  in  order  to  separate 


■f*  AM  -i-  V«  ^  -MO  r-w  in4-  ■!  f  nn-H  rtw  f  O  *7  \  \ 

I  Ox  one  CLTcig  c^uaoluii  ^D^uauiuu  \c^»  [  )  )  . 

Kp  and  K&,  the  cubic  coefficient  in  Equation  (5*8)  must  be  well  determined.  In 
other  words,  K^s  can  not  be  too  small.'  At  the  Ames  Free  Flight  Wind  Tunnel, 
which  has  a  test  section  of  15  feet,  this  determination  has  been  found  to  be 
quite,  difficult"* 

This  difficulty  can  be  avoided  by  attempting  to  measure  these  quantities 

separately.  can  be  determined  fairly  well  from  the  quadratic  tern  when 

Kg  =  0,  i.e.,  for  spinning  bodies  of  revolution  or  missiles  with  uncanted  fins. 

Next,  the  influence  of  K  can  be  eliminated  from  the  auadratic  term  bv  causinc 

p  ..  “  - 

4*  Vi  a  m  7  4*  "v.  rt  n  a  ^  ^  a  ^  aa  ^  7  7  4  v\  rr  7  n  4  V  i  f  m  a  a  ^  V  V  a  a  am  4  am  a  4  7a  a  a  V  a  a  -  r  a  J  ^  7  J  » V-  X 

oilc  mum;!  ou  i  coui  x w  iumug  vciutiuj  ilSckr  one  ociiuci'  ux  one  uubcrv eu.  xxiglll» 


path. 


x 


“5 


/c  r\  \ 

\j*y) 


This  can  be  done  by  properly  prespinning  the  model  in  a  direction  opposite  to 
steady  state  spin  through  the  use  of  a  rifled  gun. 

If  KpS  is  large,  the  exponential  character  of  Equation  (5*7)  is  manifested 
and  both  and  Kg  can  be  determined.  The  usual  least  squares  fitting  process  is 

x  - 

handicapped  by  the  fact  that  the  parameter  appears  in  a  nonlinear  fashion.  An 

■I  f  p  4"  "1  iro  f +^r»fr  r>yAr»oao  la  W  a-mo^ama  maa  ac  no  »m» 

J.OCx^*wx»i-  i.  A.  V  pj  wv-wuu  xo  |  wi  ^  j  u^ocooax  J  t 
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A  qaI  nf*  i  n-j  vfilnac  ck  fk  A  nnrl  IT  o  -r*r»  r*r>n  U  iT’Gfl  for  t‘,h  i  f*.  fJlTOCGnij  , 

o  c  p 

To  obtain  these  values,  we  differentiate  Equation  (5 .7)  and  combine  the  result 
with  Equation  (5*7)  to  eliminate  the  exponential  term. 


0*  +  K  0  -  (K  0f)s  +  (K  A  -  K  0  -  0*) 

v  p's'  v  p  p*o  r s' 


(^.10) 


,  * 

0  can  be  computed  by  first  differences  and  Equation  (5-10)  fitted  by  least 

j'  ,  ,» 

squares  since  it  is  linear  in  the  parameters  K  ,  K  0  ,  and  K  A  -  K  0  -  0 

p  p's  p  do  s 

Values  of  these  parameters  plus  the  measured  value  of  0q  provide  the  required 
set  of  initial  numbers . 

These  values  for  0  ,  0_,  A,  and  K  may  be  placed  in  Equation  (5.7)  and  0 
O  S  P'  ~  '  r 

computed  for  each  measurement  point.  The  difference  between  these  numbers  and 


r*  i^\Trrr\i  ito.)  -Pro*  r* 


the  actual  measurements  will  be  denoted  by  £0  =  0  ,  -  0  ,  .If  the 

'measured  observed 

data  is  reasonably  good,  these  residuals  of  fit,  A0’s,  should  be  less  than  ten 

degrees  (0.2  rad.).  Assuming  that  they  are  due  to  small  errors  in  0  ,  0  ,  A. 

o  s 

and  Kp,  we  differentiate  Equation  (5*7)  with  the  result; 


,  -tv  s  -tv  s 

A0  =  A0  +  sA0  +  (e  p  -  l)AA  -  Ase  p  AK 


(5.11) 


Equation  (5 . 11 ) ,  for  measured  0  and  s ,  are  linear  in  the  unknowns  A0  , 
,•  0 
£0  ,  AA,  and  AK  .  These  equations  may  be  solved  by  least  squares  and  new 
s  p 

parameters  computed  from  the  old  by  adding  the  corresponding  changes,  i.e., 


A0q,  etc.  if  the  corresponding  residuals 


iller,  the  process  is 


successful  and  may  be  repeated  until  no  further  improvement  is  obtained. 
Probable  errors  in  fit  of  less  than  0.5°  are  common. 

5 .6  Experimental  Results 

For  a  spinning  body  of  revolution,  Charters  and  Kent^"^  have  derived  a 
simple  relation  between  the  roll  damping  moment  C  .  and  the  skin  friction  drag 

Ju 

p 

coefficient  .  For  a  basically  cylindrical  body  and  a  reference  length 
"sf 

which  is  the  cylinder  diameter,  the  relation  is 


(6.1) 
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The  derivation  may  be  summarized  by  the  following  steps: 


(l)  For  a  spinning  body  of  revolution,  the  local  velocity  of  the 
surface  relative  to  free  stream  is 


y  _  \r 

V  —  V 

r 


/ 1  .  (  Pd  Y 
V  A  T  l2V  / 


It.* 2) 


where  d  is  body  diameter 

(2)  The  local  skin  friction  force  has  the  direction  of  V  and 

r 

V 

therefore  is  inclined  at  an  angle  of  arccos  —  to  the  axis  of  rotation. 

r 

(3)  The  skin  friction  drag  coefficient  and  the  roll  damping  moment 
may  be  written  in  terms  of  the  integrals 


(l/2)pV2SCD  =  (1/2)  j  p'r 
sf  Q 


iV2  Cf(nd)dx 


(6.3) 


(l/2)pV2SiCjg  =  (1/2)  ^ 


i  ,  ..2 
QP^r 


Cf(nd)dx 


(6.4) 


(4)  For  constant  d  =  l,  Equation  (6.1)  follows  directly  from 


TJ’/n  i  t  o  4"  ■>  o>"i  c  (f,  *  _  £  M 

Uu  w  x  one  y  —  U  ♦  ^  y  * 


The  roll  damping  moment  has  been  measured  for  shell  and  other  bodies 


'  no  r>Q 


-2 

> ;  K  C  .  is  negative  but  oris — thir d.  t, 

Q.  Z 

T) 

X' 

the  drag  coefficient.  Thus,  for  these  shapes  the  roll  per  distance  actually 
increases,  (K  <  0).  C  can  be  measured  to  an  accuracy  of  better  than  10  percent 

p  l-r 

XT  r  -lA  c  12 

and  relation  (6.1)  has  been  verified  by  a  number  of  investigators^  ’  '>~  . 

The  rolling  motion  of  finned  missiles  was  first  studied  on  a  ballistic  range 
5-14 

by  Bolz  and  Nicolaides  .  The  test  vehicle  was  a  cone  cvlinder  of  fineness 


rfifl  A  4on  tri  4  Vl  ennor*£»_r\l  ovwa  f 

i  G,Uxv  uun  "iwu  uuUai,  v-  pj-out  1W1  1U  Cl  UOllUilH  ilUO« 


The  fin  chord  to  body  diameter 


was  one.  As  can  be  seen  from  Figures  5*8  and  5-9 »  the  experimental  results 

verify  the  prediction  of  the  linearized  theory.  In  this  program,  and 

P 

C.  were  determined  to  an  accuracy  of  better  than  2  percent.  (Reference 

6  nd2 

area  was  body  cross  sectional  area,  — and  reference  length  was  body 

diameter,  d.) 
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SHADOWGRAPH  OF  TEST  MODEL  WITH  BOUNDARY- LAYER  TRIP 
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CHAPTER  VI 


LINEARIZED  THEORY  OF  MOTION  CF  SYMMETRIC  MISSILES 
6.1  Introduction 

All  the  early  stability  tests  on  ballistic  ranges  involved  spinning  bodies 
of  revolution  which  might  be  used  as  shell.  In  19^7,  57nm  models  of  the  General 
Electric  Dragonfly  missile  were  fired  at  the  Ballistic  Research  Laboratories  . 
Since  that  time  the  number  of  nonspinning  missiles  tested  has  increased. 

Although  static  stability  was  initially  obtained  by  the  use  of  fins,  in  recent 
times  statically  s  tat  Is  bocliss  of  revolution  which  could,  be  used  as  I  CBM  nose 
cones  have  dominated  this  area  of  ballistic  range  testing.  Traditionally,  «1 1 
stability  tests  have  implied  the  assumption  of  linear  moments,  but  this 
limitation  also  has  been  lifted  in  certain  cases . 

In  this  chapter  we  will  reverse  the  historical  development  and  consider 


tho  thonrot.-ipfll  ]  v  a  rrml  r»QQ^  nf  a  tinnqTT!  riylqq-M 


i iuulcx  a 


pin 


and  the  associated  gyroscopic  and  Magnus  moments  will  be  introduced.  This 
work  will  not  be  linearized  until  very  late  in  the  development  so  that  the 
results  can  be  used  in  Chapter  VIII  on  the  analysis  of  nonlinear  data. 

6.2  Small  Amplitude  Pitching  Motion  of  a  Nonspinning  Missile 

For  an  airplane,  the  pitching  motion  is  the  angular  motion  of  the  airplane's 
thrust  line  in  the  vertical  plane  of  symmetry.  The  yawing  motion  is  then  the 
motion  of  the  thrust  line  perpendicular  to  this  plane  of  symmetry.  When  a 
missile  in  free  flight  is  considered,  these  definitions  have  to  be  altered. 

Since  most  missiles  have  rotationally  symmetric  bodies,  the  thrust  line  can 
be  replaced  by  the  axis  of  symmetry.  Unfortunately,  a  missile  can  possess  a 
number  of  planes  of  symmetry  and  none  of  them  may  be  vertical.  In  this  case, 
the  selection  of  the  plane  of  pitching  motion  is  arbitrary. 

We  leave  the  definition  of  pitching  motion  for  a  missile  without  a  plane 
of  symmetry  to  the  engineer  who  designed  such  an  unusual  configuration.  If  a 
basically  symmetric  missile  is  either  accidentally  deformed  in  some  minor  way 
or  has  intentional  asymmetries  due  to  control  surfaces  at  a  small  angle  of 
deflection,  its  pitching  motion  can,  however,  be  defined  in  terms  of  its  original 
symmetric  shape.  The  influence  of  the  asymmetries  is  incorporated  into  the  total 
aerodynamic  force  and  moment  by  the  introduction  of  small  constant  force  and 
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Before  concentrating  on  the  simple  problem  of  pure  pitching  motion  we  mu.;t 
Tirol  consider  two  coordinate  systems:  a  missile-fixed  system  and  an  earthi- 
fixed  system.  The  missile-fixed  coordinate  axes,  X,  Y,  and  Z,  are  defined  in  the 
following  way  c':  X  -  centerline  of  symmetrical  body,  positive  forward;  Y  - 
perpendicular  to  plane  o>  symmetry,  positive  to  the  right  when  looking  forward*; 

",  -  perpendicular  to  both  the  X  and  Y  axes  and  directed  according  to  the  right 
hand  rule,  i.e.,  downward.  In  a  similar  manner,  earth-fixed  axes  may  be  defined: 

-  perpendicular  to  the  gravity  vector;  Yg  -  perpendicular  to  Xq  axis  and  the 
gravity  vector,  positive  to  the  right  when  looking  in  the  positive  Xg  direction; 

Zg  along  tne  gravity  vector,  positive  downward. 

These  two  axis-systems  are  shown  in  Figure  6.1.  The  relative  position  of 
the  axes  can  be  described  by  the  Eulerian  angles  ty,  0,  0.  The  velocity  vector, 

V,  is  aloo  shown  in  this  figure.  Its  coordinates  in  the  mis  5 i_0 - f* ix0(l  system 
are  u,  v,  w  and  its  orientation  determined  by  the  angle  of  attack,  a,  and  the 
angle  of  sideslip,  p.  Notice  that  positive  a  and  p  are  measured  opposite  to  the 
positive  directions  of  the  Z  and  Y  axes  respectively. 

The  aerodynamic  force  acting  on  a  missile  in  flight  is  usually  expressed  in 
the  missile-fixed  system  with  coordinates  F„,  F„,  F^..  From  these  force  components 
ii mens ionless  coefficients  can  be  obtained. 


-  (1/2) 

PDV  cx 

(2.1) 

=  (1/2) 

pSV^Cy 

(2.2) 

=  (1/2) 

psv^z 

(2.3) 

Tne  aerodynamic  moment  has  components  M^,  My,  My.  To  distinguish  the  corresponding 
coefficients  the  subscripts  £,  m,  n  are  used.  In  conflict  with  aircraft  design 
practice  but  in  accordance  with  our  remarks  of  Section  5.5  we  will  use  only  one 
reference  length. 


If  the  plane  of  syrmetry  is  not  exactly  vertical,  it  is  possible  to  place 
yourself  on  the  axis  of  symmetry  in  a  forward  facing  position  in  the  plane 
of  symmetry  so  that  your  head  is  above  your  feet.  Once  this  is  done  the 
selection  of  the  positive  direction  of  the  Y  axis  is  easy.  If  the  missile 
is  rolling,  this  selection  must  be  made  at  some  instant  in  time  and  the 
direction  fixed  with  respect  to  the  missile  thereafter. 
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R,  =  (1/2)  pV2SZC. 

A  L 

My  =  (1/2)  pV2S«Gm 
\  =  (1/2)  oVr2SiCn 


(2  =  1) 
(?• 5) 
(2.6) 


These  quantities  are  shown  in  Figure  6.2.  The  angular  velocity,  which  is  not 

shown  there,  has  components  p,  q,  r  with  positive  senses  which  are  the  same  as 

those  shown  for  C  C  ,  and  C  ,  respectively. 

£’  m'  n'  ^  J 


For  a  nonspinning  missile  undergoing  pure  pitching  motion,  the  functional 

- ^  n  n  n  ^4-  -u^  ^,4. 

& u  Uc  o  oaocu 

L j  U1 

the  negative  of  C 


For  small  angles,  is  essentially 


Z'  m 

q.  It  is  not  unreasonable  to  expect  the  lateral  force  and 

pitch  moment  to  depend  on  the  actual  motion  in  addition  to  the  flow  parameters, 

Mach  number  and  Reynold's  number.  Although  the  moti  on  can  be  described 

completely  by  a,  q,  and  all  their  time  derivatives  at  some  given  time,  we  will 

only  consider  a,  q,  and  A.  For  small  amplitude  motion,  a  further  assumption  of 

linear  dependence  can  be  made  and  the  exnlicit  relations  for  C  and  C  follow*. 

Z  m 
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(2.7) 
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*  r  (  4?) 

(2.8) 
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<1 

vv  / 

'"4  Vv  ! 

are 

due  to 

asymmetries. 

The  various  important  angles  are  shown  in  Figure  6.3.  The  equations  of 
motion  are  the  drag  Equation  (5.2.2),  the  force  equation*  along  the  inertial 
axis  Zg  and  the  moment  equation  about  the  Y  axis. 

mV  =  -  D  ~  Fv  (2.a) 

■A.  ' 

mz_  =  F„  cos  9  -  F„  sin  ©  +  mg  (2.10) 


This  force  equation  neglects  a  small  Coriolis  force  term  due  to  the  earth's 
rotation.  This  term  will  be  introduced  in  the  more  complete  analysis  of 
Section  6.9. 
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wherp  T  is  t.h»  moment  nf  -fnpT~ti«  about  th*3  ^  avi  g . 
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T7!  - - „  C.  ~Z.  41  4  —  —  1 - ±-  1 _ J_ 

r  x  um  r  Igure  x  u  xt>  Lietir  uuuo 


z^  =  -  V  sin  (©  -  a) 


(2.11) 


(2.12) 


This,  together  with  the  b,s sumption  of  &  sine.  11  £LngXe  of  sitsclc  &nd.  &  fist 
trajectory  (0  <  <  l),  allows  us  to  revise  Equation  (2.10). 


m 


v(e  -  a)  +  v(e  -  a) 


=  Fz  -  Fx9  +  mg 


(2.13) 


Equations  (2.7,  9)  and  the  relation  6  =  q  can  now  be  used  to  simplify 
Equation  (2.13). 
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zi 

0.1  /  *  *\  * 
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a  o 
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(2.1b) 


oSi  -3  *  * 

Since  the  density  factor  £ —  is  less  than  10  ,  C„  and  C  can  probably  be 

£-IU  La  Lam 

q  & 

neglected  in  comparison  with  unity  and  Equation  (2.14)  can  be  written  in  an 


even  simpler  form. 


¥■  -  ¥  -  (c:  ♦  eh  a .  c:  -  gw 


-2 


(2.15) 


a 


From  this  form  we  see  that  a  good  first  approximation  for  q  is  &.  The  remaining 
terms  in  Equation  (2.I5)  can  have  a  measurable  effect  and  so  should  be  retained. 
Equation  (2.15)  is  now  used  to  eliminate  9  =  q  and  q  from  Equation  (2.11)  with 


..  A  ,  A  A  A 

h  +  H,&  -  ft,  a  =  A,  +  G 


(2.16) 
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where  = 
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y  ffi  v  \£/ 


k  =  /  iS  radius  of  gyration  about  the  Y-axis. 

y  V  m  JT 

The  starred  aerodynamic  coefficients  contain  the  density  factor  and,-  hence, 
are  quite  small.  Products  of  starred  coefficients  are  omitted  in  comparison 
with  a  starred  coefficient*. 

Equation  (2.l6)  is  linear  but  with  variable  coefficients.  This  variation 
in  the  coefficients  can  be  removed  by  replacing  time  by  distance  as  the 
independent  variable. 


•  "  »  /  v  \ 

. .  a  +  Hja  -  l-^a  =  A1  +  G  ^  ~  J 


(2.17) 


%  U  / 


*  The  products  of  starred  coefficients  are  important  for  the  angular  motion 
of  lighter  than  air  ships,  submarines  and  topedoes-  For  these  vehicles, 


M'  can  be  negative  when  C 
1  in 


is  positive'.  See  Reference  6-10. 


a 


39 


where  H, 


■Sr  *  -2  ,  * 
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Z  D  y  '  m  m  ' 
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y  m 
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1  y  m 


G  •  -  (kf  c *  -  C*)  (gw;2) 
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The  solution  to  this  equation  for  a  statically  stable  missile  (C  <  0)  i 
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a 


a  and  0  are  determined  by  initial  conditions.  By  definition  the  influence  of 
o  o 

this  initial  disturbance  will  damp  out  if  the  missile  is  dynamically  stable 
(H1  >  0).  For  such  a  missile  the  angle  of  attack  will  eventually  reach  an 

-j- 

equilibrium  angle  aa  +  a  .  The  aerodynamic  part  of  this  angle,  a  ,  can  easily 

be  interpreted  as  satisfying  the  requirement  that  the  aerodynamic  moment 

coefficient  must  vanish  for  equilibrium,  i.e.,  C  a  +  C  =  0.  A  similar 

'  m  a  m 

a  o 
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interpretation  for  a  is  possible  when  it  is  noted  from  Equation  (2.15)  that 
-2  ,  6 

g*V  is  the  curvature  of  the  trajectory  due  to  gravity.  This  angular  velocity 

has  a  moment  whose  magnitude  is  determined  by  C  giV  .  In  addition  to  this 

\ 

there  is  a  change  in  the  curvature  due  to  decrease  in  velocity.  This  angular 

p  n 

acceleration  due  to  drag  gives  rise  to  the  quantity  k  C^giV  f~.  For  equilibrium 

y  u 

the  sum  of  these  two  terms  must  b^  matrhpd  bv  the  moment  coefficient  C  a  . 

mg 

-2  -5  a 

For  models  as  big  as  6  inches,  giV  is  always  less  than  6  x  10  J  and,  hence, 


this  gravity- induced,  steady-state  angle  of  attack  is  not  directly  rr 


mnnr’nvinV'l  , 


j  iiiv-ci.o  oaj.  _i_c  • 


The  aerodynamic  induced  angle,  a  ,  can  be  quite  large  and  is  measurable. 

& 

The  vertical  motion  of  the  missile  can  now  be  Obtained  from  Equation  (2.10). 
The  usual  change  of  independent  variable  is  made  to  absorb  the  velocity  variation, 
and  small  angles  are  assumed. 
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(2.19) 


r*  0  t 

According  to  Equation  (2.15),  ^  can  be  replaced  by  a  in  Equation  (2.19)  since 
the  remaining  terms  of  Equation  (2.15)  yield  products  of  starred  coefficients 
after  the  substitution. 
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(2.20) 


The  first  two  integrals  in  Equation  (2.20)  yield  linear  functions  plus  damped 
cosines.  The  cosines  do,  however,  differ  in  phase.  The  third  integral  can  he 

r\ 

approximated  by  s^  .  Thus,  the  vertical  motion  is  essentially  a  parabola  ih  s 
2 

plus  damped  cosines  which  arise  from  the  pitching  motion. 


Equation  (2.20)  indicates  a  measurable  influence  of  a^.  During  firing- 

table  tests  of  finned  bombs  during  World  War  II,  it  was  found  that  times  of 

6-3 

flight  were  higher  than  expected  =  In  fact,  negative  drag  coefficients  were 

required  to  fit  the  observations  by  a  theory  in  which  only  the  drag  force  is 

considered.  The  definition  of  a  shows  that  it  is  positive  for  C  negative  and, 

g  mq 

hence,  the  missile  will  fly  in  an  equilibrium  position  with  its  nose  slightly 
above  its  flight  path.  For  almost  all  missiles  the  aerodynamic  force  is  in  the 
direction  of  the  angle  of  attack  and,  therefore,  is  negative.  Thus 

*  *  a 
(C..  +  C_)  ccp  opposes  gravity  and  increases  the  flight  time.  For  short  flights 

in  a  ballistic  range  this  effect  is  too  small  to  measure.  Gravity- induced 
equilibrium  angles  do,  however,  cause  measurable  effects  for  missiles  with 
high  spin  rates.  These  effects  are  described  in  Section  7-2. 


6.3  Combined  Pitching  and  Yawing  Motion  of  a  Nonspinning  Missile 

Except  for  the  absence  of  gravity  the  analysis  for  the  yawing  motion  is 
almost  identical  to  that  for  the  pitching  motion: 


my  =  Fv  cos  ^  +  F  sin 
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(3.1) 

(3.2) 
(3-3) 
(3.U) 

(3-5) 


In  Figure  6.4,  the  positive  senses  of  t  and  0  are  indicated  and  the  correctness 
of  Equations  (3.1-3)  directly  verified.  Note  that  the  definition  of  0  is 
consistent  with  that  of  a  in  that  it  is  positive  when  the  missile’s  nose  is 


).  A 


away  from  the  positive  direction  of  the  corresponding  transverse  axis- 


M  rvvimn  1  P 

ilUi  Uiti_l._l.jr  y  Vy 


is  negative  like  C_  .  The  sense  of  is  such  that  if  has  a 


positive  direction  in  accordance  with  the  right  hand  rule.  Thus  for  pure 


yawing  motion,  '*  —  r. 


The  relation  between  r  and  3  is  Quite  similar  to  EQuation  (2.15) 
ri  _  pi  .  ,  *  .  p* \  q  „* 

— =  -  r T  iCy  ■  *  S' p  * 

P  0 


From  this,  an  equation  for  yawing  motion  similar  to  Equation  (2.17)  can  be 
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is  radius  of  gyration  about  Z-axis. 


It  follows  immediately  from  the  right-hand  rule  that  for  a  positive  C  ,  a 

np 

positive  p  induces  a  rotation  in  the  negative  3  direction.  Thus  yawing  motion 
is  statically  stable  when  C  is  positive. 

P 

With  the  exception  of  Chapter  X,  we  will  consider  throughout  this  book 

basically  symmetric  bodies.  This  symmetry  must  apply  to  mass  distribution 

as  well  as  exterior  shape,  i.e.,  1^  =  I .  With  the  present  restriction  to 

zero  spin  and  neglect  of  the  small  asymmetry  terms  (C^  =  Cy  =  =0) 

0000 

it  should  not  be  possible  to  distinguish  pitching  or  yawing  motion.  This 
slightly  metaphysical  argument  will  be  replaced  by  a  more  explicit  algebraic 
one  in  Section  6.5.  For  either  approach,  relations  between  various  coefficients 
can  be  obtained. 
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The  new  symbols  introduced  above  are  used  to  indicate  the  coefficients  of  a 
symmetric  body.  Equations  (2.17)  and  (3-7)  can  now  be  combined  by  the  simple 
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k^  =  k  =  k  is  transverse  radius  of  gyration, 
t  z  y  ^ 

The  solution  to  this  simple  equation  in  the  complex  variable  I  is 
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0,.  and  K,^  for  j  =  1.2  are  determined  by  the  initial  angles  and  angular 
jo  ju 

velocities . 

The  various  terms  in  Equation  (3.10)  can  represent  two  dimensional  vectors 
in  the  u{3  plane,  s,  then,  is  the  sum  of  these  vectors  (Figure  6.5).  A  sample 
motion,  which  is  shown  in  Figure  6.6,  is  a  damped  elliptical  motion;  the 
semi-major  axis  of  the  ellipse  is  +  K^,  its  semi-minor  axis  is  Ik^  -  K_,|, 
and  its  center  is  located  atK^e  '  Equation  (2.17)  for  the  pitching  motion 

can  be  obtained  from  the  imaginary  part  of  Equation  (3.9).  This  use  of  the 
complex  variable,  which  introduces  a  minor  elegance  in  the  problem  of  a 
nonspinning  missile,  forms  an  essential  feature  of  the  analysis  of  the  motion 
of  a  spinning  missile. 

6.4  Exact  Equations  of  Motion  for  a  Spinning  Missile 

The  motion  of  a  rigid  body  may  be  described  by  a  pair  of  vector  equations 

mT  =  F*  +  mg  (4.1) 

If  =  5?  (4.2) 

where  V* is  velocity  of  center  of  mass 
F* is  aerodynamic  force 
g* is  acceleration  due  to  gravity 
I?  is  angular  momentum 
vt  is  aerodynamic  moment. 

The  components  of  the  vector  derivatives  in  Equations  (4.1-2)  may  be  calculated 
by  the  well-known  relation 


(lA,  U2,  u  y  =  (Ux,  U2,  U3)  +  (Ux,  U2,  U3) 


(4.3) 
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where 


U.^  Ug  ,  Uj  are  the  components  of  the  vector  U  and 


(jd  is  the  angular  velocity  vector  of  the  coordinate  system 
with  respect  to  an  inertia  system. 

For  the  missile-fixed  axes  previously  defined,  the  vectors  in  the  above 
equations  may  be  given  in  component  form 
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F  =  (Fx,  Fy,  Fz) 

g"  =  (g  sin  0,  g  cos  9  sin  0,  g  cos  9  cos  0) 

4-  I 


f\ .  ),  \ 
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yzr’  V  +  XxzP  +  Xyz J* ^  (4 •  Y) 


(4.6) 

(4.9) 


£  =  (P,  1,  r) 

where  1,1,  and  I  are  products  of  inertia, 
xy  xz7  yz 

If  we  assume  that  the  plane  of  aerodynamic  symmetry  (the  XZ  plane)  is  also  a 

plane  of  mass  symmetry,  two  of  the  products  of  inertia  vanish.  (\y  -  Iyz.  =  0) 

and  H  reduces  to  the  form  applicable  to  an  airplane  (i  p  +  I  r,  I  q,  I  r  +  I  p) 

x  xz  y  z  xz 

For  an  aircraft  the  X-axis  can  be  defined  as  the  thrust  line  or  the  zero  lift 

line.  A  third  choice  could  be  that  direction  for  which  I  vanishes  (principal 

axis  of  inertia).  If  this  is  done  and  mass  rotational  symmetry  assumed  (i  =  I  ) 

y  z 

the  angular  momentum  vector  has  the  simple  component  form  (i  p,  I  q,  I  r). 

x  y  y 

If  the  aerodynamic  force  and  moment  can  be  defined  in  terms  of  u,  v,  w,  p, 
q,  and  r,  and  the  effect  of  gravity  neglected.  Equations  (4.1-2)  yield  a 
complete  set  of  six  first  order  differential  equations  in  these  six  dependent 
variables.  In  those  cases  where  gravity  has  a  measurable  effect,  relations 
between  the  Eulerian  angle  defining  the  orientation  of  the  missile-fixed  axes 
and  earth-fixed  axes*  and  the  components  of  angular  velocity  are 


In  ballistic  range  work,  the  influence  of  the  earth's  rotation  can  normally 
be  neglected  and  the  distinction  between  earth-fixed  axes  and  inertia  axes 
thereby  vanishes.  One  exception  to  this  is  the  deflection  of  subsonic  models 
due  to  Coriolis  force.  This  effect  is  discussed  in  Section  6.9. 
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(4.10) 

(4.11) 

(4.12) 


p  =  ^S-^sin9 

q  =  9  cos  0  +  cos  ©  sin  0 

r  =  -  ©  sin  0  +  cos  ©  cos  0 

Equations  (4.10-12)  now  increase  our  set  of  differential  equations  to  nine 
first  order  equations  in  nine  dependent  variables.  Finally,  if  the  location 
of  the  missile  (xg,  yg,  ZQ)>  which  is  indirectly  present  in  the  aerodynamic 
force  and  moment  through  air  density,  air  temperature  and  wind  structure,  varies 
enough  to  cause  measurable  effects,  then  the  velocities  must  be  integrated  and 
our  set  of  equations  becomes  twelfth  order. 

Fortunately,  in  ballistic  range  tests,  0  <  <  1  so  that  Equations  (4.10-12) 
can  be  linearized  and  absorbed  into  Equations  (4.1-2)  without  increase  in  order. 
The  axial  angular  velocity  component  of  the  roll  equations  can  be  solved 
independently  of  the  other  component  equations  and  the  rotational  symmetry 
assumption  reduces  the  system  to  a  quite  simple  fourth  order  system.  This 
simplicity  is  exploited  through  the  use  of  a  nonrolling  coordinate  system  and 
complex  variables. 


For  this  nonrolling  coordinate  system,  the  missile's  velocity  vector  is 
(u,  v,  w)  and  its  angular  velocity  vector  is  (p,  q,  r)  while  the  angular 
velocity  of  the  coordinate  system  is  od  =  (0,  7)  where 

q  =  q  cos  ^  -  r  sin  v  =  v  cos  ^  w  sin  ^  (4.13) 

r  =  q  sin  ^  +  r  cos  w  =  v  sin  +  w  cos  (4.14) 

t 


o 

The  angular  momentum  vector  has  the  simple  form  (I^p,  Iyq,  Iyr);  its 
transverse  components  transform  in  the  same  way  as  those  of  the  velocity 
and  angular  velocity  (Equations  (4.13-14)) 


(4.15) 


Hy  =  (lyq)  cos  $  -  (Iyr)  sin  $  =  I? 


(4.16) 


H~  =  (I  q)  sin  ^  +  (I  r)  cos  =  I  7 

v  y  y 


(4.17) 
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in  the  nonrolling  coordinates,  the  angular  momentum  for  a  missile  with  muss 


o  yimUc  v  y 


jimple  form  while  derivatives  of  vector's  cu 


UIIIJJU  I.UU  MJ 


Equation  (4.5)  are  much  simpler.  For  example,  the  angular  momentum  vector 
and  its  derivative  are 


HO  (i  p,  I?,  i?)  (I,.i  ) 

J  J 

=  (lxf>>  Iyt  +  IXP^  Iyr  -  Ixi*l)  (4.19) 


The 


simplicity  of  Equation  (4.19)  and  the  corresponding  equati 


iun  i or 


velocity  vector  makes  the  nonrolling  system  a  very  attractive  one. 


The  four  transverse  component  equations  derived  from  Equations  (4.1-2) 
reduce  to 


rv 
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(4.23) 

where  tilde  superscripts  denote  components  in  the  nonrolling 
coordinates . 


If  the  second  equation  of  each  pair  is  multiplied  by  i  and  added  to  the 
first,  the  independent  variable  changed  to  distance,  and  dimensionless  variables 
introduced,  Equations  (4.20-23)  assume  the  simple  form 


T’-  C*T  -  iyp  =  C~  +  iC~  +  (g~  +  i&J-e/V2  (4.24) 

u  i  Z  y  z 

-  c*p  -  IFF  =  (cj,  +  iC*)  k^2  (4.25) 
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where 


u 


/v  , 

V  +  1W 

V 

<3  -  i??i 

V 


7  = 


u 

V 


The  complex  variable  i  locates  the  plane  which  contains  the  velocity 
vector  and  the  missile’s  axis.  Its  magnitude  is  the  sine  of  the  angle  between 
the  velocity  vector  and  the  missile's  axis  and  this  angle  will  be  called  the 
resultant  angle  of  attack,  a^.  From  this  it  is  clear  that  y  is  cos  a_.  Since 
Pu  is  the  gyroscopic  moment,  P  is  a  measure  of  the  gyroscopic  effects.  When 
the  force  and  moment  are  specified  in  terms  of  T,  p  and  their  derivatives,  “p 
can  be  eliminated  between  Equations  (4.24-25)  and  a  second  order  equation  in 
4  results. 

6. 5  Linear  Force  and  Moment  Expansion 

The  proper  linear  terms  will  first  be  derived  in  the  more  natural  missile- 
fixed  coordinates  and  then  transformed  to  nonrolling  coordinates.  The  variable 
I  =  (v  +  iw)  V  1  is  quite  similar  to  the  angles  of  attack  and  sideslip.  Indeed, 
for  small  angles  5  =  p  +  ia.  Coefficients  of  this  variable  will  have  the 
subscript  a  for  this  reason. 

i  i 

Our  object  is  the  most  general  linear  expression  in  p,  i  ,  p  for  the 
force  and  moment  acting  on  a  symmetric  missile.  In  components  form,  Y-force 
expansion  would  be 

CY  =  el  (v/'V +  +  +  e4(ri/v)  +  e  (v/v)' 

+  e6(w/v)'  +  e7(qi/V)’  +  eg(r i/v)'  (5.1) 

A  similar  expansion  would  apply  for  Since  v/V  =  (I  +  I)/2  and  w/V  =  (I  -  I)/2i 

and  similar  relations  exist  for  the  other  quantities  in  Equation  (5.1),  the  force 
coefficient  expansion  can  be  written  in  terms  of  complex  quantities. 

CY  +  iCZ  =  V  +  V  +  V*  +  V*  +  Bl^  +  B2fr  +  B5r<  +  B4tr'  (5.2) 
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If  wo  now  consider  a  now  missile  system  which  differs  in  roll  angle,  0, 
from  our  original  system  and  identify  quantities  in  this  system  by  a  circiimf'Jox 
superscript,  it  is  clear  that 


Sy  +  iCz  =  (Cy  +  iCz)e^  =  aJ  +  A2u  +  AyV  + 


A(ift  +  (B,  *  +  "B, 


A  A|  A  f 

,h  +  E^-f  +  E,  (i  )c 

('j.'A 


A 

O  r  ri 
f  iy) 


r\ . 

We  now  select  0  to  be  an  angle  of  rotational  symmetry.  In  other  words,  when 

A. 

the  missile  is  rotated  through  the  angle  0  its  relative  appearance  to  a  fixed 
observer  is  unchanged.  Its  force  expansion  then  should  be  unaffected.  But  the 
coefficients  of  the  conjugate  variables  are  multiplied  by  e  ^ .  We,  therefore, 
have  the  important  result  that  if  a  missile  has  an  angle  of  rotational  symmetry 
of  120°  or  less,  the  coefficients  of  the  conjugate  variables  must  be  zero.  In 
other  words,  the  linear  force  expansion  of  a  symmetric  three-fin  missile  is  the 
same  as  a  body  of  revolution  and  terms  which  express  their  different  nature 
must  be  second  order  or  higher!  The  reasoning  is  identical  for  the  moment. 


For  a  rotationally  symmetric  missile 

t  ! 

CY  +  iC^  =  A.,  t  +  A2p  +  A^£  +  A,,!! 


Cm  *  iCn  ’  Cl5  +  C2“  +  C35'  +  V' 


(5-M 

(5-5) 


The  coefficients  in  Equations  (5.U-5)  are  complex  quantities.  The  character 
of  their  real  and  imaginary  parts  can  be  determined  by  a  consideration  of  their 
dependence  on  roll  for  a  missile  with  a  plane  of  mirror  symmetry.  If  the 
XY-plane  is  taken  to  be  the  plane  of  symmetry  and  the  coordinate  system  is 

.  A 

transformed  by  a  reversal  of  the  Z=axis  (Z  =  -  Z),  a  missile  wTith  mirror  symmetry 
would  be  unaffected  and  so  would  be  the  form  of  its  aerodynamic  force  and  moment 
expansion.  Under  such  a  transformation  the  velocity  and  force  vectors  reverse 
their  Z  components  but  the  angular  velocity  and  moment  vectors  reverse  their 
X  and  Y  components.  This  is  due  to  the  fact  that  the  sign  of  a  component  of 
angular  velocity  or  moment  is  directly  related  to  the  right-  or  left -handedness 
of  the  coordinate  system  through  their  definitions  as  cross  products.  For 
example,  a  positive  p  rotation  rotates  the  Y-axis  toward  the  Z-axis,  a  positive 
q  rotation  rotates  the  Z-axis  toward  the  X-axis  and  a  positive  r  rotation  rotates 
the  X-axis  toward  the  Y-axis.  A  reversal  of  the  Z-axis  then  has  the  stated 
effect  of  changing  the  sign  of  p  and  q.  The  results  of  this  transformation  can 

Oil  niw  rt  v< 

D  LuuLucu.  ii<cu  • 
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(5.10) 


Since  mirror  symmetry  requires  that  it  be  impossible  to  distinguish  a 
change  in  the  aerodynamic  force  and  moment  from  this  transformation  the 
coefficients  in  Equations  (5.9-10)  must  equal  those  in  Equations  (5. 4-5).  If 
the  coefficients  are  constants,  this  means  that  the  imaginary  parts  of  A1 ,  Av 
C2,  and  and  the  real  parts  of  A0,  Ak,  C ^  and  C*  must  be  zero.  The  presence 
of  non- zero  roll  should  affect  the  force  and  moment,  and  so  we  make  the  assumption 
that  the  coefficients  are  functions  of  p.  With  this  assumption  the  quantities 
which  previously  had  to  vanish  must  now  change  sign  with  p.  This  means  that 
they  are  odd  functions  of  p.  (A  power  series  expansion  in  p  would  contain  only 
odd  powers.)  The  remaining  real  and  imaginary  parts  do  not  change  sign  and, 
therefore,  are  even  functions  of  p.  When  p  is  zero,  the  odd  functions  are  zero. 

In  order  to  make  this  explicit,  we  will  write  these  odd  functions  as  products  of 
p  and  even  functions  of  p. 


The  b  and  d.  cocfl' icients  are  called  Magnus  coefficients  and  denote 

i  1  J 

interaction  effects  between  spin  and  cross-flow  induced  by  one  of  the  four 

complex  quantities.  Simple  examples  of  a  Magnus  effect  are  the  curving  of  a 

6-5 

baseball  and  the  hooking  or  slicing  of  a  golf  ball 

Equations  (5.11-12)  are  not  in  the  proper  form  for  use  in  the  equation.; 
of  motion.  We  need  to  convert  to  the  nonrolling  coordinates  by  multiplication 
by  exp  i 0.  This  will  transform  the  transverse  force  and  moment  variables  £ 

t  i 

and  p  properly  but  not  5  and  p  . 


.  r  -  1  (|i)  T 

(5.15) 

**7  -7 

(5.1^) 

This  follows  from  the  fact  that  derivatives  of  vectors  in  rotating  systems 
do  not  transform  in  the  same  way  as  the  vectors  themselves.  A  second 
complication,  which  follows  from  this,  is  that  a  spinning  missile  flying  at 
a  constant  angle  of  attack  in  the  nonspinning  system  (I  =  0)  has  a  non-zero 

i  » 

5  .  The  non-Magnus  coefficients  of  I  are  usually  taken  to  be  non-steady 


For  a  body  of  revolution 


damping  force  and  moment  coefficients  -  C«  and  Cw 

N  •  M. 

a  A 

flying  with  constant  £,  it  is  very  hard  to  see  why  a  damping  force  or  moment 
6-6 

should  be  present  and  it  is  difficult  to  explain  such  a  force  and  moment 
for  a  symmetric  finned  missile.  For  this  reason  we  will  rearrange  the  terms 
in  Equations  (5.11-12)  so  that  the  variables  are  i,  p,  and  their  derivatives 
in  the  nonrolling  coordinates.  The  coefficients  are  then  identified  so  that 
they  reduce  to  the  quantities  as  defined  in  Sections  6.2  and  6. 5  for  the  case 
of  zero  roll  and  small  angles. 
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One  last  task  remains  before  returning  to  the  equations  of  motion.  This 
is  the  determination  of  the  variation 'of  the  coefficients  defined  above  with 
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for  two  models  with  identical  motion  but  different  center  of  mass  locations. 

The  difference  in  c.g.  location  will  be  denoted  s  =  (x  -x  ) l  ^  where 

eg  eg  eg 

the  circumflex  superscript  identifies  quantities  associated  with  the  missile 
with  the  second  c.g.  location.  Under  our  assumption  that  corresponding 
points  have  the  same  motion,  the  angular  velocity  is  unchanged  but  the  linear 
velocity,  which  is  defined  to  be  the  velocity  of  the  center  of  gravity,  is 
affected. 


1  =  t  -  is  u 

eg 
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P  =  P 


P  =  P 


(5.17) 

(5.18) 


In  a  similar  fashion  the  force  is  unaffected  but  the  moment,  which  is  defined 
with  respect  to  the  c.g.  location,  does  change.  If  we  assume  that  the  change 

i*  P  i*  c  cmflll  pnpi  l  rrVi  ,t_r*  T  psotrps  ■+■  V»o  ■fr'i+ol  vpI  OP'ltv  V  Invariant  +-}-:£>  *o72f 

^  vuwu^ii  w  avm*  w  viiv  vJUwj.  *  ClUu  X  Ojr  j  V  y  Xuvai  IQUU^  Oiie  ICiatlULlb 

for  force  and  moment  coefficients  are  the  same  as  for  the  force  and  moment 

themselves.  This  means  that  the  resulting  c.g.  relations  are  good  for  small 
'angular  velocities . 


Cy  +  iCy  =  Cy  +  iC7 


(5.19) 
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Equations  (5. 15-18)  are  substituted  in  Equations  (5. 19-20)  and  the  scught=for 
c.g.  transformations  are  found.  -These  are  listed  in  Table  6-1. 
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ami 


This  table  reveals  a  number  of  interesting  facto.  Note  that  if  C 


M. 


;i  rp  ii .  'i  t 


a 


TT 

"4 

zero,  then  and.  must  be  non-zero  for  moot  c.g.  pooitiono. 


Thus  the  omission  of  these  coefficients  in  previous  articles  would  yield 

incomplete  expansion  under  c.g.  transformations.  C„  and  C  have  no  measuratJ< 

N.  M. 

q  q 

effect  on  missile  motion  and  are  omitted  for  this  reason*.  In  a  similar  way 


'N 


and  C„.  are  needed  for  consistent  c.g.  transformation  for  non-zero  C. 
M  D  r 


pq 


pq 


pa 


and  C„  but  have  no  measurable  effect  on  motion.  The  Magnus  coefficients 
M 


va 


involving  a  and  4  also  have  no  measurable  effect,  do  not  appear  in  c.g.  relation 
for  other  coefficients,  and,  therefore,  can  be  completely  dropped  at  this  timey*. 
The  final  form  of  the  linear  force  and  moment  expansion  in  nonrolling  coordinates 
is,  therefore. 
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Equations  (5*21-22)  contain  considerable  information  which  is  not 
immediately  obvious  due  to  the  unfamiliar  complex  number  notation.  It  should 
be  noted  that  |  =  (3  +  ia  has  components  whose  algebraic  sign  is  the  reverse  of 
the  corresponding  forces  and  moments.  This  can  be  seen  from  Figures  6.5  and  6.4. 


In  wind  tunnel  measurements  of  damping  in  pitch  derivatives  by  means  of 
forced-oscillations  the  coefficient  CM  can  have  a  measurable  effect  quire 

4 

similar  to  the  static  moment  coefficient  C 


** 


'M  * 
a 


The  various  statements  about  certain  coefficients  having  no  measurable  effect 
may  be  verified  by  retaining  these  quantities  in  the  analysis  of  the  next 
section.  It  would  then  develop  that  through  the  small  size  of  the  density 

n  " 


_r«  ft 

factor  Ts-jp  these  coefficients  would  have  little  effect  on  the  pitching 

and  yawing  motion  or  the  c.g.  transformation  of  coefficients  which  do  influence 
that  motion. 
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Thus,  the  normal  force,  -  T,  is  a  force  proportional  to  the  sine  of  the 

a 

angle  of  attack  and.  acting  in  the  same  direction  as  the  angle  of  attack  when 

CN  is  positive,  (With  the  exception  of  very  low  fineness  ratio  configuration 
a 

C„  is  always  positive.) 


The  Magnus  force  for  a  spinning  body  of  revolution  is  at  right  angles 

I/v»i 

£  I •  Its 

usual  direction  can  be  inferred  by  consideration  of  the  Magnus  force  which 

causes  a  baseball  or  golf  ball  to  "hook"  or  "slip".  (Figure  6.7)  From  the 

figure  we  can  see  that  the  circulation  induced  in  the  flow  field  about  the 

ball  bv  the  spin  and  viscosity  causes  an  asymmetric  pressure  distribution. 

This  pressure  distribution  then  gives  rise  to  a  side  force.  The  same 

theoretical  model  may  be  used  for  a  spinning  body  of  revolution.  For  this 

purpose  we  assume  that  Figure  6. 7  is  a  cross  sectional  view  of  the  flow  field. 

The  cross  flow  is  induced  by  the  angle  of  attack  and  has  a  magnitude  V  | £  |  . 

This  simple  picture  predicts  a  Magnus  force  with  direction  obtained  by  a  90° 

rotation  from  the  normal  force  against  the  spin.  Since  multiplication  by  ip 

yields  a  90°  rotation  in  the  direction  of  spin,  the  usual  sign  for  CN  is 

pa 

negative. 

The  interpretation  of  the  damping  force  coefficients  C  and  C  is 

_  I*  • 

q.  a 

somewhat  more  difficult.  This  is  not  unexpected  since  the  coefficient  C 

N_ 

U 

is  a  linear  function  of  c.g.  location.  Indeed  since  both  coefficients  have 
not  been  directly  measured  and  have  only  importance  for  their  effect  on  the 
dependence  of  their  moment  coefficients  on  c.g.  location,  we  will  interpret 
them  in  terms  of  the  c.g.  transformation  of  the  moment  coefficients  and  move 

t" r\  hnncirler  thocn  r* r\o-P-P-f  r*-f  OTvf  c? 

The  form  of  the  moment  expansion  was  obtained  by  multiplying  the  force 
coefficients  by  i  and  replacing  the  N  subscripts  by  M’s.  This  is  justified 
by  the  convention  the  orientation  of  a  moment  is  the  axis  about  which  it  tends 
to  cause  rotation  and  this  axis  is  perpendicular  to  the  corresponding  moment. 
As  each  coefficient  is  discussed  its  influence  on  stability  will  be  explored. 
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The  most  important  moment  coefficient  is  that  associated  with  the  normal 

force.  i.e.f  the  static  moment  coefficient,  C„  .  This  moment  causes  rotation 

M 

rv 

about  an  axis  at  right  angles  to  the  angle  of  attack.  Since  a  negative 
coefficient  corresponds  to  a  moment  opposing  the  angle  of  attack,  a  negative 
coefficient  is  stabilizing. 


The  Magnus  moment  coefficient,  ,  has  an  axis  of  rotation  in  the  plane 

pa 

of  the  angle  of  attack  and,  therefore,  causes  the  missile's  nose  to  rotate 

around  the  trajectory.  For  a  positive  Cw  this  motion  is  in  the  direction 

pa 

of  the  spin.  Its  influence  of  stability  depends  i 
moment  coefficient  with  the  others  and  will  be  examined  later. 


The  damping  moment  coefficients  C„  and  C„  cause  rotation  about  the 

9  * 

rvy  /■*-»! 

same  axes  as  the  corresponding  angular  velocities  u  and  |  .  Their  moments 
both  oppose  their  angular  velocities  when  they  are  negative  and,  hence, 
negative  values  of  these  coefficients  are  stabilizing. 

6 . 6  Pitching  and  Yawing  Motion  of  Symmetric  Missiles 

With  linear  expansion  of  the  force  and  moment  defined,  the  derivation 

and  solution  of  the  equations  for  pitching  and  yawing  motion  can  proceed. 

* 

Equations  (5.21-22)  are  substituted  into  Equations  (4.24-25)  and  C  ’s 
neglected  in  comparison  with  one. 
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(6.2) 


The  expression,  -  C^,  has  a  special  significance  which  we  can  derive  here. 

a 

Figure  6.8  shows  the  model  at  a  fixed  resultant  angle  of  attack,  a  ,  and  the 

X  " 

resolution  of  the  aerodynamic  force  into  axial  and  normal  components,  F  and 

X 

F^,  as  well  as  drag  and  lift  components . 


(6.5) 


D  =  -  cos  a  Fv  +  sin  a  F„ 
r  X  r  N 


L  =  sin  a  +  cos  a  F„ 
r  X  r  N 


(6.4) 


where  cos  a  =  y 
r 


sin  ar  =  | 5 |  =  S 

The  definitions  of  the  corresponding  aerodynamic  coefficients  can  be 
inserted  and  common  factors  cancelled 


CD  =  "  yCX  +  5  °N 


(6.5) 


CL  "  7CN  +  Cx 
a  a 


(6.6) 


Cx  can  be  eliminated  from  equation  (6.5)  by  the  use  of  Equation  (6.6) 


?CL  -CN  -CD 
a  a 


(6.7) 


Thus,  this  combination,  CN  -  C^,  is  essentially  the  well  known  lift 

a 

coefficient.  Equation  (6.5)  also  expresses  the  usual  relations  for  induced 

drag.  For  constant  C  and  C„  this  yields  CL  ?  =  C  +  (l/2)  C  =  0  +  (l/2) 

Q-  o  a  a 

Equations  (6.1-2)  can  be  combined  to  eliminate  ~  and  ~  . 


T"  +  (H  -  -  IP)  !’  -  (M  +  iPT)  T  -  G 

where  H  =  ?C*  -  C*  -  if  (C*  +  rC*  ) 

a  q  a 
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(By +±8%)*  '  f  .  -2  *  *  1  r  (fkr  +  W* 

U  =  L“^? — J  ‘  |_vcVCD  +  ipj  [  ~  v2  ~ 
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Equation  (6.8)  is  exact  and  has  no  linear  approximations  other  than  thorx 
which  are  implicit  in  the  definition  of  the  force  and  moment  expansions.  The 
retention  of  derivatives  of  the  lift  coefficient  and  the  Magnus  force 
coefficient  allows  these  coefficients  to  be  variables  and,  therefore,  could 
represent  a  nonlinear  force.  Since  the  moment  coefficients  were  not 
differentiated  in  the  derivation  their  derivatives  do  not  appear.  They  are 
not  necessarily  constants  and  can  represent  a  nonlinear  moment.  This  equation 
can  and  will  be  used  in  the  chapter  dealing  with  nonlinear  forces  and  moments . 

In  this  chapter,  however,  Equation  (6.8)  will  be  linearized  by  the  approximations 


7  =  1»  7  =  C_  =  C 


a 


N 


=  0 


(6.9) 
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g~  +  ig~  =  ig  (flat  trajectory) 

.  .  _p  #  *  1  _p 

G  =  !  P  -  ^kt  CM  "  V  g/y 
q  J 


(6.10) 
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For  high  spin  rates  P  can  be  as  large  as  .05.  The  imaginary  part  of  G  is  the 

same  as  that  given  in  Section  6.2  for  no  spin.  Since  it  was  found  there  that 

this  term  has  no  effect  on  the  motion  which  can  be  measured  during  the 

relatively  short  flight  in  a  ballistic  range,  it  will  be  neglected  and  G 

_2 

approximated  by  Pg/V 
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Although  Equation  (6.8)  is  a  fourth  order  system  in  the  variables  v  and 

6-6 

w  it  is  written  as  a  simple  second  order  system  in  the  complex  variable J 
This  symmetry  allows  us  to  express  the  solution  as  a  constant  plus  two  complex 
exponentials . 

id,  i0o 

*5*  ~  ,  ,  . 

I  =  §g  +  K^e  ~  +  Kge  (6.13) 


where  £ 


>  _  r  -P 

g  +  iFT_ 


K.  =  K4  e 
J  Jo 


$a  =  0^,  + 

J  J 


/,  r.  ._  +  r~  Z2  12  ~~7z 

A.J  +  1J0J  =  il/2)  -tt  +  It  -W4M  +  h  -  F  +  2lP(2T  -  H) 

The  actual  size  of  I  is  barely  measurable  but  its  influence  on  the  trajectory 

6 

can  be  measured  for  subsonic  tests  in  ballistic  ranges.  Since  M  usually  is 

much  larger  that  FT  and  is  positive,  £  is  a  small  negative  real  number.  Thus 

this  equilibrium  angle  due  to  gravity  curvature  of  trajectory  causes  the 

missile  to  point  to  the  right  of  its  flight  path.  The  lift  force  associated 

with  this  angle  causes  a  drift  to  the  right.  This  right  deflection  is 

characteristic  for  all  artillery  shell.  For  most  supersonic  tests,  5  will 

S 

be  neglected  in  any  consideration  of  Equation  (6.13)  but  its  associated  drift 
will  be  computed  when  the  swerving  motion  of  the  c.g.  is  discussed. 

The  other  terms  in  Equation  (6.13)  describe  the  missile’s  response  to 

initial  conditions.  The  amplitudes  of  the  disturbances  due  to  initial 

conditions  (K.  )  can  either  grow  or  decrease  exponentially  while  the  orientat  ion 
JO  | 

changes  at  constant  rates  (0,).  The  resulting  motion  is  known  as  a  damped 

J 

epicycle.  As  has  been  shown  previously  a  nonspinning  statically  stable 
missile  (P  =  0,  M<0)  has  a  damped  elliptical  pitching  and  yawing  motion 
when  H  is  positive,  and  an  undamped  elliptical  motion  when  it  is  negative.  The 
damping  rates  are  equal  and  the  frequencies  are  negatives  of  each  other.  If 
the  same  missile  is  spun  slowly,  the  frequencies  shift  a  little  and  the 
elliptical  motion  starts  to  precess.  This  is  shown  in  Figure  (6.9a).  At 
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higher  spin  rates  the  damping  exponents  change  as  well  as  the  frequencies  and 
more  complicated  motion  with  nodes  is  produced  as  shown  in  Figure  (6. 9h).  Noi 
that  maximum  angles  occur  when  the  modal  vectors  add  (K^  +  JC,)  and  minimum 
angles  occur  when  they  subtract  (Ik,  -  K_ I ) . 

X  d  - 

For  a  statically  unstable  missile  such  as  a  rifle  bullet  the  motion  is 
quite  different.  If  the  missile  is  not  spun,  the  modal  frequencies  are  zero, 

-  H  -  Jk  M  +  H2 


one  modal  amplitude  is  strongly  damped  = 


/ 2  where  the 


+  Jk  M  +  J 


/ 2.  This  motion 


other  amplitude  is  strongly  undamped  X0  =  I  -  H 

C, 

is  quite  unstable. 

For  a  spinning  missile  as  well  as  a  nonspinning  missile  the  predominant 
aerodynamic  effect  is  the  static  moment  coefficient.  If  the  other  aerodynamic 
coefficients  are  neglected,  the  frequencies  and  damping  exponents  become 


Xj  +  i0j  =  (1/2) 


iP  -  Jk M 


(6.14) 


The  presence  of  spin  in  the  square  root  indicates  the  possibility  of 
stabilizing  a  statically  unstable  missile  by  spin.  When  spin  is  large 
enough,  the  right  side  of  Equation  (6.14)  is  a  pure  imaginary  number.  The 
damping  exponents  become  zero,  and  the  envelope  of  the  motion  will  net  grow. 
The  necessary  condition  on  the  spin  rate  is  the  inequality 

2 

4m  -  P  <0  (6.15) 

For  a  statically  unstable  missile,  M  is  positive  and  Inequality  (6.15) 
may  be  written  in  the  form 


>  1  (6.16) 


The  ballistician  has  for  many  years  made  use  of  a  gyroscopic  stability  factor 

s  which  is  proportional  to  the  ratio  of  squared  evroscopic  spin  to  the  static 

2 

moment  coefficient  and  is  precisely  P  /4m.  Thus  s  must  be  greater  than  unity 

6 

for  a  statically  unstable  missile  to  perform  periodic  motion  and  not  quickly 
tumble.  This  is  the  condition  for  a  sleeping  top. 


Since  the  sign  of  M  is  not  essential  for  Inequality  (6.15),  we  will 
define  a  missile  to  be  gyros copically  stable  when  that  inequality  is  satisfied 
or  equivalently 


I  <1 

s 

e 


(6.17) 


Thus  a  statically  stable  missile  (M  <  0)  is  always  gyroscopically  stable. 

According  to  Equations  (6.13,  6.14)  the  pitching  and  yawing  motion  of  a 
gyroscopically  stable  missile  with  only  a  linear  static  moment  acting,  may  be 
represented  by  the  sum  of  two  complex  numbers  with  constant  magnitude  and 
rotating  with  constant  angular  velocity. 


*i  ■  1/2 


P  -  4M 


(6.18) 


We  will  call  these  numbers  modal  vectors.  If  CM  is  positive,  they  rotate  in 

“a 

the  same  direction  as  the  spin  and  the  usual  epicyclic  motion  associated  with 
the  motion  of  the  axis  of  a  top  results.  If  C  is  negative,  the  square  roots 


14 

a 

dominate  in  Equation  (6.13)  and  these  modal  vectors  rotate  in  opposite  directions. 

Gyroscopic  stability  like  static  stability  is  sufficient  for  oscillatory- 
motion  but  does  not  guarantee  that  initial  conditions  will  not  grow.  This 
requirement  of  dynamic  stability  which  reduces  to  a  need  for  positive  H  for  a 
nonspinning  missile  is  satisfied  when  the  exponential  coefficients  are  both 
negative.  The  pitching  and  yawing  motion  for  a  dynamically  stable,  gyroscopically 
stable,  but  statically  unstable  missile  is  shown  in  Figure  6.10.  For  a 
statically  stable  missile  the  frequencies  differ  in  algebraic  sign  and  the 
nodes  are  on  the  outside.  For  this  case  the  frequencies  have  the  same  sign 
and  the  nodes  are  on  the  inside. 

The  initial  orientation  and  amplitude  of  the  modal  vectors  are  functions 
of  initial  conditions.  Equation  (6.13)  and  its  derivatives  may  be  evaluated 
at  s  =  0  to  yield  and  |  in  terms  of  these  quantities. 
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These  relations  determine  I  and  1  in  terms  of  quantities  measurable  from 

o  o 

ballistic  range  tests.  In  design  studies  it  is  sometimes  important  to  compute 

rsj  rs t  t 

initial  values  of  the  modal  vectors  from  specified  values  of  and  £  .  This 
can  bo  done  by  inverting  Equations  (6.19-20). 
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(6.22) 


The  ajnplitude  of  these  disturbances  due  to  initial  conditions  will  grow 
exponentially  for  a  dynamically  unstable  missile  (x^  >  0  and/or  X^  >0). 

The  relations  for  X .  +  10  .  are  convenient  for  computing  the  frequencies 

J  J 

and  exponential  coefficients  when  the  aerodynamic  coefficients  are  known.  In 
ballistic  range  work  X  .  +  i0  .  is  known  from  the  data  analysis  and  the  aerodynamic 
coefficients  are  unknown.  The  necessary  equations  can  be  easily  derived  from 
the  fact  that  in  Equation  (6.12)  the  coefficient  of  £  is  the  negative  sum.  of 
the  roots  of  the  characteristic  equation  (x  ■  +  i0A)  and  the  coefficient  of  T 
is  their  product. 


P  -  0i  +  02 


(6.23) 
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(6.25) 

(6.26) 


b.7  Dynamic  Stability  Criteria 

One  of  the  primary  uses  of  range  data  is  the  determination  of  proper 
dynamic  stability  for  a  proposed  eonf iguration.  Tne  usual  dynamic  stability 
requirement  is  that  the  damping  exponents  are  nonpositive  throughout  the  flight 
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of  the  missile.  We  will  first  consider  a  generalized  version  of  this  stringent 
icquuL-iufin.  xiien  bne  more  moderate  conditions  or  a  realistic  designer  that 
the  angular  motion  does  not  become  excessively  large  will  be  analyzed.  As  we 
shall  see,  then,  und<  r  this  criterion  the  damping  exponents  may  become  positive 
for  short  portions  of  the  flight. 

Equations  (6.24)  and  (6.26)  may  be  solved  explicitly  for  the  damping 
exponent  s 
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The  indicated  divisions  in  Equations  (7.1=2)  may  now  be  performed 
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According  to  Equations  (7.3-4)  the  damping  exponents  for  a  nonspinning  missile 

Ls  not  zero,  one  of  the  damping  exponents  will  be  less 
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than  H/2  and  the  other  greater  than  H/2.  From  Equation  (6.18)  which  is  a  good 
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|  ^  -  jt/gj  bas  a  magnitude 
which  is  less  than  one  for  statically  stable  missiles  and  greater  than  one 
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Aitnougn  aquation  (0.10)  is  exactly  correct  for  H  =  T  =  0,  it  is  an  excellent 
approximation  when  aerodynamic  damping  is  present.  Use  of  this  equation  is 
equivalent  to  the  assumption  that  the  product  of  the  frequencies  is  much 
greater  than  the  product  of  the  damping  exponents  in  Equation  (6.25). 
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for  statically  unstable  missiles.  Thus,  the  h/2  in  the  second  term  of 
Equation  (7. 4)  dominates  the  first  term  for  a  statically  unstable  missile 
and  one  damping  exponent  must  be  positive  if  the  Magnus  and  lift  coupling 
term  T  is  zero. 

We  will  now  obtain  precise  stability  conditions  on  T  as  a  function  of 

spin.  The  usual  dynamic  stability  requirement  is  that  the  exponential 

coefficients,  X.,  be  both  negative.  These  conditions  are  somewhat  more 
J 

restrictive  than  necessary  for  a  missile  designer.  His  primary  concern 
is  that  the  initial  conditions  do  not  have  an  adverse  effect  during  the 
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trajectory  of  the  given  missile.  In  some  cases  we  would  want  initial  angles 
to  damp  to  half  their  values  in  1,000  ft;  in  other  cases  if  they  doubled  in 
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criteria  in  our  analysis,  we  will  introduce  a  level  of  damping,  X,  which  is 
set  by  the  designer  and  require  that  the  damping  exponents  be  less  than  this 
value.  The  resulting  generalized  criteria  reduce  to  the  conventional  one  when 
X  =  0.  If  the  frequencies  are  approximated  by  Equation  (6.18),  this  requirement 
of  generalized  dynamic  stability  assumes  the  form 


X  -  X  = 


(H  +  2X 


P(2T  -  H) 
F2  -  4m 


=  0 


(7-5) 


Inequalities  (7*5)  are  equivalent  to  the  following  inequalities 
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The  second  inequality  in  (7-6)  can  be  squared,  solved  for  spin  squared,  and 
reduced  to  a  simple  form  by  the  introduction  of  a  generalized  dynamic  stability 
factor,  s^ 
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The  relation  between  the  spin  and  the  various  types  of  stability  are 
s'munarised  in  Figure  6.11.  To  these  relations  must  be  added  the  single 
inequality  H  +  2\  =  0.  Three  comments  of  great  importance  may  be  derived 
from  an  inspection  of  this  figure. 

(1)  A  dynamically  stable  missile  must  be  gyroocopicaliy  stable. 

(2)  If  lies  in  the  interval  (0,2),  a  statically  unstable 
missile  may  be  stabilized  by  sufficiently  high  spin  and  a  statically  stable 
missile  is  always  dynamically  stable. 

(3)  If  lies  outside  this  interval,  a  statically  unstable  missile 
cannot  be  spin-stabilized  and  statically  stable  missile  may  be  made  dynamically 
unstable  by  a  sufficiently  high  spin. 

Experimental  values  of  the  various  coefficients  for  a  variety  of  configurations 
will  be  given  in  the  next  chapter.  From  these  data  it  can  be  shown  that  most 
of  the  ways  a  symmetric  missile  can  theoretically  become  unstable  have  been 
observed. 


The  use  of  Figure  6.11  to  determine  lower  bounds  on  spin  for  a  statically 

unstable  missile  and  upper  bounds  for  statically  stable  missiles  is  based  on 

the  assumption  that  s,  is  not  a  function  of  spin.  The  bounds  on  s  indicated 

d  g 

in  the  figure  may  then  be  converted  into  bounds  on  spin. 

Frequently  the  spin  variation  of  the  Magnus  moment  can  not  be  ignored. 
When  this  is  the  case  s^  is  a  function  of  spin  and  Figure  6.11  can  not  be 

u 

used  to  determine  spin  boundaries  for  dynamic  stability  for  a  given  missile. 
Inequality  (7-7) >  however,  is  still  valid.  The  equality  sign  in  (7.7)  yields 
an  equation  in  spin  whose  roots  are  endpoints  of  intervals  of  spin  for  which 
the  missile  is  either  dynamically  stable  or  dynamically  unstable.  If  the  spin 


in  n  m’  von  i  n+< 


a  -Hi 


In 


...VI«TU,r  (' 7  r7\  -f-Vl£S 

C^UQJ.1  OJ  ^  |  |  )  one 


be  dynamically 


stable  for  spin  in  that  interval  and  conversely. 


nPVio  non  in*  ■ 


c^viirGiiisnt  ’fciis’fc  *uxi0  exponential  coefficients  be  negative  throughout 
the  flight  is  much  stronger  than  necessary  in  a  number  of  applications.  This 
can  be  seen  by  the  following  example.  Consider  the  case  of  a  specific  shell 
whose  exponential  coefficients  are  strongly  negative  for  M  =  2.0  except  for  the 


Cc. 


Mach  number  interval  ( 0-9 >  1.1)  where  both  exponents  arc  positive.  Exact 
numerical  integration  showed  that  an  initial  maximum  angle  of  attack  of  four 
degrees  for  the  launch  Mach  number  of  two  will  decay  to  a  tenth  of  a  degree 
before  the  Mach  number  decreases  to  1.1.  The  dynamic  instability  associated 
with  the  transonic  velocities  then  will  cause  the  maximum  angle  to  grow  to 
approximately  one  degree  and  then  decrease  a  second  time  when  subsonic  stability 
is  established.  Thus  the  "dynamically  unstable"  shell  has  maintained  a  small 
angle  of  attack  over  the  entire  trajectory*. 


The  effect  of  slowly  varying  coefficients  can  be  treated  by  means  of  the 
WKB  method^  '3.  The  exact  algebra,  however,  will  make  considerable  use  of  the 
epicyclic  solution  associated  with  constant  coefficients.  We  will  assume  the 
solution  for  varying  coefficients  can  be  written  in  the  form 

a-  ^2 

t  -  K-^e  +  ^  (7.8) 

where  the  are  not  necessarily  constants  and  the  K.  are  not  necessarily 
exponential  functions.  Differentiating  Equation  (7.8)  twice,  we  have 
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+  l1 


+  (Xg  + 


(7-9) 


i  = 


xi  +  xi  '  W±)2  +  i^i  +  2xiV 


Kle 


X2  +  X2  ”  +  ^'^2  +  ^X2^2' 


i  fr 


K.e 
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(7 ,  in'! 
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where  X  . 
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Equations  (7-8-10)  are  now  substituted  in  the  homogeneous  form  of  Equation  (6.12) 
and  the  result  grouped  as  to  mode 


Precisely  this  behavior  occurred  during  the  development  of  an  important  shell. 
As  a  result  of  ballistic  range  tests  the  shell  design  was  stated  to  be 
unsatisfactory  by  the  experts.  Fortunately,  the  project  engineer  who  had  not 
studied  the  mathematical  theory  of  this  chapter  was  able  to  get  approval  for 
full  range  firings.  The  shell  gave  one  of  the  best  dispersion  patterns  ever 
observed  and  the  "experts"  learned  from  bitter  experience. 
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\($\f  -  P0*  +  M  -  X  (x,  +  H)  -  X  -  i  I  (20*  -  P)x,  +  0"  +  H0*  -  Hr]) 

l-  1  11  1  I.  1  -Lij.  |J 


=  -  {  (0g)2  -  P02  +  M  -  x2(X2  +  H)  -  4 


i  (202  -  P)\2  +  02  +  H02  -  FT  j  (K2/K1)e 


(7-11) 


where  0  =  0X  -  v2 

If  \  and  0  vary  slowly  in  a  period  of  %,  Equation  (7. 11)  can  only  be  satisfied 

J  J 

by  both  expressions  in  braces  vanishing.  Since  the  X,  terms  in  the  real  part  of 

J  t  p 

each  of  these  expressions  is  usually  much  smaller  than  (0.)  ,  the  frequency 

J 

equations  for  varying  coefficients  are  essentially  the  same  as  those  for  constant 
coefficients . 


D  rz 2 
*rl  1  7? 


(7.12) 


The  logarithmic  derivatives.  Xj,  are  somewhat  different  as  can  be  seen  from 
the  imaginary  parts  of  the  expressions  in  braces. 

I  M 

H0  -  FT  +0 
\  —  _  J  J 

J  20*  -  P 


x*-  jw/jf  -4 


(7-13) 


where  =  0^  -  | 


XJ  =  ' 


Xirk  TTTl 

gj  - 


For  constant  frequencies  and  spin- to  velocity  ratio,  X.  reduces  to  X  . 

J  j 

The  definition  of  Xt  can  be  used  to  yield  a  prediction  on  the  size  of 

j 

the  modal  amplitudes  for  slowly  varying  coefficients. 
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k./k.  = 

J  JO 


4m(0) 


P  -  4M(s) 


l/4 
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Tv 
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)  ds 


(7.14; 


Unce  the  major  concern  of  a  designer  is  the  growth  of  initial  disturbance;  , 


this  can  be  controlled  by  a  restriction  on  the  ratios  K./K 


jo 


The  actual 


value  of  these  ratios  may  be  computed  throughout  flight  and  compared  with  the 
restriction.  Clearly  strict  dynamic  stability  will  satisfy  this,  but  these 
considerations  allow  a  certain  amount  of  dynamic  instability  over  part  of  the 
flight.  Thus  a  much  weaker  type  of  stability  has  been  formulated.  This 
stability  includes  the  necessity  of  gyroscopic  stability  throughout  the  flight. 

An  interesting  application  of  these  equations  is  that  of  a  missile  with 
constant  aerodynamic  coefficients  and  spin  which  is  entering  or  leaving  an 
exponential  atmosphere* 


=  P 


-erz 


(7-15) 


where  z  is  the  altitude  in  feet. 

If  'l'(s)  is  the  angle  the  flight  path  makes  with  respect  to  the  vertical,  z  can 
be  related  to  our  independent  variable  s  by  the  equation 
s 

z  =  -  J'  l  cos  (s-^)ds^  (7.16) 

o 

o 

where  0  <  7  <90  -»  entering  the  atmosphere 

90°  <  ^  <  180°  -> leaving  the  atmosphere 


•  M  «  1, 

—  =  az  cos  ?  =  a 
P 

For  this  case  it  is  easy  to  show  that 


X  .  -  X  .  -  tt 


j  -j  4(1  -  sg; 


■r.rVi*: 

niiv. 


_  TTH 


(7.17) 


(7- 16) 


For  the  earth  a  good  value  for  a  is 


aAA  AX' 

C.C.,  (aaj  i^. 
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As  can  bo  scon  from  Equation  (7.18)  an  entering  statically  stable  missile 

(s  <  0)  is  stabilized  by  the  density  gradient  while  an  entering  statically 

unstable  missile  which  is  gyroscopically  stable  (s  >  l)  is  destabilized  by 

6 

the  density  gradient.  These  predictions  are  only  good  at  altitudes  for  which 

A 

the  changes  in  the  coefficients  are  small  in  a  period  of  0,  the  nutational 
frequency. 


6 . 8  Angular  Motion  of  a  Slightly  Unsymmetric  Missile 

Since  no  actual  missile  is  exactly  symmetric,  the  study  of  the  effect 

of  slight  asymmetries  is  very  important  for  a  missile  designer.  It  is  also 

quite  possible  that  lifting  surfaces  may  be  slightly  deflected  so  that  they 

may  induce  a  trim  angle  of  attack.  These  aerodynamic  asymmetries  which  may 

be  either  intentional  or  unintentional  have  the  effect  of  introducing  constant 

nonzero  lift  and  moment  when  the  angle  of  attack  with  respect  to  the  basic 

symmetric  missile's  axis  of  symmetry  is  zero.  This  introduces  a  constant 

C„  +  iC  into  Equation  (5.15)  and  a  constant  C  +  iC  into  Equation  (5.16'). 
i  ^  m  n  '  ' 

o  o  o  o 


The  corresponding  linear  differential  equation  for  the  pitching  and 
yawing  motion  is 

T*  +  (H  -  iP)T*  -  (M  +  iPl)T  =  lAe1^ 

Where  A  =  f§*  [^(C,  ♦  iCn  )  /p  (^)-  l)  (C  ♦  XCZ  ) 

'  L  o  o  \  \  x  /  /  a  o  J 

0  =  J'  0  ds^  and 


(8.1) 


d*  _  EJ 


For  constant  spin  rate  the  aerodynamic  asymmetry  introduces  an  exponential 
forcing  function.  A  particular  solution  to  this  inhomogeneous  equation  can 
be  obtained  by  assuming  a  solution  of  the  same  form 


This  particular  eolation  can  easily  be  found  by  direct  substitution  and  the 
result  added  to  the  epicyclic  solution  to  yield  the  general  solution  of 
Equation  (8.1) 

10,  iCL  i0. 

(''.3) 


!  =  K,e 


i0„ 

I  '  d  r 

+  IC,e  +  K,e  J 


10 


where  K,e 

3 


30 


-  iA 


-  F0 


+  M  -  U15  H  -  Pi1  j 


The  denominator  of  K-,e 

3 


10 


30 


is  usually  dominated  by  the  real  part.  This  real 
part  is  the  quadratic  expression  for  the  epicyclic  frequencies.  Therefore  it 
will  become  z-ero  when  the  spin  rate  equals  one  of  these  frequencies.  By 
definition  the  magnitude  of  the  slow  rate  is  less  than  (l^/I^) j 0 *  j  and,  hence, 
for  most  missiles  this  resonance  phenomena  occurs  when  the  spin  rate  equals  the 
fast  rate.  The  usual  value  of  I  /l  is  less  than  a  tenth  and,  thus,  the  P0  term 

•  2  x'  y 

is  much  smaller  than  (0  )  and  the  spin  rate  may  be  directly  compared  with  the 
zero-spin  pitch  frequency,  \l -  M,  to  determine  whether  concern  about  resonance 
is  reasonable. 


In  most  ballistic  range  tests  for  which  asymmetries  are  present,  the  spin 
rate  is  usually  not  constant.  If  it  is  not  near  resonance  Equation  (8,3)  is  a 
good  approximation  but  for  roll  rapidly  varying  through  resonance  it  is 
definitely  a  very  poor  approximation.  This  special  case  can  be  handled  by 
the  method  of  variation  of  parameters. 

If  varying  spin  is  assumed,  Equation  (8.1)  is  a  linear  equation  with 
varying  coefficients  with  a  specified  forcing  function.  Although  the  particular 
solution  for  this  forcing  function  is  quite  sensitive  to  roll  rates  near 
resonance,  the  complimentary  solution  for  the  homogeneous  equation  is  relatively 
insensitive.  Indeed,  for  roll,  rates  near  resonance  P  is  usually  quite  small. 

This  means  that  we  can  either  neglect  P  on  the  left  side  of  Equation  (8.1)  or 
replace  it  by  a  constant  average  value.  The  problem  is  now  to  find  the  solution 
of  an  inhomogeneous  linear  equation  with  constant  coefficients. 

The  basic  approach  of  the  method  of  variation  of  parameters  is  to  assume 

:ants  K.  ,  0  in  Equation  (6.13)  to  be  functions  of  s.  For  convenience, 
J°  <J°  .•  d 

ly 


the  quantities  A.  =  K.  e  are  introduced 
J  jo 


^7r\ 

I  v 


“  V 


(At  +  i0T )s  (A0  +  i0o)s 


+  Age 


(8.4) 


where  A .  are  complex  functions  of  s  and 


is  the  particular  solution  due  to  asymmetry  for 
zero  initial  conditions.  =  §  =0) 


Equation  8.4  is  now  differentiated. 

»  4  /l(  \  —  / \  .  j  rt  \  _ 

~i  ,  j  i  \*-0  +  ±v0)s 

E.  =  fl,  I  irf  U  e  1  +  (x  +  ±d  U  e  d  d 

JA  '  1  1  1  v  2  *2  2 


+  A,e 


»  (Ax  +  ,  (Ag  +  i02)s 


4.  A  <=> 
■  *2~ 


(8.5) 


If  the  last  two  texuis  ax-e  set  equal  to  zero 
,  (Ax  +  i0,)s  ,  (A?  +  i0?)s 

AjS  +  Age  =  0 


(8.6) 


and 
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;A  =  (X1  +  i?il)Ale 


(A-,  +  10.  )s 


+  (XQ  +  i0g)Age 


(ao  +  i0o)s 


(8.7) 


Differentiating  again  we  obtain  the  following: 

#*»"  ,t  p  (Xt  +  i0n  )s  ,  0  (A^  +  i0_)s 

t  ■  (li  +  i0l>  V  "  +  (O,  +  l?S2r  Age  c  2 


+  (a.  +  i0,  )A,  e 

'  J.  _L'  ± 


,  ,  (X1  +  i0*)a  ^  ,  (Ag  +  i0g)s 


+  (A  +  10o)Aoe 


(8.8) 


If  Equations  (8.4,  8.7-8)  are  substituted  in  Equation  (8.1)  and  it  is  recalled 
that  terms  without  A^  come  from  the  solution  of  the  homogeneous  equation  and 
must  cancel,  the  following  equation  can  be  obtained. 


(Ai  +  i0*)A|e 


,(xi  +  i0i)s  <x2 +  102>s  ..a 


+  (Ag  +  i0g)Age 


=  LAe 


(8.9) 
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Substituting  Equation  (8.6)  in  Equation  (8.9) 

-(X1  +  i0^)s  +  10 

Ai  -  — - ~r— rr  (s.w) 

X1  -  X2  +  i(0x  -  02) 

_(X2  +  i02)s  +  i 0 

a’  =  — - -r - -i-  (e.n) 

X2  -  Xx  +  i(02  -  01) 

Equations  (8.10-11)  may  now  be  integrated,  the  integrals  substituted  in 
Equation  (8.4),  and  the  general  solution  for  arbitrary  rolling  motion  0(s) 
derived 

i0,  i02  „ 

6  =  Kxe  1  +  Kge  2  +  ^  (8.12) 


where 


iA 


(X1  +  i0^)(s  -  sx)  (X2+  i02)(s  -  sx) 


-  e 


i0(s1 ) 
e  1'  ds. 


Xi  -  xp  +  i(01  -  02) 


As  was  noted  in  the  definition  of  |A ,  its  initial  values  are  zero  and,  hence, 
the  relations  between  the  initial  values  of  the  modal  vectors  and  the  angle 
and  angular  velocity  are  precisely  those  of  Equations  (6.21-22).  Since  the 

i03 

tricycle  mode  (K-,e  )  in  Equation  (8.3)  does  not  have  this  property,  the 

initial  values  of  the  modal  vectors  in  this  equation  do  not  have  the  same 
relationship  with  the  initial  angles  and  angular  velocity.  The  correct 
relations,  however,  may  be  quickly  derived  if  necessary, 

6.9  Swerving  Motion 

Once  the  pitching  and  yawing  motion  is  known  and  the  linear  aerodynamic 
force  is  determined,  the  motion  of  the  center  of  mass  may  be  computed.  A 
first  approximation  to  the  actual  trajectory  is  the  so  called  particle 
trajectory.  This  is  the  trajectory  determined  by  gravity  and  zero-lift  drag 
and  is  followed  by  a  missile  which  maintains  zero  angle  of  attack.  The  motion 
of  the  center  of  mass  of  a  missile  perpendicular  to  this  trajectory  is  called 
the  swerving  motion.  This  motion  is  quite  important  in  ballistic  range  work 
since  direct  determination  of  the  aerodynamic  force  is  possible  from  measurements 
of  this  motion. 
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For  most  ballistic  range  tests,  the  observed  trajectc 
approximated  by  a  horizontal  straight  line.  Thus  the  swerving  motion  is 
primarily  the  projection  of  the  three-dimensional  motion  on  a  plane  perpendicular 
to  this  line.  If  the  -  axis  of  the  earth  fixed  coordinate  is  directed  along 


this 


th 


:e  variation  01  tne  i  ana  Z  coorainates  of  the  center  of  motion 

e  e 


provide  us  with  the  swerving  motion.  The  differential  equation  controlling 
these  quantities  is 


m(ye  +  i5e)  =  Fy  +  ^7.  +  img  +  mac  (9.1) 

e  4'e 

a  ,  the  transverse  component  of  the  Coriolis  acceleration  induced  by  the 
c 

earth's  rotation,  can  be  specified  in  terms  of  the  latitude  of  the  range,  0  , 
and  the  azimuth  of  the  line  of  fire*,  9^, 

Sr 

w 

O  A 

a  =  — j-  a 
etc 


W 


fll 


(sin  9.  +  i  cos  9#  sin  9  ) 

Z  ZB. 


(9.2) 


where  ft  =  2jr  radians  per  day 
=  7.3  x  10^  rad/sec 

The  force  components  along  the  Yg  and  Zg  axes  are  given  by  Equations  (2.10) 
and  (3-1).  For  small  angles  these  equations  may  be  combined  with  Equations 
(2.12)  and  (3.3)  to  yield: 


F  +1F  =  F  +  iF  -  F 

Y  Z  Y  Z  rX 

e  e 


/  ir  +  ii  \ 

(r- 


(9-3) 


For  a  slightly  unsymmetric  missile  with  a  linear  aerodynamic  force,  then. 


The  BRL  Aerodynamics  Range,  for  example,  is  located  at  Aberdeen  Proving 
Ground,  Maryland  and  fires  due  East.  For  this  range  9.  =  39°  26*  and 
9a  =  90°. 
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F„  = 
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(1/2)  (pV'  )S  CD 


/q  I.  1 
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I  +  uz 


=  -  (1/2)  oV~S 
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l  a  '  '  pa-! 


(9-5) 


+  iCN  *  +  CN.  T’  •  +  1CZ  ^  I- 
q  a  o  o  J 

Equation  (4.24)  may  be  used  to  eliminate  ft  in  Equation  (9.5)  with  the  usual 
neglect  of  squared  density  terms.  The  derivatives  with  respect  to  time  may 
be  replaced  by  derivatives  with  respect  to  arclength  by  the  relation 


it  ti  /  i  \  2  /  V  £  \ 

ye  +  ize  =  (ye  +  i5e}  [t)  ’  (ye  +  ize}  W3  ) 


^ye+i2e)  (l)2  +  (ye  +  iZe)(|)  S 


(9.6) 


These  equations  may  now  be  used  to  write  Equation  (9.1)  in  a  much  more 
convenient  form 


y  +  iz^ 
e  e 


[<£  s  1  * 

u  a  v  '  pa-i 


,  *  *  ,  *  -  *  s/i 

-  CN  +  V'1’  +  (CY  +  1CZ  )e 
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Equation  (9-7)  car.  now  be  integrated. 
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where 


*i  ‘  I  f  ?  dsi 


o  o 


„  -1  8  S0 

P„\  r  r2 


If  (S)?d^ 


o  o 


t  _  r  r  ds 

x3  =  J  b 


i, .  =  /  /  e  y 


[4  =  /  /  e*y  dsi  ds; 


B  ,  are  complex  constants  and 


s  .p 

is  a  good  approximation  for  the  coefficients  of  g/V~ 

<-  Q 

and  a  . 


Note  that  if  the  spin  to  velocity  ratio  is  a  constant,  I ^  reduces  to  1^. 
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For  low  spin  the  integrals  provide  the  same  frequencies  that  are  present 
in  F, .  Thus  the  lateral  motion  of  a  symmetric  missile  with  epicyclic  angular 
motion  is  represented  by  a  parabola  plus  an  epievclic  motion.  The  relative 


amplitude  and  phase  of  the  two  motions  are  determined  by  C 


L  andCN  r'll. 
q 


A  slightly  asymmetric  missile  has  a  tricyclic  angular  motion  and  its  lateral 

motion  is  also  a  parabola  plus  a  tricyclic  motion.  The  relation  between  the 

two  three-mode  motions  is  determined  by  the  relative  magnitude  and  phase  of 

the  asymmetric  force  and  moment  terms  as  well  as  the  values  of  C  and  C  +  C 

L  N  H. 

a  q  a 

For  high  spin  the  lateral  motion  is  much  simpler.  This  is  due  to  the  fact  that 
the  magnitude  of  an  integral  of  a  circular  constant  frequency  motion  is  the 
magnitude  of  the  motion  divided  by  its  frequency.  Thus  the  integrations  in 
I I2>  and  I_  have  the  effect  of  reducing  or  eliminating  the  influence  of 
high  frequency  terms  like  0.,  and  0  and  retaining  only  the  slow  mode  0  .  Thus 


us 


the  lateral  motion  of  a  missile  with  high  spin  is  essentially  a  parabola  p 

a  spiral  with  the  slow  mode  frequency.  The  relative  magnitude  and  phase  is 

determined  by  CT  and  C„T 
L  N 

a  pa 

6.10  Aerodynamic  Jump 

The  study  of  the  center  of  mass  motion  is  quite  important  for  ballistic 
range  measurements  of  aerodynamic  force  coefficients.  An  even  more  important 
application  of  this  study  is  present  in  the  requirement  that  the  flight  path 
be  sufficiently  controlled  so  that  the  range  instrumentation  does  not  suffer 
damage.  A  simple  calculation  shows  that  the  periodic  part  of  the  motion  does 
not  have  sufficient  amplitude  to  cause  difficulty.  Indeed  this  motion  is 
seldom  greater  in  linear  extent  than  a  model  diameter. 

We  are,  therefore,  interested  in  the  average  direction  of  the  flight  path. 

The  parabolic  deflections  due  to  gravity  and  Coriolis  acceleration  are  easily 
predicted  and,  thus,  we  are  primarily  concerned  with  the  average  effect  of  the 
aerodynamic  forces. 

It  should  be  noted  that  there  are  other  causes  of  dispersion  of  models,  i.e., 
muzzle  whip,  blast  effects.  These  are  usually  small  with  the  exception  of  effect 
of  gun  blast  on  finned  missiles.  For  the  short  time  during  which  the  muzzle  gases 
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are  flowing  forward,  over  the  rear  fin  surfaces,  they  can  impart  a  large 
destabilizing  angular  momentum  to  the  model  and  a  much  smaller  linear  momentum 
This  initial  angular  velocity  can  generate  an  angular  motion  which  will  give 
rise  to  a  large  displacement.  This  displacement  is  precisely  the  quantity  we 
want  to  study  in  general.  It  is  of  great  importance  that  the  lift  of  fins  in 
this  reverse  flow  be  kept  as  small  as  possible. 


We  will  now  study  in  detail  this  average  effect  of  the  aerodynamic  force. 
This  effect  will,  of  course,  be  related  to  the  initial  angles  and  angular 
velocities  present  at  launch.  These  initial  conditions  may  be  caused  by  bore 
clearances,  sabot  separation,  or  by  blast  effects.  This  average  deflection  is 
called  aerodynamic  jump  and  may  be  defined  by  the  equation 
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In  order  to  take  the  indicated  averages  of  1-^  and  I*  we  can  solve  Equation  (8.1) 
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No  direct  measurements  of  C.r  +  C.  have  been  made  in  ballistic  ran(;e 

N  N* 

q  6c 

c  reasonable,  therefore,  to  omit  this  term  in  our  consideration  of 
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aerodynamic  jump. 
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Equation  (10.5)  shows  the  sensitivity  of  the  aerodynamic  jump  to  initial 
conditions.  Usually  i  is  small  and  a  good  approximation  to  the  symmetric 


part  of  the  jump  is 
much  on  the  roll  hi 


story. 


C.  /CM  I  •  The  asymmetric  part  depends  very 
a  /  MaJ  0 
For  constant  roll, 


$  =  — r  (10.6) 

0 

If  the  roll  is  varying  in  accordance  with  Equation  (5.5*6),  <J>  can  be  computed 

6-9 

as  a  function  of  two  parameters  " .  One  special  case  of  some  importance  is 
that  where  the  spin  producing  moment  is  large  in  comparison  with  the  spin 
damping  moment  over  the  one  or  two  periods  of  the  pitching  motion  and  the 
initial  roll  rate  is  zero.  Since  the  aerodynamic  Jump  is  an  average  over  this 
periodic  motion,  the  roll  equation  for  use  in  C>  assumes  the  very  special  form 
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For  this  constant  roll  acceleration, 
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It  is  interesting  to  note  that  the  aerodynamic  Jump  is  at  right  angles  to  JA 

A 

for  constant  spin  and  a  45  angle  for  constant  acceleration  from  zero  spin. 
Note  that  neither  (10.5)  or  (10.8)  show  any  special  effect  of  a  resonant  roll 
rate.  This  is  due  to  the  fact  that  although  the  resonance  angular  motion 
amplitude  can  be  quite  large,  its  average  influence  on  the  trajectory  is  not 
exceptionally  large.  The  increased  induced  drag  can,  however,  affect  the 
gravity  drop  along  the  Zg  -  axis. 

6.11  Spin  Lock-in  and  Catastrophic  Yaw 

Although  spin-pitch  resonance  can  yield  large  trim  angles,  the  probabil¬ 
ity  of  a  steady-state  spin  rate  exactly  that  of  tne  pitching  rate  is  quite 
small.  Moreover,  aerodynamic  damping  usually  keeps  the  trim  angle  to  be  a 
reasonably  small  multiple  of  asymmetry  angle.  By  the  introduction  of  roll 

6-n 

orientation  dependent  moments,  J,  D,  Nicolaides  has  demonstrated  the 

possibility  of  spin  pitch  lock-in  and  large  multiplying  factors  for  the  asym¬ 
metry  angle.  In  this  section,  we  will  outline  the  work  which  has  very  suc¬ 
cessfully  explained  the  rather  high  incidence  of  resonance  and  the  large 
trim  angles  associated  with  it. 

The  complex  angle  of  attack  in  missile  fixed  coordinates,  §,  can  be 
described  by  an  amplitude,  6,  and  an  orientation,  ».  G  is  the  angle  between 
the  plane  of  the  angle  of  attack  and  a  fixed  plane  on  the  missile  such  as  one 
containing  a  fin.  Nicolaides'  model  then  was  that  there  is  both  a  roll  moment 
and  a  side  moment  depending  on  angle  of  attack  and  0.  The  revised  form  of 
Equation  (5«5.2)  is 
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(11.1) 
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If  the  missile  has  an  angle  of  rotational  symmetry,  0  ,  then  C  „  is  periodic 

*0. 

witn  Deriod  fl  .  Now  the  nrpsenpp  nf  different-inf  f'apt  ft  pvr-iudp^  tip  no^cH- 
~s  ~  * ~~  ~  f 

Liiit-  of  a  plane  of  mirror  symmetry.  If  the  missile  with  no  differential 
cant  (&£  =  0),  has  a  plane  of  mirror  symmetry,  the  roll  moment  due  to  orienta¬ 
tion  must  be  zero  when  the  plane  of  the  angle  of  attack  lies  in  the  plane 
of  symmetry.  This  is  clear  from  the  fact  that  if  we  measure  0  from  a  plane 
of  mirror  symmetry  and  transform  coordinates  by  a  mirror  reflection,  the 
missile  maps  into  itself,  roll  moments  are  reversed  but  otherwise  unchanged. 


(e,  s)  =  -  c,  (-e,  &) 


(11.2) 


a  a 

A  rotationally  symmetric  missile  with  n  similar*  fins  has  a  symmetry  angle 
of  2rt/n  and  is  an  odd  periodic  function  of  2it/n.  This  means  that  it  can  usually 
he  represented  by  a  Fourier  sine  series. 
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The  simplest  form  of  for  an  n  fin  missile  is,  therefore, 

a. 


C  .  =  a  sin  n0 

1 
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(11.4) 


Maple  and  Synge  assume  that  all  forces  and  moments  are  analytic  func¬ 
tions  of  the  transverse  velocities  in  missile  fixed  coordinates  (v,w).  This 
means  that  these  forces  and  moments  must  be  at  worst  infinite  power  series  in 
I  and  I. 

But 
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In  order  for  5C„  to  be  analytic  in  v  and  w.  then,  a.  Bsin  nkfl  must  be  a  polv- 

x,  -  '  '  k  '  * 

a 

nomial  or  a  convergent  infinite  series  in  i  and  T. 
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(11.6) 


aKSsin  nk0  =  (a/  ‘  “*)  (i-)  (S'*  -  f*) 

The  last  term  in  Ea.  (11.6)  is  a  polynomial  in  4  and  I  and,  hence,  a  pol/nomic 

1  _  2  2\l/2  / 

in  v  and  w.  a  5  is  a  function  of  6  =  (v  +  w  ;  '  /V  and  would  satisfy 

^  p 

the  Maple -Synge  requirement  if  it  were  represented  by  a  power  series  in  6  . 
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and 
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Thus,  the  simplest  expression  for  an  n  fin  missile  under  the  Maple -Synge 
analyticity  assumption  is 


C£  =  b105n  "  1  sin  na 


(11.9) 


For  the  general  roll  moment,  (ll.l),  the  roll  equation  of  the  last 
chapter  (5.5*6)  assumes  the  form 
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where  K  =  ECii  go . 
a  21  l 

X  a 


0  =  0  +  0  is  orientation  of  plane  of  angle  of  attack  with  respect  to 
a  fixed  plane. 

During  conventional  resonance,  a  missile  performs  a  circular  yawing  motion 
at  constant  frequency  0  .  For  a  missile  in  circular  constant  frequency  motion 
Equation  (11.10)  predicts  a  constant  rolling  motion  when  0-0=0,  a  con¬ 
stant,  and 
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When  this  occurs,  the  roll  motion  is  "locked"  to  the  circular  pitching  mo- 

* 

tion  and  the  missile  performs  a  lunar  motion  .  This  can  occur  whenever  tnere 
*  /  \ 

exists  a  0  which  satisfies  Equation  (11.11)  and  this  equilibrium  value  of 
* 

6  is  a  stable  value. 


,  /dK  \ 

r  (sTy1 
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0  =  0 


Thus,  for  many  pitching  motions  which  are  near  circular,  lock-in  can  occur 

and  a  resonance  situation  can  exist.  The  dccurence  of  lock-in  depends  on 

the  roll  frequency  and  the  pitch  frequency  being  close  to  each  other  at  some 

point  of  the  flight  since  widely  different  frequencies  will  make  vary 

a 

rapidly  with  a  zero  average  and,  hence,  it  will  have  no  influence  on  the 
motion . 

If  we  now  consider  the  possibility  of  a  side  moment  dependent  on  orienta¬ 
tion,  a  Magnus  like  moment  could  exist  for  very  small  roll  rate.  If  this 
term  and  an  asymmetric  moment  term  are  added  to  Equation  (5.22),  we  have 
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where  CqM  is  a  function  of  0  -  $  and  6,  and  has  the  same  symmetry- 


properties  as  those  of  C .  . 

L 

a 


*  This  motion  is  so  named  since  it  is  quite  similar  to  the  moon's  rotating 
about  the  earth  in  such  a  way  that  it  presents  the  same  face  to  our  view. 
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Oar  differential  equation  (O.l)  becomes 
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For  circular  motion  uhich  is  "lock-in"  with  the  rolling  motion: 
(01  )2  -  P0'  +  M  =  0 


(11. IS) 
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where  0  =  0  s  +  0,„ 
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If  we  have  an  unfortunate  value  of  M  which  makes  the  denominator  of  tne 

5 

expression  for  quite  small,  a  very  large  amplification  factor  is  possible 
giving  rise  to  a  situation  which  is  called  "catastrophic"  yaw  by  Nicolaides. 
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FIG. 6.2  MISSILE  AXES  SHOWING  DIRECTION  AND  SENSE  OF 

FORCES  AND  MOMENTS 
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FIG.  6.3  MISSILE  IN  PITCHING  MOTION 
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FiG.  6.5  GEOMETRICAL  CONSTRUCTION  OF  PiTCHiNG  AND  YAWING  MOTION 


FIG.  6.8  RELATION  BETWEEN  LIFT  AND  NORMAL  FORCE 
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FIG. 6.9  ROLLING  STATICALLY  STABLE  SYMMETRIC  MISSILE 
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CHAPTER  VII 


RANGE  MEASUREMENT  OF  LINEAR  COEFFICIENTS 
7.1  Geometrical  Considerations 

The  basic  information  acquired  from  a  ballistic  range  test  is  the 
location  of  the  missile's  center  of  gravity  (x^,  x0,  x^)  and  the  direction 
cosines  of  its  axis  of  symmetry.  These  measurements  are  taken  in  the  range 
coordinate  system  for  which  the  x^  axis  points  downrange  along  the  intersection 
of  the  horizontal  plane  and  the  vertical  plane  containing  the  gun,  the  xg  axis 
lies  in  the  horizontal  plane  pointing  to  the  left  and  the  x7  axis  points  up. 

The  theory  of  the  preceding  chapter  made  use  of  an  earth-fixed  system,  a 
missile-fixed  system  rotating  with  the  missile,  and  a  missile-fixed  system 
with  zero  roll  rate.  The  range  system  is  essentially  the  same  as  the  earth- 
fixed  system  after  a  rotation  of  180^  about  the  x.j  axis 

*  x  =  xi  (1.1) 


ye  =  -  X£ 


(1.2) 


Z  =  -  X, 

e  3 


(1.3) 


In  order  to  relate  the  direction  cosines  to  the  angles  of  attack  and 
sideslip  in  the  non-rotating  frame  of  reference  yet  another  coordinate  system 
must  be  introduced.  (yn ,  yp,  y^)  This  coordinate  system  is  so  defined  that 
the  y.,  axis  is  along  the  missile's  axis,  the  v_  axis  is  in  the  horizontal  plane 
pointing  to  the  right  and  the  y^  axis  points  down.  Since  the  yg  axis  is 
constrained  to  lie  in  a  plane,  these  coordinates  are  called  fixed  plane 
coordinates.  They  are  approximately  the  same  as  the  missile-fixed  non-rolling 
coordinates  but  are  much  simpler  to  use.  To  verify  this  statement  we  will  first 
compute  the  roll  rate  of  these  coordinates. 

This  can  be  done  through  the  introduction  of  unit  vectors,  e^,  e^,  e^, 
along  the  fixed- plane  axes.  Since  the  y^  axis  lies  along  the  missile's  axis, 
its  components  in  the  range  axis  system  are  precisely  the  direction  cosines 
of  the  missile's  axis: 


^  =  (nx,  ng,  n,.) 
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(1.4) 


The  requirement  that  e*0  be  in  the  horizontal  plane  forces  its  third  component 
to  vanish.  Since  it  is  a  unit  vector,  perpendicular  to  e*,  and  points  to  the 
right,  its  remaining  components  are  completely  determined. 
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According  to  the  right  hand  rule, 
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The  angular  velocity  vector  (oi^,  co2,  oi^)  for  the  fixed  plane  axes  is 
defined  by  the  relations 


Equation  (1.7) 


can  be  reduced  to  the  following  relation 


(1.7) 


(1.8) 

(1.9) 
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(1.10) 


1  -  ru 


Thus  the  fixed  plane  coordinates  do  have  a  non-zero  roll  rate.  For  small 
pitch  and  yaw  angles  and  small  pitch  and  yaw  angular  velocities, 

iii 


n, 


'2 }  aa-^  y  aai>  hg , 

Therefore,  a>_,  i/V  is  a  second  order  quantity  and  for  small  amplitude  motion 
the  fixed  plane  coordinates  are  equivalent  to  the  missile  fixed  non-rolling 
system. 


n^  are  first  order  quantities  and  n^  is  essentially  unity. 
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The  components  of  the  velocity  vector  in  the  non-rolling  system  may 
be  computed  by  the  relations 


u  =  =  e^  .  V 


t-\ 

(,x. 


«*»  .  ^  ,7 

V  =  u2  =  e2  .  v 


(1.12) 


w  =  u-.  =  e,  .  V 

3  3 


(1.13) 


where  un ,  u?,  u^  are  the  components  of  V  in  the  fixed  plane  system  and 


V  is  the  velocity  vector  for  the  missile’s  center  of  gravity. 

If  (x^,  x2,  x^)  are  the  coordinates  of  the  missile’s  center  of  gravity,  its 
velocity  vector  in  range  coordinates  is  (x^,  x_.  x_) 


n2Xl  "  niX2 


nl  +  n2  ^  ^Xl^  +  ^X2^  +  (x3^ 


(1.14) 


^  nin3xi  +  n2n3x2  ~  ^nl  +  n2^  x3 

v  "  /  2  2  r~jk  .2 — 72 

V  ^  +  n_  v  +  Xg  +  x_ 

For  small  amplitude  motion  V  i  x^,  n^  =  1,  and  n2n^  can  be  neglected 


(1.15) 
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The  slope  of  the  trajectory  in  the  x^Xg  and  x-jX^  plane  can  be  easily  computed 
from  measurements  and  |  obtained  by  use  of  Equation  (l.l6).  For  large 
amplitude  the  exact  Equations  (1.14-15)  must  be  used. 


Equation  (1.16)  may  seem  obvious  to  the  reader  as  a  small  ancle 
approximation  and  he  may  wonder  about  the  introduction  of  the  fixed  plane 
coordinates  and  the  lengthy  vector  analysis.  In  answer  we  can  only  say  that 
a  careful  encineer  should  always  be  suspicious  of  an  "obvious"  result.  As 
we  shall  see  in  the  next  chapter  although  this  "obvious"  result  is  correct, 
a  similar  intuitive  feeling  that  the  fixed  plane  coordinate  oscillates  about 
the  non- rolling  coordinates  (i.e.,  (cul)ay  =  0)  is  not  true. 


7.2  Analysis  of  Motion  of  Symmetric  Missiles 

According  to  the  theory  of  Chapter  VI,  the  pitching  and  yawing  motion  of 
a  symmetric  missile  acted  on  by  linear  coefficients  is  described  by  the  equation 


i0.  id- 

rs*  rsj  rl  r  O 

4=4+  K^e  +  K^e 


where  4  1  -  (p/M)gW 
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(2.1) 


0  .  =  0  ,  +  0  .s 

Jo  ’j 

4  is  measurable  for  missiles  with  high  spin  rates  and  low  velocities. 

Gravity- induced  curvature  of  the  trajectory  usually  places  a  lower  bound  of 
500  ft/sec  on  range  tests  although  special  tests  have  been  made  at  lower 
velocities  over  a  fraction  of  the  total  range  length.  With  this  velocity 
bound,  a  statically  stable  missile  would  have  to  have  an  unusually  high  spin 
rate  so  that  4  could  be  measured  in  a  ballistic  range  test.  These  spin  rates 
are,  however,  required  for  stabilization  of  bullets  and  shell.  Since  p/M  may 
be  expressed  as  4s  /P,  we  need  only  consider  this  ratio  to  determine  the  effect 
of  model  shape,  mass  distribution,  and  spin  on  4  •  For  most  spin  stabilized 


projectiles 
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1.1  <  s  <4 
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1/30  <  <  1/10 

.25  <  &  <  .65 


(2.2) 


where  l  is  maximum  misfiile  diameter, 
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High  values  of  the  stability  factox-  occur  with  high  values  Ox  the  moment  of 
inertia  ratio  and  spin  rate  and  conversely.  Thus  the  best  maximum  estimate 
of  p/m  may  be  found  by  combining  the  low  values  of  the  above  inequalities. 


hOO  (g/v"2) 


t  a 

\c-‘D) 


According  to  Inequality  (2.3)  a  spinning  missile  with  a  diameter  10  1  ft  and 

a  velocity  of  500  ft/sec  could  have  a  £  of  six  milliradians .  This  angle  of  a 

6 

third  of  a  degree  can  be  measured  but  decreases  very  rapidly  with  increased 
velocity.  A  transonic  test  of  the  same  missile  would  produce  one-fourth  this 
value  of  C  .  A  larger  missile  would,  of  course,  have  a  proportionately  larger 

rs* 

£  .  In  any  event  |  is  quite  small  and  can  be  neglected  in  the  analysis  of  the 
angular  motion  of  supersonic  missiles.  If  it  is  necessary  to  consider  £  ,  it 

ru  g 

can  be  estimated  and  .subtracted  from  the  measured  £»s  so  that  the  motion  to 
be  fitted  is  a  simnle  epicycle 


£  -  I  =  K^e 

g  1 


(2.4) 


After  the  frequencies  of  the  epicycle  have  been  determined,  a  better  value  of 

£  may  be  computed  by  use  of  Equations  (6.6. 2-3,25) 
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-  (0X  +  02)(g*V  ) 


(2.5) 


If  the  estimate  was  quite  poor,  the  better  value  of  £  may  be  subtracted 

g 

from  the  measurements  and  the  epicycle  fit  repeated.  At  most,  one  such 
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Turning  now  to  consider  the  epicycle  of  Equation  (2.4)  we  see  that  there 
are  in  general  eight  unknown  parameters.  Four  of  these  are  the  initial 
amplitudes  and  orientations  of  the  modal  vectors  and  four  are  the  modal 
frequencies  and  damping  exponents  which  are  related  to  the  aerodynamic 


n.n^'P'Pi  ri  <=»rvh  ■hVrrmicrVi  Ulnno  +  t rtno  A 


(&  A  t-p  +v, 
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frequencies  and  damping  exponents  are  essentially  the  same  and  there  are  only 


six  unknowns.  For  this  case  the  angles  of  attack  and  sideslip  are  uncoupled 
and  each  are  described  by  damped  sine  waves  which  differ  only  in  amplitude  and 
phase.  The  analysis  of  this  special  motion  is  quite  simple  and  will  not  be 
discussed  here.  The  more  general  eight  unknown  data  reduction  will,  of  course, 
handle  this  case  and  presence  of  zero  spin  only  introduces  two  consistency 
checks  on  the  derived  frequencies  and  damping  exponents. 


The  data  reduction  is  a  combined  graphical-analytic  method  which  can  be 
divided  into  two  parts: 

(1)  No  damping  reduction  of  short  sections  of  the  observed  trajectory. 
Initial  determination  of  frequencies  and  damping  exponents  from  plots  of  phase 
and  amplitude  at  these  sections  versus  position  of  the  section. 

(2)  Refinement  of  these  initial  estimates  by  an  iterative  least 
squares  technique  which  can  be  done  by  a  digital  computer. 

Over  a  short  flight  path,  the  damping  may  be  neglected  and  the  frequencies 
determined  from  estimates  of  the  spin  and  static  moment  coefficient.  The 
approximate  slow  rate  may  then  be  used  to  reduce  the  epicyelic  motion  to  a 
much  simpler  circular  motion  by  the  transformation 

$2 

=  Kje  *  +  I^e  (2.6) 
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where  0  .  =  0  .  +0  s  +  (0 .  -  0_  )  s  and 

j  jo  ca  o  J  da 


s  is  the  coordinate  of  the  center  of  the  interval 
o 

under  analysis. 

If  the  aDDroximate  0_  is  close  to  0- .  0_  is  a  constant  and  the  motion  is 

r'da  r 2'  r 2 

i*20 

confined  to  a  circle  with  radius  and  c snt sr  at  (Figure  7  •  1 )  •  Since 

this  transformation  is  the  rotation  of  each  measured  complex  angle  of  attack 

■ 

through  the  angle  -  0o(s  -  a  ) ,  it  can  be  easily  done  by  use  of  a  compass  and 

A  i  d  0 

better  values  of  0^  determined  by  trial  and  error.  Note  that  0^  can  be  found 
from  the  angular  location  on  the  circle  of  the  complex  angle  for  s  =  s  and 

i  ° 

01  from  the  angular  location  of  the  other  transformed  angles.  The  frequencies 
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are,  however,  much  better  determined  from  a  comparison  of  phases  at  consecutive 
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amplitudes  at  the  midpoint  of  each  section  are  plotted  versus  position  on 

semi  log  paper  and  the  slopes  used  to  determine  the  \.fs. 

J 

If  these  graphical  values  are  felt  to  be  too  inaccurate,  they  may  be 

improved  by  a  least  squares  technique =  A  direct  least  squares  is  not  possible 

due  to  the  fact  that  Equation  (2.4)  is  nonlinear  in  the  0.  and  X.’s.  The 

J  J 

relations  of  the  differentials  of  these  quantities  is,  of  course,  linear 
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(2.7) 


The  differentials  of  5  are  selected  to  be  the  residuals  of  the  fit  of 
Equation  (2.4)  to  the  experimental  data 


^observed  ”  ^computed 


(2.8) 


Thus  we  have  n  equations  of  the  form  of  Equation  (2.7)  for  the  differentials  of 
the  unknown  parameters  where  n  is  the  number  of  observations.  This  set  of  n 
linear  equations  in  eight  unknowns  can  then  be  solved  by  the  usual  least  squares 
process  and  the  resulting  differentials  adding  to  their  corresponding  first 
estimates.  If  the  new  residuals  for  the  improved  set  of  parameters  are  smaller 
the  process  has  converged  and  may  be  repeated  until  the  residuals  reach  a 
minimum.  Usually  no  more  than  two  or  three  iterations  are  required. 

For  some  tests,  the  frequencies  are  observed  to  vary  between  successive 
sections  of  flight  path.  This  change  in  frequency  is  usually  due  to  varying 
P  for  a  spin-stabilized  missile  although  a  strong  Mach  number  variation  of  the 
static  moment  coefficient  sometimes  is  the  culprit.  The  frequencies  do  vary  in 
at  worst  a  linear  fashion  so  that  this  phenomena  may  be  handled  by  the 
introduction  of  quadratic  expressions  for  the  phase  angle 

01  =  +  0‘s  +  (1/2)  0';  s2 

U  U*-'  J 

and  the  use  of  a  ten  unknowns  least  squares  routine. 
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(2.9) 


When  it  is  necessary  to  use  Equation  (2.9),  the  relations  between  the 
frequencies,  damping  exponents  and  aerodynamic  coefficients  must  be  revised 
in  accordance  with  the  results  of  Section  6.7*  According  to  that  section, 
the  frequency  relations  are  unchanged  and  apply  for  the  local  values  of  P  and  M. 
Experimentally,  the  frequencies  are  best  determined  at  the  midrange  point 
(s  =0)  and  hence  values  of  P  and  M  and  their  derivatives  can  be  obtained  for 
flight  conditions  at  this  point. 


p  =  0  +0 

r  W10  *20 

(2.10) 

M  =  ^10^20  “  X1X2 

(2.11) 

p'=  +  02 

(2.12) 

M*=  0^0"  +  0^ 

(2.13) 

The  expressions  for  the  damping  exponents,  however,  are  changed. 
According  to  Equation  (6.7*13) 
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The  X  .  at  the  midrange  point  can  easily  be  computed  from  the  range 
J 

of  frequencies  and  damping  exponents  and  H  and  T  obtained  from  the 
usual  relations . 


measurements 
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Frequently,  the  0"'s  are  poorly  determined*  If  it  is  felt  that  variations 

,  J 

in  F  and/or  M  are  large  enough  to  affect  the  damping  exponents,  their  derivatives 

often  can  be  independently  measured  or  estimated  and  used  in  Equations  (2.12-1J) 

to  obtain  sufficiently  accurate  0'!*s  for  use  in  Equation  (2.14).  This  approach 

J 

naturally  depends  on  the  judgement  of  the  program  engineer. 


Although  the  angular  motion  can  be  analyzed  to  obtain  directly  the  static 
moment  coefficient,  it  does  not  yield  values  of  the  damping  and  Magnus  moments 
but  only  combinations  of  these  coefficients  with  the  drag  and  lift  coefficients, 
i.e.,  H  and  T.  The  drag  and  lift  coefficients  must,  then,  be  measured  or 
estimated  so  that  the  damping  and  Magnus  moments  may  be  obtained.  The  drag 
coefficient  may  be  easily  obtained  while  the  lift  coefficient  may  be  measured 
in  two  ways:  (l)  indirectly  from  the  static  moment  coefficients  of  two  models 


at  the  same  Mach  number  but  with  different  centers  of  gravity  (Table  6-1  and 
Equation  (6.6.7)  );  and  (2)  directly  from  the  swerving  motion  (Equation  (6.9.8)  ), 


The  analysis  of  the  swerving  motion  of  a  symmetric  missile  is  quite  simple 
because  Equation  (6.9.8)  is  linear  in  its  two  complex  unknowns,  B^  and  B1  ,  and  its 
three  real  unknowns,  CT  ,  C„  ,  and  C,_  +  C„  .  Since  the  coefficient  of 

Lct  v  \  "a 

Cw  +  C.r  is  a  single  integral  and  the  coefficients  of  CT  and  C„  are  double 
q  a  a  “pa 

integrals,  +  C„  has  a  much  smaller  effect  on  the  swerving  motion  which  has 
q  & 

not  as  yet  been  observed.  Thus  the  swerving  motion  of  a  spinning  missile  is 


linear  in  essentially  six  unknowns  &ud  the  measured  values  of  y^  and  z  may  be 
fitted  by  a  direct  use  of  the  least  squares  method.  If  the  missile  has  small 


or  zero  spin,  the  Magnus  force  has  no  measurable  effect  and  the  problem  reduces 
to  five  unknowns. 


For  a  nonspinning  missile  the  components  of  the  swerving  motion  are 
essentially  damped  sine  waves  about  a  line.  The  amplitude  of  the  wave  is 
proportional  to  the  lift  coefficient  and  could  be  obtained  graphically.  If 
spin  is  non-zero,  the  integrals  and  least  squares  solution  are  much  more 
suitable  for  calculation  by  digital  machine  than  by  an  individual  engineer. 
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7.3  Analysis  of  Motion  of  Slightly  Asymmetric  Missiles 

For  a  rolling  motion  with  a  constant  or  "nearly"  constant  roll  rate, 
the  angular  motion  of  a  slightly  asymmetric  missile  is  described  by  the 
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When  K,.  is  large  enough  to  affect  the  measured  angular  motion,  Equation  (3*1) 

3 

must  be  fitted. 


Figure  7*2  shows  a  tricyclic  angular  motion  for  which  is  20  milliradians . 
Analysis  of  this  quite  complex  notion  is  almost  impossible  without  knowledge  of 
the  roll  orientation.  In  some  range  tests  is  intentionally  induced  through 
control  surface  deflections  in  order  to  measure  control  effectiveness  but  in 
most  tests  it  appears  through  the  presence  of  slight  manufacturing  errors.  Thus 
it  is  good  policy  to  determine  the  rolling  motion  of  all  finned  missiles 
undergoing  ballistic  range  tests  so  that  the  effect  of  manufacturing  error  may 
be  eliminated  and  the  desired  static  and  damping  moment  coefficient  may  be 
obtained. 

The  roll  history  can,  then,  be  used  to  transform  the  motion  to  missile 
fixed  coordinates  by  rotation  of  each  angle  through  its  corresponding  roll 
orientation.  (Figure  7*3) 
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The  result  should  be  an  epicyclic  motion  about  the  point  K^e 
are  easy  to  determine  since  they  are  the  polar  coordinates  of  the  center  of  this 
epicyclic  motion.  The  origin  may  be  moved  to  this  point  and  the  epicyclic 
graphical  reduction  described  in  Section  7*  2  performed.  The  point  could  be 
rotated  back  through  the  angle  0  about  the  new  origin  and  a  somewhat  simpler 
epicycle  obtained  (Figure  7.4). 
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The  least  squares  iteration  may  be  performed  for  the  differentials  of 


the  unknowns . 
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Equation  (3. 3)  can  be  simplified  by  the  fact  that  the  trim  amplitude  decays 
rapidly  with  increased  spin  rate  and  at  low  spin  rates  the  sum  of  the  frequencies 
is  almost  zero.  For  moderate  Magnus  moment  coefficients,  the  difference  of  the 
damping  exponents  is  also  quite  small  for  low  spin.  If  the  graphical  analysis 
indicates  that  both  of  these  quantities  are  small,  the  differential  corrections 
can  be  reduced  to  an  eight  unknown  problem  by  requiring  them  to  be  exactly  zero. 

d|  =  (e  *  dK10  +  IK^q)  e  x 


+  (e  lSaK,0  ♦  UC.1^)  e 


+  (dKj  +  iK3d05o)  e  3 

i0n  i0o 

+  s^e  *  +  K^e  “)  dX1 


+  is^e  -  Kgf 


(3.4) 


The  assumption  of  constant  or  "almost”  constant  spin  rate  is  a  grave 
limitation  to  all  of  the  tricyclic  reductions  discussed  up  to  this  point. 

7=o 

In  1954  an  entirely  new  approach  to  the  reduction  process  was  investigated 
In  this  approach  the  actual  differential  equation  for  t  (Equation  6.8.1)  with 
T  =  0  was  programmed  for  an  analog  computer.  The  four  initial  conditions  and 
four  real  parameters  contained  in  H,  M  and  A  could  be  varied  by  changing 
potentiometer  settings.  Experimental  values  of  3  and  0  were  plotted  versus 
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s  and  a  scries  of  solution  curves  for  different  values  of  the  eight  parameters 
cuicuiatcd.  A  trial  and  error  process  could  then  be  developed  to  obtain  quite 
good  fits  of  the  data  (Figure  7-5)*  An  experienced  operator  could  obtain  good 
results  in  three  or  four  hours. 

After  our  only  experienced  operator  was  lost  due  to  the  intervention  of 
the  stork,  the  possibility  of  a  completely  mechanized  reduction  was  re-examined. 
In  Section  6.8,  a  closed  form  solution  was  derived  under  the  minor  approximation 
that  the  variation  of  P  in  the  i  of  Equation  (6.8.1)  could  be  neglected. 
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The  corresponding  equation  for  the  differential  corrections  is 
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The  initial  values  of  the  K,  and  04  s  may  he  obtained  from  estimates  of 
initial  angles  and  angular  velocities  by  use  of  Equations  (6.6.21-22).  The 
damping  exponents  can  be  initially  taken  to  be  zero  and  the  frequencies 
calculated  from  the  spin  and  estimated  static  moment  coefficient.  For  each 
guessed  value  of  E  a  sequence  of  iterations  may  be,  then,  computed  by  a 
digital  computer,  and  convergence  sought.  Usually  not  more  than  three  guesses 
of  the  complex  number  E  are  necessary  to  obtain  successful  convergence.  A 
comparison  of  these  two  methods  showed  better  fits  (one  milliradian  standard 
error  instead  of  two  milliradians)  but  essentially  the  same  values  of  M,  H  and 
A,  The  differential  corrections  method  has  the  principal  advantage  of  eliminating 
the  need  for  a  highly  trained  machine  operator. 

The  swerving  motion  can  be  handled  in  a  fashion  quite  similar  to  that  for 
a  symmetric  missile.  Two  more  unknowns  must  be  considered  in  Equation  (6.9.8), 

C  and  Cz  and  the  integrals  I  computed.  Either  Equation  (5.2)  or  Equation  (5.5) 

may  be  used  for  5  and  indicated  integrations  taken.  These  integrations  and  the 
least  squares  solution  are  extremely  laborious  for  hand  computation  but  can 
easily  be  programmed  for  a  digital  computing  machine. 

7. ^  Criteria  for  Quality  of  Results 

The  various  aerodynamic  parameters,  M,  H,  and  T,  are  related  to  the 
frequencies  and  damping  exponents  ox  the  angular  motion.  These  in  turn  may 
be  assigned  statistical  weights  through  the  least  squares  differential 
corrections  process.  If  A^'s  are  elements  of  the  inverse  matrix  obtained 
from  the  least  squares  normal  equation,  the  error  in  a  function  of  two 
variables  z  =  f(x. ,  x.)  may  be  computed  by  the  statistical  formula 
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k  is  the  number  of  unknowns  in  the  reduction. 
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Although  Equation  (4.1)  has  the  attractive  feature  of  giving  a  precise 
standard  error  for  the  aerodynamic  parameters,  it  is  a  relatively  blind 
calculation  based  on  statistical  assumptions  which  are  hard  to  verify.  In 
the  drag  analysis  an  alternate  approach  was  used  which  considered  the  size 
of  measurable  effect  of  the  drag  force  in  comparison  with  the  experimental 


accuracy.  A  similar  bus  much  more  approximate  criteria  may  be  derived  for 
the  angular  and  swerving  motion. 


Since  the  static  moment  coefficient  may  be  obtained  from  the  frequency 
of  one  modal  vector  when  the  spin  rate  is  known,  essentially  circular  motion 
with  amplitude  at  least  three  times  the  measurement  error  is  sufficient. 


When  the  spin  rate  is  not  known  both  modal  amplitudes  must  exceed  three  times 
the  measurement  error.  These  criteria  assume  that  at  least  two  revolutions 
of  the  given  modal  vector  are  observed  and  at  least  ten  individual  angles 
have  been  measured.  Good  values  of  the  static  moment  coefficient  may  need  less 
data  but  the  amplitude  requirements  must  be  made  larger.  These  amplitude 
requirements  must  be  upped  for  reasonable  determination  of  the  damping 
exponents.  With  the  amount  of  data  specified  above,  both  modal  amplitudes 

r>j 

should  exceed  five  times  the  experimental  accuracy  in  £  for  good  measurements 
of  the  damping  rates. 


The  swerving  motion  of  a  symmetric  missile  is  usually  a  spiral  motion 


with  amplitude  determined  by  C  and  phase  by  ~  C 

a  Kpa 

term  adds  another  spiral  with  amplitude  proportional  to 
average  amplitude  of  these  terms  may  be  defined  by  the  equations 
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If  the  average  radius  exceeds  three  times  the  distance  error,  the  corresponding 
force  coefficient  should  be  reasonably  well  determined. 
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5  Experimental  Results 


Since  the  194-7  firings  of  the  G.  E.  Dragonfly,  a  wide  variety  of 
configurations  have  been  tested..  These  extend  from  simple  geometrical  shapes 
of  spheres,  cones  and  cone  cylinders  to  precision  scaled  models  of  missiles 
such  as  the  Nike  Zeus,  from  20mm  bullets  to  models  of  the  Jupiter  and  Atlas 
nose  cones  and  from  the  Bolz  -Nicolaides  finned  cone-cylinders  to  five  inch 
models  of  a  delta  wing  jet  fighter.  From  this  wealth  of  experimental  work 
we  will  select  four  programs  which  not  only  illustrate  the  capabilities  of 
ballistic  range  stability  measurements  but  also  contain  results  of  intrinsic 
aerodynamic  importance. 


(l)  Static  Stability  of  Ring  Airfoils. 

Since  the  static  moment  and  normal  force  are  usually  easily  measured 
by  wind  tunnel  tests,  few  ballistic  range  programs  are  initiated  with  the  static 
stability  as  the  primary  objective.  Of  course,  when  the  Reynolds  number  and/or 
the  Mach  number  can  be  better  simulated  by  a  ballistic  range  test,  exceptions  to 
this  rule  are  made.  In  addition  to  this,  some  shapes  are  difficult  to  mount 
properly  in  a  wind  tunnel  and,  therefore,  may  require  range  tests.  A  recent 
test  of  a  series  of  ring  airfoils  was  made  for  this  reason. 


These  airfoils  had  a  simple  double  wedge  cross  section  and  the  wedge 
angle  was  varied  from  20°  to  55°  (Figure  7-6).  For  all  models  the  center  of 
gravity  was  located  at  the  midchord  point  and  the  maximum  diameter,  a,  was 
b.75  in.  The  other  dimensions  were 


9  deg. 

c  in. 

t/c 

20 

1  17 

J..  J-l 

r\  iQ 
U  •  -L  V 

25 

1. 8l 

0.23 

35 

1.70 

0.35 

These  models  were  tested  at  supersonic  Mach  numbers  up  to  3-7  and  were  found 
to  be  statically  stable  for  Mach  numbers  in  excess  of  1.7.  Below  this  Mach 
number,  the  internal  flow  choked  and  the  rings  tumbled. 


The  linearized  supersonic  theory  for  thin  ring  airfoils  is  quite 


ucv  doped  by 


similar  to  that  for  thin  two-dimensional  airfoils  and  was  d 
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H.  Mirels  in  1948.  According  to  this  theory,  the  center  of  pressure  is 
located  at  the  midchord  point  and  the  lift  coefficient,  CT  ,  is  2/f5  where 

rz — :  u<x 

P  =  VM  -  1  and  the  reference  area  is  one  half  the  wetted  area.  Thus  the 
theory  predicts  that  the  models  actually  tested  should  be  neutrally  stable. 

This  discrepancy  from  the  thin  airfoil  theory  may  be  due  to  the 
finite  thickness  of  the  models.  The  static  moment  coefficient  is,  therefore, 
plotted  versus  thickness  ratio  for  M  =  3.1.  As  can  be  seen  from  Figure  7.7, 

the  moment  coefficient  is  a  linear  function  of  t/c.  These  values  of  CM 

a 

actually  imply  a  very  large  static  margin.  If  CT  is  estimated  to  be  50  percent 

o  “ 

greater  than  its  thin  airfoil  value,  the  35  wedge  would  possess  a  static  margin 
of  almost  30  percent.  (Reference  length  was  the  chord  of  the  airfoil.) 

(2)  Dynamic  Stability  of  Cone  Cylinders . 

With  the  current  emphasis  on  the  glamorous  hypersonic  speeds,  the 
appearance  of  severe  subsonic  or  transonic  problems  is  especially  unpleasant 
for  the  aerodynamic  engineer.  Unfortunately,  a  manned  re-entry  vehicle  must 
fly  at  these  pedestrian  velocities  and  the  blunt  ICBM  nose  cones  have  subsonic 
impact  velocities.  Ballistic  range  tests  of  a  variety  of  low  fineness  ratio 
bodies  of  revolution  similar  to  shapes  used  for  re-entry  have  revealed  an 
important  subsonic  dynamic  stability  problem.  The  linear  damping  moment 
coefficients  have  been  positive  at  high  subsonic  velocities  for  certain  center 
of  gravity  locations. 

This  interesting  result  was  first  obtained  for  cone  cylinders  with 

7_e 

fineness  ratio  of  five'  (Figure  7*°)  The  actual  dependence  of  the  damping 

moment  coefficients  on  c.g.  location  was  given  in  Table  6-1. 
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According  to  Equation  (5.1),  the  sum  of  the  damping  moment  coefficients  is  a 
quadratic  function  of  c.g.  location.  This  equation  may  be  simplified  for 
ballistic  range  use  by  the  introduction  of  a  modified  damping  moment  coefficient 


A  A 

CM  +  CM. 
q  a. 


=  ^ 


\  +  \  ^  °Ma 


b-i) 


This  quantity  is  directly  computable  from  the  directly  measured  quantities 

C  +  C  and  Cw  and  the  distance  to  the  c.g.  from  some  reference  c.g. 
MM.  M 

q  a  CL 

location.  For  this  modified  damping  coefficient  sum  the  c.g.  dependence  is 
linear  and  the  damping  force  coefficient  sum  at  the  reference  point  is  the  slope 

of  a  line  fitting  experimental  points  in  a 
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Three  c.g.  locations  were  used  in  the  tests  of  the  cone  cylinder  model. 
These  were  the  centroid,  0.40  caliber  forward  of  the  centroid  and  0.45  caliber 
rear  of  the  centroid.  The  centroid  was  taken  as  the  reference  point  and 
modified  damping  moment  coefficients  were  plotted  against  c.g.  location  for 
two  Mach  numbers  (Figure  7*9) • 

At  M  =  1.26  the  scatter  of  data  is  rather  poor  but  the  subsonic  points 

have  very  little  scatter  and  are  definitely  linear  with  c.g.  location.  Since 

the  modified  damping  moment  coefficients  at  the  centroid  reduce  to  the  unmodified 

coefficients,  we  see  that  CL,  +  C..  is  definitely  positive  at  the  centroid  for 

Mq  Md 

M  =  .8.  Indeed  the  further  information  that  Cw  at  the  rear  c.g.  is  5.7  can  be 

Ma 

used  to  show  this  sum  to  be  positive  at  the  rear  c.g.  as  well. 

These  values  of  CL.  +  CL,  are  sufficient  to  make  H  negative  and  the 
Mq  Md 

missile  with  the  middle  or  rear  c.g.  location  must,  therefore,  be  dynamically 
unstable.  To  add  insult  to  injury,  it  wets  found  that  spinning  cone  cylinders 
at  M  =  .8  were  dynamically  unstable  for  all  c.g.  locations.  This  state  of 
affairs  is  caused  by  the  Magnus  moment.  In  Figure  7*10,  the  Magnus  moment 
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coefficient  is  plotted  versus  c.g.  location.  As  can  be  seen  from  this  figure, 
the  center  of  pressure  of  the  Magnus  moment  changes  radically  with  Mach  number 
while  the  Magnus  force  coefficient  which  is  the  slope  of  the  lines  in  this  plot 
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and,  t-lms ,  the  dynamic  s !  ability  factor,  s^,  will  be  negative  when  H  is  positive. 
According  to  the  stability  analysis  of  the  preceding  chapter,  dynamic  instability 
will  therefore  occur  for  any  spin  rate.  This  simple  conf iguration  displays  two 
of  the  three  ways  in  which  a  staticallv  unstable  missile  can  not  be  made 


dynamically  stable 
third  way  (s ,  >2). 


by  spin.  The  fourth  illustrative  program  will  include  the 


(3)  Magnus  Characteristics  of  Finned  Missiles. 

A  symmetric  finned  missile  with  a  reasonable  stability  margin  does 
not  require  spin  for  stability  and,  therefore,  a  Jmowledge  of  its  Magnus 
characteristics  is  unnecessary.  Unfortunately,  perfectly  symmetric  missiles 
are  impossible  to  manufacture  and  the  aerodynamic  asymmetries  possible  with 
usual  manufacturing  tolerances  are  sufficient  to  cause  undesirable  dispersion 
when  the  spin  rate  is  zero.  Usually  this  spin  rate  required  to  reduce  dispersion 
due  to  asymmetries  is  not  large  enough  to  encounter  Magnus  instability  although 
care  is  required  to  avoid  resonance.  It  is  important  to  know  what  are  the  spin 
bounds  for  Magnus  instability  if  any  so  that  too  high  a  spin  rate  is  not  induced. 


One  case  where  high  spin  rates  may  not  be  avoided  is  that  of  a  finned 
rocket  with  short  burning  time  with  respect  to  the  pitching  period.  In  order 
to  avoid  dispersion  due  to  rocket  jet  misalignment,  several  revolutions  during 
burning  should  he  induced  and  quite  large  spin  rates  are  possible.  A  knowledge 
of  the  Magnus  characteristics  of  a  proposed  rocket  shape  of  this  type  would  be 
essential  to  prove  the  feasibility  of  its  satisfactory  design. 

A  brief  examination  of  the  theory  of  Section  6.7  shows  that  Magnus 
instability  is  impossible  if  the  dynamic  stability  factor  is  between  zero  and 
two.  This  bright  possibility  of  no  spin  upper  hound  is  considerably  dimmed  by 
the  fact  that  for  every  finned  missile  whose  Magnus  characteristics  have 


been  measured  lias  been  found  to  be  outside  this  interval, 
on  spin  for  dynamic  stability  is  always  present. 


Tnus,  an  upper  bound 
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A  study  of  the  properties  of  two  finned  missiles  is  sufficient  to 
obtain  characteristic  values  of  the  coefficients  and  a  feeling  for  the  effect 


of  these  coefficients  on  the  spin  bounds.  These  missiles  —  the  Basic 

f-J  ^  * **7  Q 

Firmer  ’  and  the  127/60nsn  Anti-Aircraft  Missile  —  are  shown  in 


Figures  7*11  and  7*12  and  some  of  their  physical  properties 


are  listed  below. 
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The  Basic  Firmer  has  a  ten  caliber  lone  cone-cvlinder  bodv  with  20° 

w  v  t/  — - — 

total  cone  angle.  Its  four  simple  wedge  fins  are  square  in  planeform  with 

chord  equal  to  the  body  diameter,  d,  and  have  a  thickness  ratio  of  .06.  This 

missile  has  been  tested  with  various  c.g.  locations  and  at  various  Mach 

7-9 

numbers  both  with  and  without  fins  . 

The  127/60mm  missile  was  the  longest  firmed  missile  tested  in  a 
** 

ballistic  range.  Its  127mm  span  fins  are  mounted  at  the  end  of  a  59  inch  long 
body  whose  maximum  body  diameter  is  60mm.  This  missile  is  fired  from  a 
smoothbore  127mm  gun  by  use  of  a' center  ring  sabot. 


The  various  aerodynamic  coefficients  for  these  missiles  are  plotted 
versus  Mach  number  in  Figures  7.15-15*  Since  several  c.g.  locations  were  used 
in  the  flight  tests  of  the  Basic  Firmer,  the  normal  force  coefficient  could 
he  measured  both  directly'  from  the  swerving  motion  as  well  as  indirectly  from 
the  static  moment  coefficient  for  different  c.g.  locations.  for  the 

127 / 60mm  missile  could  only  be  obtained  from  the  swerving  motion.  The  normal 
force  is  essentially  the  same  for  both  missiles.  The  greater  effective  area 
of  the  Basic  Firmer  fins  does  provide  slightly  more  lift.  These  fins  actually 
deliver  four  times  the  lift  of  the  body  as  can  be  seen  from  the  body  only  curve. 


* 

From  base  in  body  diameters, 

**  Recently  the  45-5  inch  long  five  inch  high  altitude  gun  probe  was  tested  in 
the  Transonic  Range  and  established  a  new  Xsiig'wii  r0corcl  7-15 


115 


The  damping  moment  and  Magnus  moment  data  for  the  12  y/ 60  mm  missiJe 
are  much  poorer  than  those  of  the  Basic  Finner.  'The  12'7/60mm  fins  were  usually 
at  least  slightly  damaged  by  sabot  fragments.  This  damage  introduced 


aerodynamic  asymmetry. 


Since  the  spin  vels  noi  luc&siirsd.  during  "this 


rimr/Torri  . 


a  tricycle  reduction  could  not  be  employed.  The  epicycle  reductions  which  were 
used  had  average  to  poor  fits  and  corresponding  poorer  determination  of  the 


damping  and  Magnus  coefficients. 


The  much  larger  Cw  u  C„  for  the  127/60mm  is  due  to  the  iact  that 

m  M. 

q  a 

its  fins  have  twice  the  lever  arm  that  those  of  the  Basic  Finner  have  and, 
therefore,  the  same  angular  velocity  will  induce  twice  as  large  a  cross  velocity. 


Thi 


o  xcau 


us  to  expect  damping  moments  foux-  times  lax'ger  and  this 


expectation  is  surprisingly  well  verified  by  Figure  7.14. 


The  Magnus  moment  coefficient  behavior  for  these  two  missiles  is 


completely  unexpected.  The  moments  are  even  opposite  in  algebraic  sign!  A 
simple  physical  model  can  be  constructed  to  explain  the  positive  Magnus  moment 
but  not  the  negative  moment. 


The  127/ 60ram  missile  increases  its  spin  rate  throughout  its  flight 
tlirough  the  range.  Thus,  each  fin  has  a  associated  force  in  the  direction  of 
rotation.  Normally,  these  forces  cancel  each  other  and  provide  no  yaw  moment 
but  only  a  roll  moment.  The  fins,  however,  are  partly  in  the  boundary  layer 
separated  flow  of  the  body  and  it  is  not  unreasonable  to  expect  the  lee  side 
fin  of  a  missile  at  angle  of  attack  to  have  less  effective  area  and  hence  less 
lift  force  than  that  of  the  windward  side  fin.  This  unbalanced  force  will  cause 
the  missile's  nose  to  rotate  in  a  direction  of  the  spin  and,  thereby,  generate  a 
positive  Magnus  moment  coefficient. 


Since  the  Basic  Finner  does  not  have  its  fins  in  a  separs 


has  its  spin  rate  at  the  steady  state  rate  for  most  of  its  flight,  it  is  not 
surprising  that  its  Magnus  moment  is  quite  different.  An  explanation  for  its 
negative  coefficient  is  not  possible  at  this  time. 


These  radically  different  Magnus  moments  provide  essentially  different 
dynamic  stability  factors,  for  the  Basic  Finner  is  negative  while  for  the 
127/60n®i  it  is  positive  and  greater  than  two.  The  dynamic  stability  factors 
and  their  associated  upper  bounds  for  spin  are  plotted  in  Figure  7.16. 
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Immediately  after  World  War  II,  the  linearized  theory  of  projectile 
motion  was  reasonably  well-understood  but  few  actual  values  of  the  damping 
and  Magnus  coefficients  were  available.  A  joint  Army- Navy  program  was 
initiated  to  fill  this  void.  Since  at  that  time  these  measurements  could 
only  be  made  by  a  ballistic  range  and  the  only  available  range  was  the  BRL 
Aerodynamics  Range,  the  Army  agreed  to  construct,  measure,  and  launch  the 
models  while  the  Navy  provided  its  data  analysis  facilities  at  MIT.  The 
bodies  tested  were  secant  ogive  cylinders  with  two  caliber  heads  and  three 
different  cylinder  lengths.  The  radius  of  the  secant  ogive  was  twice  that  of 
a  tangent  ogive  and  was  so  selected  for  minimum  drag.  The  total  lengths  of 
these  models  were  five,  seven,  and  nine  diameters.  Three  center  of  gravity 
locations  were  used  for  each  body  length  and  tests  were  made  of  each 
configuration  at  three  different  supersonic  Mach  numbers.  Thus,  twenty-seven 
sets  of  linear  coefficients  were  obtained  from  this  ambitious  program^-1. 

From  this  mass  of  data  we  will  extract  several  interesting  results. 

In  Figure  7*17  the  Magnus  moment  coefficient  for  the  various  total 
lengths  and  center  of  gravity  locations  is  plotted  for  a  Mach  number  of  1.8. 

The  scatter  of  the  data  is  clearly  small  enough  to  allow  a  good  determination 
of  the  Magnus  force  coefficient  and  its  center  of  pressure. 

A  rather  simple  model  for  the  Magnus  force  on  a  cylinder  has  been 
obtained  by  J.  C.  Martin  ~  .  According  to  this  theory  the  boundary  layer 

thickness  on  a  cylinder  at  angle  of  attack  grows  as  one  moves  from  the 
windward  side  to  the  lee  side.  If  the  cylinder  is  spinning,  the  spin  will 
cause  the  plane  of  symmetry  of  this  boundary  layer  configurations  to  rotate 
to  a  position  out  of  the  plane  of  the  angle  of  attack.  The  pressure  distribution 
for  the  resulting  effective  shape  can,  then,  be  found  by  the  usual  slender  body 
approximations  and  integrated  to  yield  a  force  component  normal  to  the  plane 
of  the  angle  of  attack.  For  a  laminar  boundary  layer  Martin  shows  that  the 
Magnus  fores  coefficient  has  the  form 
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A  *  , 

C  =  AL5  ( ;-3) 

pa 

A 

where  L  is  the  cylinder  length  in  calibers 
* 

5  is  the  boundary  layer  displacement  thickness  in  diameters 
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The  actual  boundary  layer,  however,  was  partly  turbulent  and  the  actual  chape 
was  not  a  cylinder.  The  displacement  for  either  a  laminar  or  turbulent  boundary 
layer  can  usually  be  expressed  as  a  power  of  the  length  whose  exponent  is  l/2 
for  laminar  flow  and  approximately  l/5  for  turbulent  flow.  If  we  assume  a  nth 


power  variation,  the  c.p.  for  the  Magnus  force 
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<  9/14  for  l/2  >  n  >  l/p. 


(5*6) 


The  presence  of  the  nose  may  be  incorporated  in  the  theory  by 

replacing  the  total  model  length.  L,  by  an  effective  cylinder  length  L  -  L 

n 

where  is  the  same  for  all  models  at  the  same  Mach  number.  Equation  (5.6) 
assumes  the  form 


c-p-M-(fHr)  a-  L„>  (5.7) 

The  possible  range  of  the  coefficient  in  Equation  4.7  is  between  . 600  and  .643 
and,  hence,  this  C,P.M  dependence  on  type  of  boundary  layer  is  quite  small. 

The  C.P.M  is  plotted  against  length  in  Figure  7.18  and  lines  fitted  with  slope 
.64.  The  accuracy  of  this  predicted  dependence  of  C.P.  on  length  is  much  better 
than  our  approximations  deserve. 


In  Figure  7-19,  the  modified  damping  coefficients  are  plotted  versus  c.g. 
location.  This  figure  also  shows  that  the  force  coefficient  at  the  centroid 
and  damping  moment  coefficient  of  the  centroid  are  reasonably  well  determined. 
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For  bodies  of 


7-11 

Slender  body  values  of  these  coefficients  have  been  computed 
revolution  whose  reference  length  is  base  diameter  and  reference  area  is  the 
base  area  then  these  slender  body  values  reduce  to  simple  expressions 
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where  x  is  the  distance  from  the  base  to  the  center  of  gravity 
in  diameters 

A 

V  is  body  volume  in  cubed  diameters 

and  xc  is  distance  from  the  base  to  the  centroid  in  diameters. 

According  to  these  relations  CM  +  C  at  the  centroid  is  a  complicated  function 

1  Q 

of  added  cylinder  length,  but  CM  +  CM  at  the  centroid  has  an  exceptionally 

q.  a 

simple  form. 


a.  2 


(CM  +  CM  )  =  -  2  xr 


iU*  v- 

q  a 


(5.12) 


a2 

In  Figure  (7.20),  the  ratio  of  (C„  +  C.,  )  to  x  is  plotted  for  the  three  body 

C  C  ” 

a 


lengths  and  three  Mach  numbers  and  compared  with  the  predicted  value  of  -2.  The 
predicted  dependence  on  length  is  good  but  the  coefficient  should  be  decreased 
from  -2  to  -3.5- 

A  study  of  the  dynamic  stability  of  the  nine  caliber  long  models  shows 
that  the  third  possible  form  of  linear  instability  is  actually  physically 
realizable  (s^  >  2).  At  a  Mach  number  of  1.8  the  forward  center  of  gravity 
models  have  a  of  2.5.  Thus,  all  three  combinations  of  aerodynamic 

coefficients  for  which  a  body  of  revolution  can  not  be  made  dynamically  stable 
by  spin  (s^  <0,  >  2,  H  >  0)  have  been  observed. 
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are  twice  those  of  this  book  due  to  the  use  of  i/2  instead  of  l  in  reducing  the 
angular  velocities  to  dimensionless  quantities . 
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CHAPTER  VIII 


RANGE  MEASUREMENT  OF  CUBIC  COEFFICIENTS 


8 .  1  Introduction 

In  Section  50>  we  considered  the  effect  of  a  drag  coefficient  which  was 
a  linear  function  of  the  squared  angle  of  attack.  In  that  case,  it  was  found 
that  the  usual  drag  analysis  yielded  a  "range"  value  of  the  drag  coefficient 
which  corresponded  to  an  "average"  5^,  i.e.  “3^. 

The  analysis  of  the  pitching  and  swerving  motion  up  to  this  point  has  been 
limited  to  linear  forces  and  moments.  Thus,  the  corresponding  coefficients  were 
assumed  to  be  independent  of  angle.  The  large  angle  of  attack  program  involving 
the  20mm  bullet  which  was  mentioned  in  Section  5*3  included  complete  reductions 
of  the  pitching  and  swerving  motion.  Although  angles  up  to  twenty-five  degrees 
were  used  and  these  angles  were  outside  the  linearity  range  for  the  forces  and 
moments,  reasonably  good  fits  of  the  observed  motion  were  obtained.  The 
derived  values  of  the  various  coefficients  did  vary  with  the  amplitude  of  the 
pitching  motion.  Thus,  it  should  be  possible  to  correlate  these  "range"  values 
of  the  coefficients  with  the  proper  effective  angles. 

In  this  chapter,  the  angular  motion  of  missiles  acted  on  by  a  cubic  static 
moment  will  be  analyzed  in  detail.  This  will  be  done  first  by  neglecting  the 

.  i 

influence  of  geometric  nonlinearities  (7  =  0,  7  =  l).  Then  moderately  large 

geometric  angles,  a  cubic  Magnus  moment  and  linear  lift  force,  and  damping 
moments  will  be  inserted  and  relations  between  the  "range"  values  and  the  proper 
average  angles  derived.  Finally,  a  similar  analysis  will  be  developed  for 
swerving  motion  with  cubic  lift  and  Magnus  forces.  The  utility  of  these  relations 
will  be  demonstrated  by  comparison  with  certain  exact  solutions,  and  application 
to  actual  range  flight  tests. 

8.2  Cubic  Static  Moment;  Quasi-Linear  Solution 

The  simplest  possible  nonlinearity  in  the  aerodynamic  moment  is  a  cubic 
static  moment.  The  usual  symmetry  consideration  requires  that  it  assume  the 
following  form 

C  +  iC  =  -  iCw  i 
m  n  M_ 
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attack  in  an  expansion  of  the  static  moment,  not  the  derivative'  of  i.be  static 
moment.  Thin  distinction  vanishes  for  a  linear  moment  but  in  iiruxjrfani 
otherwise. 

The  as  sump  i  ion  that,  only  a  cubic  static  moment  in  act-inn  reduces 
Equation  (  . t  .12)  to  a  very  simple  form 
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If  vanishes,  the  solution  to  Equation  (2.2)  would  be 

i0n 
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(2.i) 


Our  objective  is  to  find  an  approximate  solution  to  the  nonlinear  Equation  (2.  ) 
which  has  a  form  similar  to  that  of  Equation  (2.3).  This  solution  wo\ELd  then 
supply  relations  between  the  observed  frequencies,  the  observed  damping 
exponents,  the  static  moment  coefficients,  cq,  c^  and  the  modal  amplitudes, 

Kj_,  Kg. 

The  technique  used  will  be  an  extension  of  the  Kryioff-Bogoliubof f 

p  _  o 

equivalent  linearisation  method  ’  .  This  method  in  a  generalization  jf  the 
method  used  in  Section  o-7  for  variable  coefficients  and  can  be  used  for 
nonlinear  damping  as  well.  Somewhat  more  accurate  analyses  of  the  cubic 
static  moment  motion  are  possible  from  a  consideration  of  the  exact  elliptic 
integral  solutions  but  these  analyses  lead  to  essentially  the  same  relations 

8-3  4 

when  applied  to  ballistic  range  data  ’  .  We  will  make  use  of  the  exact 
solution  to  determine  the  error  implicit  in  our  approximations. 
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Wo  assume  that  the  K.  and  0.,s  of  Equation  (2.3)  urn  function:;  of  r,  and 
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differentiate  that,  equation  twice  • 
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Equations  (2.3-2. 5)  may  now  be  substituted  in  Equation  (2.2)  and  the  result 
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For  a  linear  static  moment,  the  left  and  right  sides  of  Equation  (2.6) 

are  quite  similar  in  form  with  the  exception  of  the  presence  of  factor  K7^e~  . 

|  *  -1- 

If  0  ,  and  K ,  vary  slowly  in  a  period  of  m,  the  linearized  Equation  (2.6)  could 
J  J 

only  be  satisfied  by  requiring  that  both  sides  vanish.  Another  way  of  saying 
this  would  be  that  the  periodic  term  on  the  right  of  the  linearized  Equation  (2.6) 
has  no  effect  on  the  almost  constant  terms  on  the  left  side.  Indeed,  the  average 
of  this  term  over  a  period  of  %  is  essentially  zero.  The  influence  of  the 
fluctuating  nonlinear  term  will,  therefore,  be  obtained  by  computing  its  average 
over  a  period  of  $. 
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Equation  (2.6)  may  now  be  averaged  with  the  aid  of  Equation  (2.8),  the  small 
damping  term  in  the  real  part  may  be  neglected  in  comparison  with  the  squared 
frequency,  and  a  pair  of  real  equations  formed  from  the  real  and  imaginary  pare.: 
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Similar  relations  for  the  other  mode  may  be  obtained  by  multiplying  Equation  (2.cj 
by  (K^/Kg)e  '  and  averaging  the  result. 
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According  to  Equations  (2.10)  and  (2.1?),  damping  of  the  modal  amplitud'  can 
be  caused  by  frequency  changes.  These  changes  may  be  caused  fy  Mach  number 
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exiting  or  entering  missile  or  changes  in  the  spin  to  velocity  ratio,  P.  The¬ 
re  lat  ions  for  these  variations  are  exactly  the  same  as  those  derived  in  Section 
6.7.  The  nonlinear  frequency  equations.  Equations  (2.9)  and  (2.il),  show  a 
fourth  cause  of  a  frequency  change.  This  would  be  caused  by  a  change  in  modal 
amplitude  and  would  entail  a  coupling  between  the  frequency  and  damping  equation 
For  constant  coefficients  and  no  aerodynamic  damping,  the  modal  amplitudes  and 
frequencies  are  constant  and  this  complication  is  not  present.  Whenever  linear 
aerodynamic  damping  is  considered,  this  coupling  effect  must  be  considered. 

For  the  analysis  of  range  data  acquired  from  missiles  with  cubic  static 

moments,  the  frequency  equations  with  the  average  or  mid-range  values  of  the 

2 

proper  equivalent  squared  amplitude,  8  .  need  only  be  considered.  Modified 

^  J 

relations  for  the  sum  and  products  of  the  frequencies  now  follow  by  first 
eliminating  Mq  between  Equations  (2.9)  and  (2.11)  and  then  eliminating  P 


between  these  equations. 
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According  to  Equation  (2.14),  the  product  of  the  measured  frequencies  should 

2 

be  plotted  versus  the  effective  squared  angle,  8g  and  a  line  fitted.  The 
intercept  of  the  line  is  proportional  to  the  linear  term  in  the  static  moment 
while  its  slope  is  proportional  to  the  cubic  term. 
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8.  dubic  Static  Moment;  Exact  Solution 

In  order  to  obtain  the  exact  solution  of  Equation  (2.2),  it  is  necessary  to 

»s»  t 

eliminate  the  coupling  term  iPg  »  For  the  case  of  no  aerodynamic  damping  and  a 
constant  roll  to  velocity  ratio,  this  can  be  done  by  transforming  coordinates  to 
a  coordinate  system  which  is  rolling  with  roll  rate  p/2 

t  =  lei(i/2>Pe  (3.1) 
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Equation  (3*2)  has  a  very  nice  property  of  possessing  an  exact  solution. 
This  property  is  not  shared  by  the  other  nonlinear  equations  to  be  studied  in 
this  chapter.  Since  a  comparison  of  the  exact  solution  with  the  approximate 
solution  of  Section  8.2  will  establish  confidence  in  this  approximate  method, 
we  will  now  derive  this  solution. 

First  we  write  the  conjugate  of  Equation (3. 2) 
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^8  )g  = 
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Multiplying  Equation  (3.2) 
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Equation  (3-3)  by  I 


and  adding,  we  have 
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is  twice  the  sum  of  the  kinetic  and  potential  energies  and  is  constant  for 
the  conservative  cubic  moment.  A  second  first  integral  may  be  obtained  by 
eliminating  the  static  moment  between  Equations  (3.2)  and  (3.3). 
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(3.5) 


In  order  to  interpret  this  constant  of  the  motion  polar  coordinates  are 


introduced. 
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Thus  C9  is  twice  the  flight  path  component  of  angular  momentum.  These  equations 
may  be  combined  to  yield  simple  relations  for  8  and  & 


(s2)'  =  -  7-  cj;  +  i+C^2  +  +  2R3Sf 


(3.9) 
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The  solution  to  Equation  (3-9)  is  an  elliptic  function  of  the  first  kind®-^. 

A 

Of  the  four  possible  combinations  of  signs  for  Mq  and  Mg,  three  correspond  to 
periodic  motions.  These  three  are  illustrated  in  Figure  8.1  and  may  be  identified 
by  the  following  types : 

(a)  Stable  at  small  angles;  more  stable  at  larger  angles 
(Mq  <0,  Mg  <  0) . 

(b)  Stable  at  small  angles;  less  stable  at  larger  angles 

(Mo  <0,  Mg  >  0). 

(c)  Unstable  at  small  angles;  stable  at  larger  angles 
(Mq  >0,  Mg  <  0). 

These  periodic  solutions  may  be  written  in  terras  of  the  sine  amplitude 


function. 


type  (a)  moment 
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where  61  is  minimum  value  of  6 


50  is  maximum  value  of  6 
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The  inequalities  associated  with  Equations  (j. 12-13)  are  quite  important 
in  themselves.  Three  interesting  observations  may  be  made: 

1.  For  type  (a)  static  moment,  periodic  motion  of  any  amplitude  is 

possible . 

2.  For  type  (b)  static  moment,  symmetric  planar  periodic  motions 
(p  =  0,  6^=0)  are  possible  for  all  amplitudes  for  which  the  moment  does  not 
change  sign;  circular  motion  ( 5^  =  bg),  however,  is  possible  only  when  the 
nonlinear  part  of  the  moment  is  not  greater  than  two-thirds  of  the  linear  part. 


3.  For  type  (c)  static  moment,  possible  periodic  motions  are  those 

2 

for  which  the  median  value  of  6  yields  a  stable  moment.  The  variables  and 
&p  may  be  approximately  related  to  the  modal  amplitudes  and  Kg  by  the 
equations 


&2  =  (Kx  +  Kg)2 


(3.14) 
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(3.13) 


The  period  of  the  elliptic  sine  function  in  Equations  (J. 11-12)  corresponds  to 

a. 

half  the  period  of  0  in  Equation  (2.  7)  since  that  equation  can  be  put  in  the 
form 
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types  (a)  and  (c) 
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where  K(k)  is  the  complete  elliptic  integral  of  the  first  kind. 
(The  period  of  sn(tus,  k)  is  4K/oo. ) 

%  may  be  computed  from  the  quasi-linear  frequencies  of  Equations  (2.9) 

2  2 

and  (2.11)  for  two  special  motions;  circular  (&^  =  K^IC>  =  ^  °r  P-*-anar* 

(62  =0;  Kx  =  Kg). 

For  circular  motion 
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and  for  planar  motion 
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m  is  the  ratio  of  the  cubic  part  of  the  static  moment  to  the  linear  part. 
Equations  (3.18)  and (3. 19)  are  compared  with  the  corresponding  exact  equations 
for  different  moment  types  in  Figures  8.2  and  8. 3-  These  figures  show 
amazingly  good  agreement  of  the  approximate  curve  and  the  exact  one. 


This  motion  is  called  planar  although  except  for  zero  spin  it  is  not  planar 
in  range  coordinates.  In  other  words,  planar  epicyclic  motion  is  motion  which 
includes  zero  total  angle  of  attack. 


'  .  •<  EfiYvl  of  Large  Angles  on  Epicyclic  Frequencies 

When  large  angles*  are  considered,  two  quite  different  effects  b^corr/ 
important.  First  the  simple  geometric  relations  of  Section  7-1  between  the 
direct  ion  cosines  and  the  angles  of  attack  and  sideslip  must  be  re-examined 
and  then  the  derivative  of  the  cosine  of  the  total  angle  of  attack,  y1 ,  must 
be  retained  in  the  equations  of  pitching  and  yawing  motion. 

Although  Equations  (7.1.14=15)  must  be  used  to  calculate  the  transverse 
components  of  the  velocity  vector  in  the  fixed  plane  coordinates,  the  exact 
transformation  between  these  coordinates  and  the  non-rolling  coordinates  seems 
to  be  quite  formidable,  Qur  geometrical  intuition  may  provide  some  help.  The 
motion  of  the  2  and  3  axes  of  the  fixed  plane  system  is  directly  related  to  the 
missile’s  periodic  pitching  and  yawing  motion  and,  therefore,  the  average  roll 
rate  of  the  fixed  plane  system  should  be  zero  and  frequencies  measured  in  the 
fixed-plane  coordinates  the  same  as  those  measured  in  non-rolling  coordinates. 
To  verify  this  conjecture  for  our  cubic  analysis  we  must,  therefore,  show  that 
at  least  the  second  order  approximation  in  average  angular  velocity  of  the 
fixed  plane  axes  is  zero. 

Since  as  given  bv  Equation  (7  =  1.10)  is  at  most  second  order  we  can 

use  the  first  order  approximations  for  and  (Equations  (7.1.16)) 
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When  the  higher  order  terms  in  Equation  (7.1.10)  axe  omitted,  it  reduces  to 
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Angles  are  usually  thought  to  be  large  when  the  approximations 
0  =  sin  0  =  tan  6,  cos  ©  =  1  are  not  very  well  satisfied.  A  glance  at  a 
table  of  trigometric  functions  shows  that  angles  in  excess  of  fifteen 
degrees  cause  these  relations  to  be  in  error  by  more  than  %  and  could  be 
considered  to  be  large . 


Equations  (4.i-i)  may  now  be  combined. 
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The  trajectory  terms  in  Equations  (4.1,  2,4)  arise  from  the  lift  force  and 
are  much  smaller  than  the  first  terms  in  each  of  those  equations  and,  therefore, 
will  be  neglected.  A  good  approximation  for  the  transverse  velocities  is  oar 
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According  to  Equation  (4.7)  co^i/V  has  a  non-zero  average!  This  set  back 
to  our  geometrical  intuition  may  be  overcome  by  the  realization  that  Equation  (4.7 
allows  us  to  state  the  relation  between  frequencies  in  the  range’s  fixed  plane 

,«  t 

coordinates,  0jr ,  and  frequencies  in  the  theory’s  non-rolling  coordinates,  0 j . 
Since  the  measured  frequencies  are  actually  averages 
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The  differential  equation  for  a  cubic  static  moment  and  large  geometrical 
angles  can  be  obtained  from  Equation  (6-6.8) 
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Note  that  the  definition  of  VL  differs  from  that  following  Equation  (2.2).  This 
difference  arises  from  the  retention  of  the  cosine  in  the  definition  of  M.  If 
the  numerical  effect  of  this  difference  is  large,  this  means  that  the  non¬ 
linearity  in  M  is  primarily  a  geometrical  nonlinearity  rather  than  an 
aerodynamic  nonlinearity. 

2  l/P  t 

Since  y  -  (l  -  6  )  '  ,  the  y  term  can  be  quickly  computed 
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From  Equations  (2.7)  and  (k.13)  we  see  that  if  damping  is  neglected  in  comparison 
with  frequency 
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These  equations  may  now  be  substituted  in  Equation  (k.12),  solved  for  the 
damping  and  frequency  of  the  first  mode,  and  averaged  over  a  cycle  of  the 
difference  frequency.  The  resulting  equation  for  damping  is  the  same  as 
before  but  the  frequency  equation  for  the  first  mode  becomes 
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Similarly,  the  frequency  equation  for  the  second  mode  become. 
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Equations  (4.15-4.16)  differ  from  Equations  (2.9)  and  (2.11)  by  the  presence 
of  the  last  terms.  Equations  similar  to  Equations  (2.13-14)  may  now  be  derived 
by  first  eliminating  M  between  Equations  (4.15-16)  and  then  eliminating  P. 
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1  . 

Thus  the  y  /y  term  does  have  a  noticeable  effect  on  the  equations  for  sum 
and  product  of  the  frequencies  in  the  non- rolling  coordinates.  What  is 
measured,  however,  are  the  frequencies  in  range  coordinates.  If  Equations 
(4.10-11)  are  used  to  eliminate  the  non-rolling  frequencies  and  fourth  powers 
of  K-i’s  omitted,  we  have  the  following 

'J  rs  r\ 


$lx  +  ^2r 


*'lr  ■  4  (1  •  (?)  <4  +  4 >  >  =  Mo  + 


/  v2  ^2  \ 

V  '  ir  ‘  dv  > 

(4.19) 

(4.20) 

The  second  term  in  Equation  (4.l8)  comes  from  the  cosine  of  the  total  angle 

•  1 

of  attack  and  could  be  quite  large  in  comparison  with  0  .  02  since  it  involves 


the  square  of  the  higher  frequency  and  this  frequency  could  be  much  larger  than 
the  lower  frequency.  The  modification  of  Equation  (2.14)  due  to  range 
frequencies  and  the  cosine  terms  is  Equation  (4.20)  and  it  can  be  seen  that 
this  modification  is  quite  small.  The  two  effects  mentioned  at  the  beginning 
of  this  section,  therefore,  almost  cancel  each  other  and  should  be  either  both 
retained  or  both  omitted  from  the  relation  for  the  product  of  the  frequencies. 
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The  left  side  of  Equation  (4.20)  contains  an  additional  effect  of  thr 
cosine  of  the  total  angle  of  attack  through  the  definition  of  MU,. 
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If  Equation  (4.20)  is  multiplied  by 
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the  modal  amplitudes  neglected,  the  following  very  convenient  form  results. 
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Thus  the  range  value  of  the  static  moment  coefficient, 
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,  is  essentially 


the  product  of  the  range  measured  frequencies  and  a  set  of 
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»s  for  the  same  Mach  number  and  a  cubic  static  moment  should  yield  a 


straight  line  where  plotted  versus  6g.  The  intercept  is  cq,  the  zero  amplitude 
static  moment  derivative,  and  the  slope  is  c2,  the  cubic  coefficient  in  an 
expansion  of  the  static  moment. 

6.5  Cubic  Magnus  Moment 

As  could  be  seen  in  Section  8.4,  the  effect  of  large  geometric  angles  is 
rather  small  and  will  be  neglected  in  the  remainder  of  this  chapter.  Equation  (4.20) 
is  quite  simple  to  apply  and  will  be  used  in  the  data  analysis.  It  now  remains 
to  include  the  aerodynamic  damping  terms.  The  only  damping  nonlinearity  which 
has  been  measured  in  a  ballistic  range  has  been  a  cubic  Magnus  moment  and  we 
will  restrict  our  study  to  this  nonlinearity.  It  is  a  simple  exercise  for  the 
reader  to  incorporate  any  other  type  of  nonlinearity  he  may  prefer.  The  exact 
expressions  for  arbitrary  nonlinearities  will  be  derived  in  the  next  chapter. 
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For  cubic  static  and  Magnus  moments  Equation  (6.6.12)  assumes  the  form 
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The  quasi=linear  analysis  of  previous  section  results  in  the  same  frequency 
equations  [Equations  (2.9-11)]  but  different  damping  equations 
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One  unfortunate  feature  of  Equations  (5.2-5)  is  the  way  they  are  coupled- 
The  damping  equations  are  affected  by  the  derivatives  of  the  frequencies  but 
these  derivatives  in  turn  involve  the  damping.  If  Equation  (5-2)  is  differentiated, 
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It 

Equation  (5»b)  and.  a  similar  equation  for  0g  t>e  use<^  to  eliminate  the 
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202  -  P 


(5.7) 


(5-8) 


can  he  obtained  from  the  measured  frequencies  and  the  left  sides  of  above 

equation  can  then  be  computed  for  each  test.  In  most  cases,  however,  the  M„ 

term  is  quite  small  and  will,  therefore,  be  omitted  in  the  remainder  of  this 

discussion.  Equations  (5.7-8)  show  that  for  a  series  of  tests  of  the  same 

configuration  at  the  same  Mach  number  the  measured  damping  exponents  should 

2 

be  plotted  versus  5ei  .  If  the  points  fall  on  lines  whose  slopes  are  negatives 
of  each  other,  the  Magnus  moment  is  cubic  with  cubic  coefficient  determined  by 
the  slope. 


In  most  range  tests,  this  plot  is  not  made  but  the  linear  formulas 
(Equations  6.6-25-26)  are  used  to  compute  range  values  of  the  damping  in 
pitch  and  Magnus  moment  coefficients.  If  these  formulas  are  applied  to 
Equations  (5-7-8),  and  the  static  moment  is  assumed  to  be  linear,  the 
following  relations  result 


156 


M 


pa 


A  A  'd 

C  +  C0So 

o  2  e 


(5.9) 
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+  C 


M. 


=  C 
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CM.  +  C2A 
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(5.10) 


where  A  =  - 

x 


*t"  /  A  .  A' 

xy  /  0l  +  ^2 


x0‘  -  0*  , 


(K^  -  K^) 


It  is  interesting  to  note  the  reappearance  of  6  .  This  should  not  be  too 
surprising  since  the  Magnus  moment  is  the  imaginary  coefficient  of  I  while  the 
static  moment  is  the  real  part  of  the  same  coefficient.  Equation  (5-9)  contains 
a  warning  that  a  dependence  of  the  range-measured  damping  moments  coefficient  on 
amplitude  of  motion  does  not  necessarily  imply  a  nonlinear  damping  moment.  A 
nonlinear  Magnus  moment  can  have  this  effect. 

8.6  Cubic  Lift  and  Magnus  Forces 

The  swerving  motion  for  linear  aerodynamic  forces  has  been  discussed  in 
some  detail  in  Section  6.9.  For  simplicity,  we  will  neglect  damping  force, 
the  asymmetry  force,  gravity  and  Coriolis  acceleration.  These  quantities  may 
be  easily  incorporated  in  the  nonlinear  analysis  as  additive  terms.  In  range 
coordinates  then,  Equation  (6. 9*7)  reduces  to 


x2  + 


=  f  £ 


(6.1) 


where  f  =  CT  +  1  C.T 
L  V  N 

a  pa 


If  the  lift  and  Magnus  forces  are  cubic  in  5.  f  has  the  form  f  +  f-8  . 

'  o  2 

Integrating  6.1  with  the  excellent  approximation  s  =  -j-  ,  we  have* 


x9  +  lX^ 


s  s^ 

-  Bo  +  V  +f  f  <fo  +  f2s2)  T 


ds^ds-. 


o  0 


(6.2) 


Note  that  of  this  equation  are  the  negatives  of  of  Equation  (6.9.8), 
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If  f  were  linear,  this  reduces  to  the  form 


x2  +  xx. 


=  B.+B,.s+fL 
ol  1/  o  i 


(6.3) 


where  1^  ~  J  J  '  ^s2~si 

o  o 

Now  motion  caused  by  nonlinear  forces  has  been  well  fitted  by  Equation  (6.3) 
and  provided  a  range  value  of  f. 


f  =  f  +  f_8“ 
r  o  2  es 


(6.4) 


2  2 

where  5  is  an  effective  value  of  6  for  the  swerving  motion. 

6S 

2 

It  remains  to  determine  this  effective  value  of  6  .  To  do  this  we  will  assume 
that  I  is  well  approximated  by  an  epicycle  without  damping. 


h  ~  Co l  +  CliS  +  1 1 
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Therefore  the  equation  for  the  cubic  forces  becomes 


x_  +  ix_ 

2  A  A 

- j — -  =  B  +  B,s  +  f  I.  +  U 

l  OX  O  Jt  d.  C 


(6Ji) 


where  B .  =  B .  +  f  C. .  +  f _C . 

1  1  o  it  2  1 


In  order  to  reduce  Equation  (6.8)  to  the  form  of  (6.6),  the  effective  5^  would 
have  to  satisfy  the  eauation 


5  1.  -  I  =  0 

es  i  c 


(6.9) 


Except  for  certain  special  cases,  I  is  not  proportional  to  I,  and,  hence,  a 
2  c  Xt 

constant  5  will  not  satisfy  Equation  (6.9)*  Our  experience  in  data  analysis 

es  /  2 

implies  that  it  should  be  possible  to  select  a  constant  value  of  5  .  The  key 

es 

to  the  solution  lies  in  the  fact  that  the  data  analysis  is  based  on  a  least 

2 

squares  technique.  We  will,  therefore,  define  5  to  be  that  number  which 

es 

makes  the  integral  over  the  observed  trajectory  of  the  squared  magnitude  of  the 
left  side  of  Equation  (6.9)  a  minimum.  This  approach  is  quite  similar  to  that 
used  for  a  quadratic  drag  coefficient  in  Section  5*3* 

2 

6gs  is,  therefore,  defined  by  the  requirement  that 


L/2 


-L/2 


2  2  — 

(5  -  O  ^  is  a  minimum 

6S  L  C  es  X/  c 


We  differentiate  the  above,  set  the  result  equal  to  zero  and  solve  for  5 
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For  large  L  the  numerator  is  the  average  value  of  II. +11,  and  the 

p  C  AJ  AS 

denominator  is  twice  the  average  of  I.I4.  6. _  can  be  quickly  computed  and 


i  t 


es 


reduced  to  the  fore 


«OM6ei +  (<>444 


• 1  r  >i  •  »s  r> 

(<OX  +  X)  kJ 


(6. ii) 


For  the  special  case  of  a  spin  stabilized  missile  0?  <  <  0  and  the 
-03  to  the  simple  fux’in 


(6.12) 


This  is  one  of  the  special  cases  which  directly  satisfies  Equation  (6.9)*  A 

I  » 

second  important  case  is  that  of  a  nonspinning  missile  (0  =  -  0g)-  For  this 

condition 


8 


2 


es 


(6.13) 


In  any  event  the  range  values  of  CT  and  CM  should  be  plotted  versus  o  and 

L/  N  _  ss 

a  pa 

the  existence  of  a  cubic  force  dependence  on  5  determined. 


8.7  Experimental  Results 

As  can  be  seen  from  the  preceding  analysis,  the  cubic  coefficients  can 
only  be  obtained  from  a  set  of  tests  of  the  same  configuration  at  the  same 
Mach  number  but  flying  at  different  effective  angles  of  attack.  These 
considerations  lead  to  a  need  for  at  least  three  and  preferably  five  tests  in 
each  set.  Although  the  analysis  has  been  applied  to  the  data  obtained  from 
finned  missiles  these  programs  were  not  large  enough  to  allow  cross  checks  of 
the  results  and  no  good  wind  tunnel  data  was  available. 

8-6 

The  ogive- cylinder  program  of  Section  7-*+  (U)  included  tests  of 

identical  models  with  three  different  center  of  gravity  locations  at  the  same 

Mach  number.  This  redundancy  of  data  can  and  will  be  used  to  determine  the 

8  T 

validity  of  the  cubic  analysis,  In  a  second  program  the  large  angle 
behavior  of  a  20mra  shell  was  specifically  studied.  Some  of  the  results  of 
this  program  can  be  used  to  increase  further  our  confidence  in  the  cubic 
analysis. 


Witul  tunnel  measurements  of  the  Magnus  force  and  moment  for  the  seven 
caliber  ogive  cylinder  and  static  force  and  moment  for  the  20mm  shell 
are  available  and  will  be  compared  with  the  results  of  the  nonlinear  reduction 
at  the  same  Mach  numbers  with  different  center  of  gravity  locations.  In 
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angle  of  attack  for  the  nine  caliber  long  models  at  Mach  number  of  1.8.  As 
can  be  seen  from  this  figure,  the  points  for  each  center  of  gravity  location 
fall  on  lines  with  well  determined  slopes. 


Of  the  twenty-seven  possible  combinations  of  length,  Mach  numbers,  and 
center  of  gravity  location,  the  data  was  sufficiently  good  to  determine  sixteen 
values  of  the  cubic  coefficient,  c (Table  8-1)  For  some  of  the  Mach  numbers 
and  length,  wind  tunnel  measurements  of  the  static  moment  were  available.  These 
were  fitted  by  cubics  and  the  cubic  coefficients  are  listed  in  Table  s-  1.  At 
all  seven  points  of  comparison  the  agreement  is  good. 


Since  models  with  three  different  center  of  mass  locations  were  tested, 
a  second  check  is  possible.  If  the  normal  force  is  expanded  as  a  cubic  section 
of  6,  the  usual  center  of  gravity  relation  (Table  6-1)  provided  that 
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where  f!  =  a  +  a  a  and  r  7  s  ^ 

N  o  2W  ‘  l  eg' 

a 

are  the  moment  coefficients  for  a  center  of  mass  located  s  calibers  forward 

eg 

of  that  for  the  c^’s. 
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In  A 


or  differ  ‘6 n't  center  of  mass  locations  ax’tj  linear 

X 

functions  of  location.  In  Figure  8.5>  the  slopes  and  intercepts  of  Figure  8.4 
are  plotted  versus  center  of  gravity  location.  The  fact  that  the  cubic 


coefficients  as  well  as  the  linear  coefficients  fall  on  a  straight  line  is 
another  point  in  the  favor  of  the  cubic  analysis.  The  slopes  of  these  lines 
are  a 2  and  respectively.  The  c.g.  intercept  of  the  cn  line  is  the  location 


l6l 


of  the  center  of  pressure  and  the  c.g.  intercept  of  the  line  is  the  c.g. 
location  for  which  the  cubic  coefficient  is  zero  and  the  moment  is  exactly 
linear.  This  point  is  rear  of  the  rear  center  of  gravity  tested. 


For  a  total  of  seven  of  the  nine  combinations  of  length  and  Mach  number 
it  was  possible  to  determine  a_  from  the  variation  of  c2  with  c.g.  location. 
The  wind  tunnel  measurements  of  normal  force  may  be  analyzed  to  yield  cubic 
coefficients  as  well  as  linear  coefficients.  These  are  tabulated  in  Table  8-2 
and  the  agreement  with  the  range  values  is  seen  to  be  reasonably  good. 


If  the  lift  force  is  cubic  in  angle  of  attack,  it  can  be  written  in  the 

form 


LO 


a!26 


(7-3) 


relation  connecting  the  lift,  di-ag,  and  normal  forces  (Equation  (6.6.7)), 
the  coefficients  of  the  normal  and  lift  forces  satisfy  the  following  equations. 
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In  Section  8.o,  it  is  shown  that  the  range  value  of  the  lift  coefficient  of  a 

2 

spinning  missile  should  be  plotted  versus  6o2  and  the  cubic  coefficient  can  be 
obtained  from  the  slope  of  a  line  through  the  data.  This  technique  is 
•illustrated  in  Figure  8.6.  Values  of  aQ  and  a2  could  be  computed  for  eight 
combinations  of  length  and  Mach  number  and  these  are  compared  with  values  obtained 


from 
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Finally,  in  the  20mm  shell  program  it  was  possible  to  study  the  influence 


of  quite  large  angles.  In  Figure  8.7,  the  range  static  moment  coefficient  is 
plotted  versus  effective  squared  angle  of  attack  and  we  see  that  the  data  is 
essentially  bilinear.  Each  line  corresponds  to  a  cubic  segment  in  the  moment 
plans.  If  the  parameters  of  each  cubic  are  calculated  from  the  slope  and 
intercept  of  its  corresponding  line  in  Figure  8.7,  they  can  be  pieced  together 
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to  form  a  smooth  moment  plot.  (Figure  8,8),  An  examination  of  the  spark  shadow¬ 
graphs  reveals  that  flow  separation  occurs  at  about  21°  and  this  explains  the 


sudden  change  in  the  moment  curve  at  this  point.  In  Figure  8.9, 


a. 


is  plotted 


versus  8  0  and  the  corresponding  lift  force  is  shown  in  Figure  8.10. 


Since  wind  tunnel  measurements  were  available,  a  further  verification  ic 
possible.  The  data  were  divided  into  that  for  angles  less  than  21°  and  that 


for  angles  greater  than  21  and  pairs  of  cubi 


X  X.  u  uC 
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coefficients  are  given  in  Table  8.3*  The  agreement  for  such  different 
experimental  conditions  is  quite  good*. 

In  view  of  this  success  with  cubic  static  moments,  our  interest  turns 
quite  naturally  to  cubic  Magnus  moments.  According  to  Equations  (5. 9-10)  the 
range  values  of  damping  in  pitch  and  Magnus  moment  coefficients  are  both  affected 
by  the  presence  of  a  cubic  Magnus  moment.  Thus  values  of  this  cubic  coefficient 


may  be  computed  from  a  plot  of 


pa 


and 


r 


'IA 


+  C. 


M, 


versus  8  and  A 
6 


respectively.  (See  Figures  8.11-12)  Eight  cubic  Magnus  coefficients  were 

determined  for  the  ogive  cylinder  and  these  results  appear  in  Table  S-4.  In 

four  cases  A  had  a  wide  enough  range  to  determine  c0  from  C„  +  C.  and 

d  MM. 

L  8  GM 

A 

the  pairs  of  values  of  are  all  within  their  combined  probable  errors. 
Unfortunately,  with  the  exception  of  the  one  value  for  the  nine  caliber  long 
models  these  probable  errors  are  quite  large.  Thus,  the  measurements  of  cubic 
Magnus  moments  are  quite  delicate. 

The  swerving  motions  of  the  ogive  cylinder  were  examined  to  determine  any 
influence  of  a  cubic  Magnus  force.  In  only  two  out  of  the  six  combinations  of 
length  and  Mach  number  was  this  effect  large  enough  to  measure  and  these  two 
cases  are  shown  in  Figure  8.13.  The  cubic  Magnus  moment  coefficients  given  in 
Table  8 -4  for  different  center  of  gravity  locations  could  be  analyzed  to  yield 
four  indirect  values  of  the  cubic  Magnus  coefficient,  a^.  The  agreement  of  these 
coefficients  is  reasonably  good. 


Additional  examples  of  the  analysis  of  cubic  static  moments  are  supplied 
by  References  8-10  and  8-11. 


163 


TABLE  ft- 1 


CutflC  GtAtIC  MOMENT  COEEFlCIEwr,  c 2,  FOR  OGIVE  CYLINDERS 

s  =rt d2/4;  i  =  a 


Mach 

Forward  c.g. 

Middle  c.g. 

Rear  c.g. 

Number 

Range  Wind  Tunnel 

Range  Wind  Tunnel 

L/d  =  5 

Range 

Wind  Tunnel 

1.3 

c 

-  0 

-8-2  0 

~  ^  0 

2.5  -5  0 

l.S 

-  36 

-  18 

8-2 

0 

L/d  =  7 

i  * 

id  +  c 

T  Z 

-  ^ 

■"U  —  J 

■O 

1.8 

-31  -  3 

-  28 

-  8 

CO 

+  1 

O 

8 

2.5 

-71  -  3 

-  66 

-  28 

L/d  =  9 

1.3 

KA 

+  1 

ON 

QO 

• 

0.3-  .03 

1.8 

-76  -  10 

-41  1  5 

-9.5-  -8 

2.5 

-76  -  5 

-2.5  -  1.5 

T  ^  )i 
J.W*f 


TABLE  8-2 


LINEAR  AND  CUBIC  NORMAL  FORCE  COEFFICIENTS  FOR  OGIVE  CYLINDERS 


Number 

0 

2 

c.g. 

Swerve 

Wind  Tunnel 

T  /  O  C 

li/U  =  J 

c.g. 

Swerve 

Wind  Tunnel 

1.3 

2.50 

2.52 

2.34 

17 

13-0.5 

12 

1.8 

2.88 

2.88 

2.80 

- 

r,r-  +  ,, 

25  -  2 

17 

1.3 

2.60 

2.50 

I'¬ 

ll 

'll 

s, 

1 

51 

25 

+ 

3 

1.8 

2.88 

2.88 

2.74 

31 

23 

+ 

10 

25 

2.5 

3.08 

3.06 

3.08 

46 

69 

j. 

10 

46 

L/d  =  9 

1.3 

2.70 

2.72 

43 

43  i  3 

1.8 

2.90 

3.13 

36-1 

38-5 

2.5 

3.33 

3.06 

87 

76  i  5 
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TABLE  8-3 


BI- CUBIC  STATIC  MOMENT  AND  LIFT  COEFFICIENT  FOR  20mm  SHELL 


c 

o 

C2 

aL0 

&TO 

±J- 

Wind  Tunnel 

2.10 

-2.8 

2.57 

10.2  'j  before  separation 

Range 

2.11 

-2.5 

2.52 

9.2  J  6  =  sin  21° 

Wind  Tunnel 

1.86 

0.0 

3.69 

0.0  ^  after  separation 

Range 

1.83* 

0.3 

3.64* 

f 

1.0  ->  6  =  sin  21" 

Evaluated  at  separation  angle  of  21°.  The  20mm  shell  flying  at  an  angle  of 
25  with  separated  flow  is  shown  in  Figure  8.14. 
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TABLE  8-4 


CUBIC  MAGNUS  FORCE  AND  MOMENT  COEFFICIENT  FOR  OGIVE  CYLINDERS 


T  /A  -  7 


°2  C2  C2  a2  a2 
Mach  No.  f.c.g.  m.c.g.  r.c.g.  c.g.  swerve 

1.3  fcM  1  40  -  13  13  -  3  —  -  36  -  42  ±  3 

L  r 

K  +C  1  50  -  36 

ri  w# 
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1.8 

CM 
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2.5 

1  1 

a 

8s 

1  1 

90-13  —  15-8 

-  48 

— 

CM  +  CM. 
q.  a 

L  J  r 
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1.8 

\  1 
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CM  +  CM. 
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FIG.  8.2  A 
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FIG.  8.2 B 


171 


PLANAR  MOTION 
TYPE  (c)  MOMENT 


EXACT 


CIRCULAR  MOTION 
TYPE  (a)  MOMENT 


o 


5 


m 

175 


FIG.  8.3  A 
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FIG.  8.3  B 
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CHAPTER  IX 


FflEDICTION  OF  NONLINEAR  MOTION 

9.1  Introduction 

In  the  last  chapter,  we  considered  the  measurement  of  certain  nonlinear 
forces  and  moments  by  analyzing  the  motion  induced  by  these  quantities  an  j 
observed  on  a  ballistic  range.  In  this  chapter,  we  consider  the  invars* 
problem  of  predicting  the  motion  caused  by  known  nonlinear  forces  and  moments. 
The  actual  expression  for  this  force  and  moment  may  be  either  measured 
experimentally  in  a  ballistic  range  or  wind  tunnel  or  computed  from  a 
theoretical  model. 


The  problem  will  first  be  attacked  by  a  generalization  of  the  quasi- linear 
method  of  Section  8.2.  This  method  provides  effective  values  of  frequencies 
and  damping  exponents  in  terms  of  "averages"  of  the  linear  and  nonlinear 
aerodynamic  terms  but  is  limited  by  the  assumption  that  the  motion  for  ui  least 
a  number  of  nutations  is  epicyclic.  The  behavior  for  a  large  number  of  cycles 
is  determined  by  the  consideration  of  the  modal  amplitudes  alone  and  the  use 
of  an  "amplitude  plane" .  Finally,  the  limitation  to  epicyclic  motion  is 
removed  for  the  special  case  of  a  strongly  nonlinear  cubic  static  moment.  In 
this  case,  the  trigonometric  functions  of  the  epicycle  are  replaced  by  elliptic 
functions  and  use  is  made  of  "quasi-constants”  of  the  motion. 


9-2  The  Quasi -Linear  Analysis 

Quasi-linear  analyses  such  as  that  of  Kryloff  and  Bogoliuboff  usually 
become  extremely  complex  when  applied  to  general  real  fourth  order  differential 
systems.  We  have  already  seen  that  the  linearized  motion  of  a  rotationally 
symmetric  missile  can  be  described  by  a  very  simple  fourth  order  system  --  one 
that  can  be  written  in  the  form  of  a  second  order  equation  in  one  complex 
variable. 


The  nonlinear  analysis  will,  therefore,  be  limited  so  that  this  symmetry 
is  preserved.  This  is  the  reason  symmetric  direction  cosines  (u/V,  v/v,  w/v) 
and  not  nonsymmetric  Eulerian  angles  were  used  to  locate  the  velocity  vector 
with  respec  t  to  the  mi  ssils  oXis  system.  In  this  chapter y  the  agiOuyu&mi  c 
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coefficients  are  assumed,  to  be  functions  of  the  complex  i,  and  its  derivative, 

,  which  are  independent  of  roll  orientation.  Clearly,  the  magnitude  of  these 
two  two-dimensional  vectors  and  the  orientation  angle  between  them  possess  this 
property.  Now  the  angular  velocity  can  be  resolved  into  a  radial  component 
along  the  complex  angle  of  attack  and  circumferential  one  perpendicular  to  the 
complex  angle  of  attack,  and  one  of  these  plus  the  magnitude  of  the  angular 
velocity  is  sufficient  to  determine  this  orientation  angle.  With  these  remarks 
in  mind,  we  will  introduce  the  rotationally  invarient  quantities 

B2  =  j  I  j2  =  I  I  ,  |  I’  I"  =  6*  I',  and 

(52)'  =  266*  =  (IT*  +  S'!). 


and  assume  that  aerodynamic  coefficients  are  at  most  functions  of  these 
quantities . 


The  actual  quasi-linear  analysis  initially  parallels  that  of  Section  6.7 
for  varying  coefficients.  (Equation  (6.7.8-10)  ).  As  in  that  section  we  assume 
that  the  solution  can  be  written  in  the  form  of  an  epicycle  whose  frequencies  are 
not  necessarily  constants  and  whose  amplitudes  are  not  necessarily  exponentially 
varying. 
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b 


(2.1) 


Differentiating  Equation  (2.1)  twice  we  have 
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(2.5) 


Since  Equation  (2.l)  is  the  solution  to  the.'  homogeneous  linearized 
equation  of  pitching  and  yawing  motion,  we  will  apply  this  quasi-iinear  i;h  '.} 
to  the  homogeneous  exact  equation: 


£  +  (H  -  -  iP)  V  -  (M  +  iPT)  7=0 


(2.4) 


where  H,  M,  and  T  may  be  functions  of 
S  ,  | I  !  ,  and  (5  )  . 

Inclusion  of  the  small  inhomogeneous  gravity  term,  G,  in  Equation  (2.4)  imposes 
great  difficulties  on  the  nonlinear  analysis  and,  hence,  will  not  be  done. 

Equations  (2.1-3)  are  now  substituted  in  Equation  (2.4)  and  the  result 
rearranged. 

(0nr  -  -  X  (X,  +  H  )  -  X* 
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where  0  =  0.,  -  0^ 


For  the  linear  equation,  the  right  side  of  equation  (2.5)  is  fluctuating 
with  an  average  very  close  to  zero  and  the  left  side  is  at  most  slowly  varying. 
In  Section  6.7,  the  frequency  and  damping  exponent  for  the  one-mode  were 
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obtained  by  setting  the  left  side  of  Equation  (2.5)  equal  to  zero.  In 

Section  .2,  the  quasi-linear  solution  for  no  damping,  small  geometrical 

angles  and  cubic  static  moment  was  derived  by  aver  aginfc  the  light  side  of 

Equation  (2.5)  with  T  =  H  =  7  =  T  =H  =0.  In  order  to  incorporate  these 

o  o 

previously  excluded  terms,  we  will  take  the  average  of  all  the  terms  on  the 
right  under  the  quasi-linear  assumption  that  these  averages  are  good  approximation: 
of  the  effect  of  "small"  nonlinearities.  The  small  terms  involving  X.  in  the  real 
part  of  Equation  (2.5)  will  be  omitted  in  comparison  with  M_ . 
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It  was  shown  in  Section  8.4  that  7  /y  is  well  approximated  by  an  odd  function 
of  0  and,  hence,  onJfy  affects  the  frequency.  With  this  in  mind,  we  separate 
Equation  (2.6)  into  real  and  imaginary  parts. 
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Similar  expressions  apply  for  the  other  mode  except  the  algebraic  sign  in  front 

A  * 

of  the  sine  0  term  in  is  reversed  due  to  the  oddness  of  the  sine  function* 


f 

Equation  (2.7)  can  be  solved  for  0..  and  differentiated  to  obtain  0  . 

t  i  if  -L 

Since  0  depends  on  amplitude  when  nonlinearities  are  present  0  will  contain 
X,'s.  Thus,  Equation  (2.8)  is  not  solved  for  X-,  since  0^  also  has  X^  as  veil 
as  X2  in  it.  An  important  simple  relation  for  frequency  is  that  derived  in 
Section  8.2  for  cubic  static  moment,  small  geometrical  angles,  and  no 
aerodynamic  damping.  The  same  relation  may  be  obtained  from  Equation  (2.7)  if 
the  no  aerodynamic  damping  condition  were  replaced  by  the  requirement  that  the 

A 

average  of  (H  -  FT)  sin  0  is  zero.  This  would  occur,  for  example,  if  H  and  T 

2 

were  only  functions  of  6  .  Under  these  conditions  then 
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Equation  (2.9)  can  be  differentiated  and  substituted  in  Equation  (2.8). 

2.  ~\  r  _2 

mgKi  J  2m2K2 

1  +  f  -1  X„ 


2^1  +  +  2K^) 


1  +  m2(K^  +  2KJ 


* 

X1 


Iff  |  p 

M0  -  P  ^  *^<4  +  2KJ) 

“  +  - - 


r  -T1 

[l  +  ^(K^  +  2K^)J 


(2.11) 


191 


A  similar  equation  can  be  derived  for  the  other  mode. 
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Equations  (2.11-12)  may  be  solved  simultaneously  for  X-,  and  . 
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where  the  a.,  !s  are  defined  in  Table  Q-i* 

J4 

One  important  occasion  for  which  Mq  and  M2  Eire  variable  is  that  of 
missile  entering  or  leaving  the  earth’s  atmosphere.  (See  Section  6.7). 
this  missile’s  static  moment  coefficients  are  constants. 
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^  is  the  angle  the  trajectory  makes  with  a  downward  pointing 
vertical. 
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Two  distinctive  properties  of  a  nonlinear  homogeneous  system  are  the 


dependence  of  the  type  of  motion  on  initial  conditions  and  the  existence  of 
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essentially  the  same  motion  for  any  initial  condition  since  any  solution  is  a 
linear  combination  of  two  linearly  independent  solutions  and  the  only  possible 
limit  motion  is  the  trivial  one  of  the  zero  amplitude  motion  for  dynamically 
stable  systems. 

For  a  second  order  system  with  coefficients  which  are  independent  of  the 
independent  variable,  the  phase  plane  is  the  best  method  for  studying  these 
properties.  In  a  phase  plane,  the  dependent  variable  is  plotted  versus  its 
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exception  of  singular  points,  every  point  has  only  one  curve  through  it.  The 
initial  values  of  the  variable  and  its  derivative  locate  the  particular 


trajectory  which  is  the  solution  for  those  initial  vali 
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the  phase  plane  immediately  reveals  the  relationship  between  initial  conditions 
and  the  various  possible  motions  and  yields  boundaries  on  these  conditions  for 
particular  motions.  The  existence,  size,  and  shape  of  limit  motions  are  also 
revealed  by  this  two  dimensional  phase  space. 

A  fourth  order  system  has,  in  general,  a  four  dimensional  phase  space. 

One  selection  of  these  four  coordinates  would  be  the  magnitudes  and  orientation 
angles  of  the  initial  complex  angle  of  attack  and  initial  complex  angular 
velocity.  Since  rotationally  symmetric  missiles  should  have  a  rotationally 
symmetric  phase  space,  their  phase  space  is  essentially  three  dimensional.  Its 
coordinates  could  be  the  magnitudes  of  the  initial  angle  of  attack  and  the 
initial  angular  velocity  and  the  angular  orientation  of  one  with  respect  to  the 
other. 


A  three  dimensional  phase  space  is,  however,  quite  difficult  to  use. 
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indicate  a  possible  two  dimensional  space  which  might  he  used.  This  simpj/''ity 

in  dimensionality  is  paid  for  by  the  assumption  that  the  actual  nonlinear 

motion  can  be  described  by  epicycles  whose  frequencies  and  amplitudes  vary 

throughout  the  missile's  flight.  Since  the  general  behavior  of  the  raotior.1  9nd 

the  existence  of  limit  motions  are  revealed  by  Equations  (3 .  1),  the  frequency 

ec nations  may  be  neglected  and  only  the  two  equations  for  the  modal  amp.jtu«f$S 
* 

considered. 

As  was  previously  noted,  phase  planes  are  used  for  equations  which  dc  not 
contain  functions  of  the  independent  variable  other  than  the  dependent  variable 
appearing  in  the  phase  plane.  For  this  reason,  we  require  that  the  X  1  -  be 

»  i  i 

functions  on  only  the  modal  amplitudes  themselves,  i.e.,  M  =  =  P  =0,  etc. 

From  basic  symmetry,  we  see  that  the  squares  of  amplitudes  should  appear  in  the 

X.'s.  We,  therefore,  rewrite  Equation  (3.1)  in  the  form 
J 

(K^)‘  =  2K^X1  (K*,  K^)  (3, V 


(k|) *  =  2K^x2  (k^,  y*) 


(3.3) 


The  independent  variable,  s,  may  now  be  eliminated  by  dividing  Equation  (3*3) 
by  Equation  (3.2) 


(K^,  K^) 


hh  (? 


(5.4) 


Equation  (3.4)  describes  the  angular  motion  by  means  of  a  point  moving  in  the 
plane  which  we  will  call  the  amplitude  plane.  For  any  point  in  this 
plane,  the  quasi-linear  frequencies  and  damping  exponents  may  be  computed  from 
Equations  (2.7-8)  for  the  one-mode  and  similar  equations  for  the  two-mode  and 
except  for  the  phase  angles,  the  motion  is  completely  determined. 

Although  Equation  (3.4)  is  a  nonlinear  equation,  it  is  first  order  and, 
therefore,  much  simpler  than  the  fourth  order  Equation  (2.4).  The  amplitude 
plane  could  be  generated  by  numerically  integrating  Equation  (3*4)  for  a 
variety  of  initial  modal  amplitudes.  Equations  of  this  form,  however,  have 

*  A  rigorous  justification  for  the  quasi-linear  treatment  of  limit  motions  is 
given  in  Reference  9-13* 
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9-2 

been  extensively  studied  by  H.  Poincare  who  showed  that  the  basic  behavior 
of  the  amplitude  plane  trajectories  is  fixed  by  the  location  and  type  of 
Equation  (3.4)  singularities.  (Singularities  are  points  for  which  both  the 
numerator  and  denominator  of  the  right  side  of  Equation  (3.4)  vanish.)  A 
consideration  of  Equations  (3*4)  indicates  the  existence  of  three  categories 
of  possible  singularities: 

(1)  the  origin  (K^  =  =  0) 

(2)  the  intercepts  of  the  one-mode  zero  damping  curve  (K^  =  0,  =  0) 

(3)  the  intercepts  of  the  two-mode  zero  damping  curves  (K^  =  0,  =  0) 

(4)  the  intersection  of  the  zero  damping  curves  (x^  =  0,  =  0)  . 

According  to  Poincare’s  classification,  there  are  four  possible  types  of 
first  order  singularities:  nodes,  saddles,  spirals,  and  centers.  In  Figure  9-2, 
a  node  and  a  saddle  are  shown  while  a  spiral  appears  in  Figure  9-3c  and  a 
center  in  Figure  9-3d.  The  particular  curve  along  which  the  actual  motion 
moves  is  determined  by  the  initial  modal  amplitudes  as  specified  by  initial 
conditions.  If  the  coordinates  are  translated  so  that  the  singularity  is  at 
the  origin,  Equation  (3*4)  has  the  form 

& 

dx 

where*  ad  -  be  ^  0  and  P  and  Q  vanish  to  at  least  the  second  order  at  the 

origin.  The  criteria  for  the  type  of  singularity  may  now  be  stated  in  terms  of 

2 

the  coefficients  a,  b,  c,  d  and  their  discriminant  D=(b-c)  +4  ad.  (See 

page  44  of  reference  9-3-): 

(1)  The  singularity  is  a  node  If  (a)  D  >  0  and  ad  -  be  <  0  or  (b)  D  =  0. 
(Note  that  if  ad  =  0,  it  is  a  node  if  be  >0.) 

(2)  The  singularity  is  a  saddle  if  D  >  0  and  ad  -  be  >  0. 

(3)  The  singularity  is  a  spiral  if  D  <  0  and  b  +  c  /  0. 

* 

If  ad  -  be  =  0  and  the  origin  is  a  singularity,  a*b*c=d=o  and  the 
singularity  is  at  least  second  order. 
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(V)  The  singularity  may  he  a  center  if  L  <  0  and  b  i  c  -  0;  ot.hc.-rv;  , 
it  is  a  spiral*. 

i 

For  the  linear  case,  the  X.  are  constants  and  the  origin  is  the  onJ j 

J 

singular  point.  If  the  spin  is  zero,  \  =  \  ,  and  all  the  trajectories 

either  toward  or  away  from  the  origin  which  is  a  node.  (Gee  Figure  9-x.;, ; 

The  particular  curve  which  describes  a  given  motion  is  specified  by  initia 

values  of  K.  and  the  direction  of  motion  along  the  curve  by  the  sign  of  the 

i  J  » 

X.o.  For  positive  X.s,  the  motion  undamps  away  from  the  origin.  When  spin 
J  J 

is  not  zero,  the  damping  exponents  are  not  necessarily  equal  and  may  even  u ! -'Tor 
in  sign.  The  origin  can,  therefore,  be  either  a  node  (Figure  9~2a)  or  a  saddle 
(Figure  9-2b).  All  these  figures  display  an  important  characteristic  of  linear 
systems  --  i.e.,  the  type  of  motion  is  indenendent  of  initial  conditions. 


AwitvI  i  r\l  o? 


now  be  constructed  for  two  siuiple  nonlineari  f  x  -s 


and  one  fairly  general  nonlinearity.  In  the  first  example,  the  effect  ox 
quadratic  dependence  in  5  of  C„  +  C  on  the  motion  of  a  nonspinning  missile 

X'l  XM, 

q  & 

will  be  considered  while  in  the  second  the  effect  of  a  cubic  Magnus  moment  on 
the  motion  of  a  spinning  missile  will  be  described.  In  both  of  these  exump  ... , 

t 

small  geometrical  angles  (y  =  1,  y  -  0)  and  a  linear  static  moment  will  be 
assumed.  In  the  third  example,  the  effect  of  a  cubic  static  moment  and  a  very 
general  nonlinear  aerodynamic  damping  moment  on  a  nonspinning  missile  will  be 
studied  in  detail. 

..  _  „  ,  9-11 

9oa  Zero  Spin  -  Quadratic  Cw  +  C, 


For  this  case,  the  damping  in  pitch  coefficients  have  the  fc 


CM  +  CM.  =  do  +  d25 
a  n 


O.fc) 


i  +  (H0  +  HgS^)! '  -  Ml  =  0 


(  -j ,  7 ) 


If  P  =  Q  =  0,  the  singularity  is  a  center.  In  reference  9-4,  the  conditions 
for  a  center  are  given  when  P  and  Q  are  quadratic  functions  of  x  and  y. 

Since  equalities  seldom  occur  in  practice  --  only  inequalities,  the  distinction 
between  a  center  and  a  very  slow  spiral  is  purely  academic. 


where  H 


pSl 

"m 


a 


cn  "  kt2  d 
D  t  o 


*2 


pSJ  v“2  ,, 
2m  \  d2 


From  Equation  (2.8)  and  the  corresponding  equation  for  the  two-mode, 

/i\  L  .  - 

)  I  no  +  “2 


*1  ~ 


(3.8) 
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\z)  L“° 


F  £ 

“2c! 


/  7  \ 

(3.9) 


The  amplitude  plane  for  the  interesting  case  where  the  damping  changes 
sign  (HqH2  <  0)  is  shown  in  Figure  9-lb.  The  zero  damping  curves  are  vertical 
and  horizontal  straight  lines  with  intercepts  (-HQ/l 0)  and  (0,  -  Hq/h  ) .  The 
tests  of  this  section  show  that  both  of  these  singularities  are  saddle  points. 
These  lines  intercept  at  a  planar  singularity  point  (-Hq/ E^,  -Eq/e^)  which  is  a 
node. 

If  Hq  is  positive,  the  origin  is  a  stable  node  and  Figure  9- lb  shows  a 
characteristic  nonlinear  dependence  of  the  type  of  motion  on  initial  conditions. 
For  initial  modal  amplitudes  which  are  inside  the  square  formed  by  the  axes  and 
the  zero  damping  lines,  the  motion  will  damp  to  zero  amplitude;  for  modal 
amplitudes  outside  this  square,  the  motion  will  grow  into  an  undamped  spiral 


If  Hq  is  negative,  the  origin  becomes  an  unstable  node  and  the  planar  node 
becomes  a  stable  node.  For  any  initial  conditions,  the  motion  will  become  a 

.  T  /o 

planar-  motion  with  maximum  amplitude  =  2  | Hy/H,,  | x/ c .  This  limit  motion 

is  another  property  of  nonlinear  systems.  It  is  interesting  to  note  that  had 

id 

we  restricted  ourselves  to  planar  motion  I  =  be  o,  the  corresponding  differential 
equations  would  be  that  for  the  van  der  Pol  oscillator. 


4-  fU 


+  HgB2)  5*  -  M5  = 


/  T  «  M  \ 

VP.xu; 
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This  planar  singularity  is  precisely  the  limit  cycle  predicted  by  the  one--!  gr< 
of- freedom  theory.  This  example  is,  therefore,  a  generalized  van  der  Pol 
equation. 

9- 3b  Cubic  Magnus  Moment. 

For  this  example,  spin  is  constant  but  non-zero,  damping  in  pitch  ant. 
static  moments  are  linear,  ana  the  Magnus  moment  is  cubic  in  6. 


_  A  A  2 

CM  *  Co  +  =26 

pa 

i  +  (H  -  iP)l  -  I  M  +  i 


P(T  +  T06  ) 
O  d 


I  =  0 


(>  u) 

(3.12) 


Equation  (3-H)  is  precisely  Equation  (8.5-1)  with  Mp  =  0.  The  quasi -linear 
expressions  for  the  damping  exponents  are  those  given  in  Section  (8.5). 


h  =  ho +  + 


.  > • x3 ) 


^2  =  ^20  +  X22  +  Kg) 
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where  X .  = 

jo 


-  FT 
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•  XJ2  • 
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±M2 


20.  -  P 
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f  •  -V 
\  •  • 
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n'22/ 


The  zero  damping  curves  in  the  amplitude  plane  are,  therefore,  straight  lines 
with  slopes  -  1/2,  —  2. 

When  spin  is  not  zero,  the  origin  is  either  a  node  (l10^2o  >  °)  or  a 
saddle  (^0^20  ^  H1656  two  possibilities  are  shown  in  Figures  9-2a  and 

9-2b  for  a  linear  Magnus  moment.  For  a  cubic  Magnus  moment,  the  intercepts  and 
intersection  of  the  zero  damping  lines  must  be  studied. 

X10  +  X12(K?  +  2K|)  =  0 


X20  ”  X12^2K1  +  K§)  =  0 


(5.15) 

(3-lb) 


199 


Singularities  exist  at  the  intercepts,  (-  ^qA^-.  °)  ^^(0,  ,  and  at 


the  point  of  intersection 


t  X10 


+  2X, 


?n 


2Xiri  +  X?n 


-  3x, 


The  type  of 


12  12 
l0/^i2>  °)  can  now  * 

Equations  (3.13-14)  in  Equation  (3.4)  and  translating  the  origin  to  this  circular 


s  xn£^u.X&. 


rity  at  the  intercept  (-  X^/X-^,  0)  can  now  be  found  by  substituting 


singularity. 


UA 


2X.^  +  X^  -  X-,«  (2x  +  y) 
lu  .u: 


(p  “  ^1Q  ) 


"/V'-V 


(3.17) 


where  x  =  ^ 


X10 


v12 


A 

=  *6 


D  =  (b  -  c)  +  4  ad 


ad  -  be  = 


^X10  +  X20^  -  ° 


Xi0(2Xi0  +  >.20) 


(3.18) 

(3.19) 


For  these  relations  and  the  conditions  given  after  Equation  (3. 5),  we  see  that 
we  will  have  a  1  -  mode  circular  singularity  in  the  first  quadrant  if  X-^X-jp  <  0. 
This  singula  rity  is  s.  sstd.ci.ls  if  X-^q  ^  Ion )  ^  0  £tnci  is  8.  nods  ot* 


20 


XJC1  W  X£»C  ■ 


Similarly,  a  2  -  mode  circular  singularity  exists  if  420^22  =  "  ^20^1 2  <  °*  This 
singularity  is  a  saddle  if  Xor,(X1n  +  -Xon)  >  0  and  is  a  node  otherwise. 

fcU  xv/  t—\y 

The  character  of  singularity  at  the  point  of  intersection  (a  mixed  mode 
singularity)  can  be  found  by  transferring  the  origin  to  that  point. 


dv  (2 Xln  +  XorJ(2x  +  y)  -  3Xlp(2x  +  y)y 

vxjf  _  iu  c~\J  Xc. 

dv  ~  {\  +  Hv  4-  £-y)  +  [x  +  2v)x 

—  20  -'-'■TO'-  -*/- 


(3.20) 
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2  2 
.'.  D  =  33X.~  +  78x.^X„^  +  33X^ 

■  ■  XU  XU  d  U  d  U 


33 


v20 


"10 


von 


..812 


"10 


j 

T  • 


(3.21) 


ad  -  be  -  3(Xln  +  2\0  )(2\^q  +  Xon) 


(3.22) 


-2  />  *»  \ 
°  +  c  =  ^10  +  *20; 


(3.23) 


The  existence  and  type  of  singularity  can  now  be  determined  in  the  usual  way 


■t.r-?  ■f'Vj  ’f’Vjo  "Po  1  1  rM.ri*  r\  rr  v»aen  1  f  o  • 

n  x  Uil  uiu.  ivxj.wnxiig  i  u  l>  ux  u  u  , 


(1)  A  mixed  mode  singularity  exists  if  X,„  has  the  opposite  sign  to 

Xin 

X  and  -  2  <  —■  <  -  .5. 

'20 

(2)  This  singularity  is  a  node  if 


-  2  <  -i£  -  1.812  or  -  .352  ^  <  -  .5 
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^0 


(3)  This  singularity  is  a  spiral  singularity  if 


1  flia 

u.  •  j_n_ 


ID 

^?n 


or  -  1 


10 

^?n 


-  .552 


(4)  This  singularity  is  a  center*  if 


)10 

^20 


* 


To  prove  that  it  must  be  a  center, 
required . 


the  relations  of  Reference  9-1*-  are 
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Most  of  the  possible  amplitude  planes  are  shown  in  Figure  9 -3a-d.  In 
Figure  9-3 a,  the  curve  terminating  at  the  singular  point  is  called  the 
separatrix  and  for  a  missile  with  stable  small  amplitude  motion  it  divides 
all  initial  conditions  into  two  families:  those  which  yield  a  motion  tending 
to  zero  amplitudes  and  those  which  yield  motion  which  grows  without  bound. 
Figures  9- 3b  and  9- 3c  demonstrate  the  existence  of  circular  limit  motion  and 
epicyclic  limit  motion  respectively  while  Figure  9_3d  shows  the  very  special 
case  of  epicyclic  motion  with  periodic  variation  of  the  amplitudes.  The 
predictions  of  the  separatrix,  the  limit  circular  motion  and  the  limit  epicyclic 
motion  have  been  tested  by  numerical  integration  of  the  exact  differential 

.9-5 


equations  with  great  success 
the  quasi- linear  theory. 


This  success  reinforces  our  confidence  in 


9-12 


9-3c  Cubic  Static  Moment  -  Nonlinear  Damping  -  Zero  Spin. 

For  this  case,  we  will  study  a  very  general  form  of  the  aerodynamic  moment 
with  the  approximation  q  =  -  i£ 


C  +  iC 
m  n 


=  -  i  (c  + 

'  r\ 


c^5^  +  c*)£  +  d£ 


(3.2U) 


2  2 
where  d  =  d(5  )  is  a  function  of  6  and 

c*  =  c*((&  )  )  is  a  function  of  ( 6 c~) 

To  define  cQ  uniquely,  we  will  require  that  c*(o)  =  0.  The  presence  of  ( 5  )  in 
c*  indicates  that  this  term  probably  affects  the  damping  as  well  as  the  more 

i 

conventional  amplitude  dependent  damping  term  d£  . 

The  simplest  nonlinear  form  of  d  and  c*  can  be  obtained  by  assuming  them 
to  be  linear  in  6  and  (6  )  respect ively. 


C  +  i' 
ra 


2  *,  2.  • 

Co  +  C2S  +  C2^6  ^ 


£  + 


do  +  a2& 


|  (3.25) 


These  two  nonlinear  damping  terms  have  different  directions  and  so  depend  on 


the  eccentricity  of  the  quasi-iinear  elliptic  motion.  For  circular  motion, 
i9  » 

£  =  &Qe  and  £  is  perpendicular  to  £. 


C  +  iC  = 
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.  _  i© 
i8_e 
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2  t  2 

c  +  c^6  *f  iQ  (d.  +  5  ) 
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(3.26) 
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For  planar  motion,  however,  £  =  &e  and  i  has  the  same  direction  as  | 

ie  r  ?  ~ 

C_  +  iC_  =  -  ie  0  (c  +  c_8  )&  +  [d  +  (d0  +  20^^16'  (5*27) 

m  u  I  ■  O  2  u  c.  d.  ' 


Thus  the  damping  moment  depends  on  the  shape  of  the  motion.  This  result  casts 
some  doubt  on  the  standard  wind  tunnel  method  of  measuring  aerodynamic  damping 
by  forcing  planar  oscillation  and  determining  the  net  energy  required  to  maintain 
the  oscillation.  If  the  actual  damping  moment  has  the  form  of  Equation  (j.25), 
these  wind  tunnel  measurements  could  lead  to  quite  erroneous  predictions  of 
actual  flight  behavior. 

Returning  now  to  the  more  general  Equation  (3.24),  we  derive  the  corresponding 
equation  of  pitching  and  yawing  motion  for  small  geometrical  angles 


"  o  t  *  9  I  O 

t  J  TT/oMe  w  U  .  /u  .  M  _  n 
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where  H(B2)  =  -  CD  -  kt2  d 

„*  _  PSi  k-2  * 

M  '  55”  kt  c  * 

For  Equation  (3=28),  the  quasi-linear  damping  exponents  can  be  obtained 
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where 


d  =  h  +  l4nu(K?  +  K?)  +  nu(l2K?  +  21K?K^  +  12k!J) 
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=  (1  +  (K^  +  21^)  +  J  -Mo  (1  +  m2(2<  +  K^) 


Equations  (3-29-30)  are  symmetric  in  the  modal  amplitudes  and,  hence,  only 
half  of  the  amplitude  plane  need  he  studied  because  the  line  is  a 

line  of  mirror  symmetry.  We  will,  therefore,  consider  only  the  lower  pie  shape 
part  of  the  amplitude  plane  bounded  by  the  -  axis  and  the  line  and 

obtain  the  upper  half  by  symmetry. 


It  is  important  to  note  that  singularities  on  the  boundaries  of  this 

region  represent  simple  types  of  motion:  circular  and  planar.  We  will, 

therefore,  study  the  implications  of  Equation  (3-29-30)  in  tne  vicinity  of 

these  boundaries,  i.e.,  for  almost  circular  motion  and  for  almost  planar 

* 

motion.  Throughout  this  work  we  will  assume  H  and  M  to  be  differentiable 

•Tnnotl  nns . 


With  this  in  mind,  H  and  M  will  be  expanded  about  a  circular  motion 
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Equations  (3.31-32)  may  now  be  substituted  into  Equations  (3.29-30)  and  the 
results  simplified  for  almost  circular  motion  (K-  <  <  K^). 
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A  singular  point  exists  on  the  axis  at  =  Sc  when 


H  =  0 
"c 


nntinnc  (  W 


oiiLcixxd'  cii'L'uxoLL  iiiuu-LUiio  ^ K-^  <C  6^ )  grow  "t OwtLi'u  ttild  larger  circular 

motions  decay  when 
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(5.36) 


This  singularity  represents  a  true  limit  cycle  if  almost  circular  motion  tends 
to  become  more  circular  near  the  singularity,  ( X <  Q) 
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<  0  for  Kg  <  <  Kx  =  6c  (3.37) 


Inequalities  (3.56-37)  are  precisely  the  conditions  that  the  circular  singularity 
be  a  stable  node.  A  somewhat  more  convenient  form  of  these  inequalities  is 
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J  1  +  m 


A  positive  Xg  would  make  the  singularity  a  saddle.  If  we  define  an  unstable 
saddle  as  one  for  which  trajectories  on  the  origin  side  of  the  saddle  are 
directed  away  from  the  origin,  we  see  that  the  singularity  is  an  unstable 
saddle  if  Inequality  (3-36)  and  not  Inequality  (3.58)  are  satisfied.  Finally, 
a  negative  has  the  effect  of  reversing  the  stability  of  the  singularity. 


d8~ 

i-  ->c 

These  conditions  are  summarized  in  Figure  9-b. 

The  prediction  of  a  circular  limit  motion  when  conditions  (3.35-36,  38)  are 

satisfied  can  be  verified  by  direct  numerical  integration  of  the  actual 

differential  equations  of  pitching  and  yawing  motion.  This  has  been  done  for  M* 
2  1  2 

and  H  linear  in  (5  )  and  5  respectively  and  motions  like  that  of  Figure  9-5 
have  been  obtained. 
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One  very  interesting  result  can  be  obtained  from  the  Figure  and  the 

various  Inequalities.  If  the  damping  moments  are  only  functions  of  &  (M  -  0), 

every  circular  singularity  must  be  a  saddle,’  This  was  shown  in  our  first 

examples  to  be  the  case  when  H  was  quadratic  in  5  and  we  now  see  that  this  is 

2 

true  whenever  H  is  a  differentiable  function  of  o  and  the  aerodynamic  moment 

is  a  function  of  only  5.  Under  these  conditions,  a  circular  limit  motion  is 

impossible.  Since  circular  limit  motions  have  been  observed  in  ballistic  range 

* 

tests,  a  more  complicated  moment  such  as  one  containing  terms  like  c  £  must  be 
present. 

Turning  now  to  the  question  of  almost  planar  motion  (K^  =  K^)  we  will 
consider  motion  in  vicinity  of  the  planar  motion  K,  =  Kp  =  K^. 


6^  =  2K*(1  +  cos  0) 

P  P 


(3.39) 


(62L  = 


Q  A  .  A 

2 0  sin  0 
P  P 


(3.40) 


where 


0’  =  J- mo(4  +  3m J 


m 


. 


M 


The  nature  of  the  almost  planar  motion  near  this  planar  motion  can  be 
conveniently  studied  by  use  of  the  variables  and  which  are  defined  by 
the  equations: 


x r  t  ^  —  \ 

Kj  =  vx  +  V 


A  T  t 

J  =  1 


(3.41) 


Since  onlv  motions  near  this  planar  motion  are  of  interest,  the  e.’s  will  be 
„  -  j 

assumed  to  be  small  in  comparison  with  unity  and,  hence,  higher  powers  of  the 

e  ’s  will  be  neglected  in  comparison  with  their  linear  powers  or  unity. 

3 

2  a  f  2  t 

m 4- V.  «  +  4  /n-P  C  1^  ao  vn  "Ka  a/1  «  a  n  nmn  1  1 

liiciciuicj  one  vai  xauxouo  0  y  'p  cuui  ^  u  )  uuxi  ue  oouij^uucu  ao  a.  S  uku.  J _ L 

2  a  |  2  t 

perturbation  of  the  corresponding  variations  in  6p>  0p  and  (B  )  . 
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(3.42) 


6  = 


5* 

0 


5P  +  +  008  ^€1  +  e2^ 

<-0; 


1 

5 


/ 


f-M 
v  o 


4  +  m  (3  +  2e1  +  4e2) 


/_U 

*»  /"n 

V  O 


4  +  m  (3  +  4t.  + 

p  1 


o  N 


=  d' 

rv 


3  mU. 

1  +  — P— x- 


_  \ 
*oJ 


8  +  8m 


,2.'  ,  2 .  t  Pa  1  a/®  +  \ 

is-)  -  (»') ?  i  m^l  (£i+e2> 

■  p' 

H  and  M  can  now  be  expanded  about  the  planar  motion. 


(3-43) 


(3-44) 


K  =  Hp  + 

"  &H  ‘ 
“5 

d6  J 

i — 1 
a> 

■  ro 

1 

S] 

P 

*  _ 

ib 

— 1 

s: 

—I 

M  =  M  + 

2v  •  1 

(s 

P 

La^°  )  J 

P 

- 

t~>  niV\  0  n  ■y'  "5  r\ 
J J  OUboCI  Ip 

t  on 

a  function 

of 

.2s* 


p  t 

n 

■  'p 


/.2n» 


(3.45) 


(3.46) 


2  2  2  *  2  * 

evaluated  for  =  6^  or  (&'“)  =  (6* “ )p.  Relations  (3.42-46)  together  with  the 

proper  expansion  of  02 K^/0  and  W  8e  placed  in  the  definition  of  X^ 

with  the  result 
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resu. 


Ow 

C.31 


i  r 


r_.  ( 4  +  2mn  \ 


H-K-K  I  <e!'  V  co6? 

yj  '-  -  \  -  -n  / 


+  2 


dH 


LdS^J 


K^(e,  +  e„)  sin^  0 
P‘  i  Z' 


2Mp  +  mp)£l  -  +  mp^e2 


M 

sin  0 


0* (4  +  3m_) 


r-  *  -i 

AM 


.  77Z 2T7. 

La^°  )  J 


?  f  8  +  ^  ?  *,  I  a 

^  [imrJ  i£i +  e2}  sin_  ^ 

■  p 


(3.47) 


<ijt 

*  i  r 
vhere  ip--EJ 

O 


!  2  sin  3  \ 


(H„)(l  -  cos  (if)  +  (Mn)  (  J 


■n '  \  'S 1 

-  N  0  , 


d0. 


The  terms  involving  ■ — 


r^rl  and  M  nay  he  integrated  hy  parts  and  the 

Klp 

1  fe  r*ATnV\-!  t.H  f  Vi  +Vi£»  a4-Viqv«  4-o  vwir 

xuu  uwiuL/xuvyU  irx  ua  vuu  uumioi 

,  ,  0],  rr 

h  ■  XP  ■  m  ;>'  j  \  L(8  +  5mp)ei +  "pe2. 


m  p  4-  flf)  +  11m 

“p~l  '  '  - P'~2 
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2n 

i  r  a,  a 

[H]  =  -  I  cos  <fi  dfo  . 

1  *  J  P 

o 

A  similar  relation  for  Xg  may  be  obtained  by  interchanging  and  in 
Equation  (3.48). 


The  numerical  subscript  on  the  outside  of  the  bracketed  expressions  in 
X1  identifies  that  term  as  a  particular  Fourier  cosine  coefficient.  It  is 
quite  surprising  that  the  influence  of  H  on  X,  is  completely  determined  by 
its  first  order  Fourier  cosine  coefficient.  The  influence  of  M  ,  however, 
is  specified  by  the  zeroth  and  second  order  Fourier  cosine  coefficients  of 
its  first  derivative.  These  coefficients  are  computed  for  fixed  modal  amplitudes. 


Kj,  and  are  thus  functions  of  these  amplitudes. 

* 

For  a  planar  singularity  X  vanishes  and  the  X.  are  linear  in  the  e  . 

P  0  J 

only  planar  value  of  the  coefficients  of  Equations  (3.29-50)  are  needed  to 


Thus, 


U  I  ft  (j  C 


*  # 

X,  =  a  X  +  a  X 
J  l  1  Jo  2 

2(4  +  3m  )  (8  +  701^) 

Where  a^  =  ag2  =  («  +  9m^(8  +  5o  ) 

-  4m  (4  +  3m  ) 

&12  “  a21  *  TB  '+  5m J 

P  P 

[H]  1  [-1  + 

*1  “  "  2(B  +  5"b)(B  +  9®n) 

[HJ  [»«i  +  b£2' 

X2  x  “2(8  +  5ml) (8  +  9nt_) 

p  p' 


(3-49) 


f  *  cxoV 
W  ^>  ) 


(3.51) 
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p 

where  a  =  -  m  (8  +  3m  )  +  (128  +  200m  +  75m  )  r 

P  P'  K  P  p'  o 


-j.~.  I Q  i  c—  N  _ 

pm  t  ?Ui  J  I  o 
p  p  2 


b  =  64  +  96m  +  33m  +  3m  (8  +  5m  )  r 


P  o 


+  (128  +  248m  +  105m'  )  r_ 

p  '  p'  2 

The  differential  equation  for  solution  curves  in  the  vicinity  of  a  planar 
singularity  in  the  amplitude  plane  is 


de  X 

2  .  2 


Q  /•  X  Vv/» 
“£1  T  uc2 

bel  +  ae2 


f  X  CO  \ 


According  to  the  usual  criteria,  the  singularity  must  be  either  a  saddle  or  a 

2  2  2  2 
node.  It  is  a  node  if  a  -  b  is  negative  and  a  saddle  if  a  -  b  is  positive 


a2  -  b2  =  (a  +  b) (a  -  b) 

=  -  4(8  +  5mp)(8  +  9m_) 


(3.53) 


x  |  2(2  +  m  )  +  (r  +  r0)(8  +  5m.) 

L  P  0  ^  P' 

x  (4  +  3mp)(l  -  2ro)  +  4(2  +  5mp)r2 


Note  that  if  rQ  and  rg  both  vanish  and  mp  is  outside  the  interval  (-  2,  -  8/9), 
a2  -  b  is  negative.  Thus,  if  the  aerodynamic  moment  coefficients  are  function; 

2  2  1 

of  S  alone  and  not  functions  of  )  ,  all  planar  singularities  are  nodes  and 
almost  planar  motions  near  a  planar  singular  motion  will  tend  to  that  planar 
singular  motion  if  neighboring  planar  motions  tend  to  the  planar  singular 
motion. 


Another  interesting  special  case  is  that  for  which  H  and  M  are  linear  in 

2  ,2s’ 

5  and  (5  ;  respectively. 


(3-54) 


H  =  H  + 

O  c. 

*  */  2v 1 

M  =  M^B  ) 

A  simple  calculation  shows  that 
[H],  =  2 I^K2 

**  ^  F 


(3-55) 


(3-56) 


r 


n 


(3.57) 


r2  =  ° 


(3.58) 


The  amplitude  of  the  planar  singular  motion  is  fixed  by  the  condition  X  =0  to 

P 

be 

Kp  =  -  Ho(K2  ‘  2^)_1  =  “  HoH2’1(1  "  2ro)_1  >  0  (3.59) 

2  2 

For  r2  =  0,  a  -  b  reduces  to  a  simple  form.  A  planar  singularity, 
therefore,  will  be  a  node  if 


(°  +  5mp)(°  +  9mp)(1+  +  3mp)(l  -  2rQ) 


U  +  2m_  +  r  (8  +  5 )  >  0 

v  °  F 


(3.60) 


When  m^  is  excluded  from  the  interval  (-  2,  -  8/9)  and  rQ  is  replaced  by 


its  actual  value,  this  Inequality  can  be  written  in  a  very  concise  form. 


4  +  2m 

_ P 

8  +  5m 

^  r\ 


*2 


<  H0  <  2 


(3.61) 


According  to  Inequalities  (3.61),  (1  -  2rQ)  has  to  be  positive  and,  hence. 
Equation  (3-59)  shows  us  that  Hq  and  Hg  must  have  opposite  algebraic  signs. 
A  negative  Hq  causes  small  amplitude  motion  to  grow  and,  hence,  the  planar 
node  will  be  a  stable  node  if 

H0  <  °»  >  0 

i.ne  VCU1ULU)  1U111  tlCU  tti'tr  m  uiiumir  l/,wi  ±11  T  J-gllTe  ^-O  • 
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(3.62) 


9.4  The  Perturbation  Method 

The  work  of  this  chapter  has  assumed  that  epicyclic  motion  which  is 
essentially  a  linear  combination  of  trigonometric  functions  is  a  good 
approximation  of  the  actual  motion.  The  value  of  this  approximation  should 
diminish  as  the  static  moment  becomes  strongly  nonlinear.  In  Section  8.3, 
the  exact  solution  for  a  cubic  static  moment  was  derived  in  terms  of  elliptic 
functions.  In  this  section,  we  will  use  this  solution  rather  than  the  epicyclic 
solution  as  our  basic  solution  and  consider  the  effects  of  both  linear  and 
nonlinear  damping  to  be  perturbation  of  this  elliptic  function  solution*.  The 
discussion  will  not  contain  all  the  necessary  details  but  will  be  more  in  the 


For  simplicity,  we  will  neglect  the  geometrical  terms  7  / 7,  assume  all 
coefficients  to  be  independent  of  arclength**,  and  require  constant  spin. 

The  approach  will  be  to  manipulate  Equation  (2.4)  to  obtain  differential 
equations  for  the  energy  and  angular  momentum.  First,  we  transform  Equation  (2.4) 
by  the  relation 


Ai  (1/2  IPs 
i  =  i 


(4.1) 


.  a"  Ta  p"l  a  ai  r  *  ,  hIa 

. .  i  -  Mq  +  1^5  i  =  -  Hi  +  M  +  iP(T  -  p  I 


(4.2) 


where  M  =  M  -  -p- 
o  o  4 


M*  =  M  -  M  -  M„62 

O  ^ 


The  conjugate  of  Equation  (4.2)  is 


a11  a  2a  ai  *  ,  n.j  * 

I  -  Mq  +  Mg5  !  =  -  H!  +  M  -  iP(T  =  |)  l 


f  \.  -7  \ 


This  perturbation  method  is  an  extension  and  modification  of  the  work  of 


Reference  9”6  and  9“7 


The  effect  of  varying  coefficients  is  considered  in  Reference  9-8. 
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Multiplying  Equation  (4.2)  by  §  ,  Equation  (4.3)  by  i  .  and  adding,  we  have 
an  equation  for  the  change  of  total  energy. 


C1  =  -  H 


a  2  4 

2C1  +  2MqS  +  MpB 


+  (M*)  (&2) 


+  P(T  -  |)  C2 


tu  h) 

V  _r*  "T  / 


.  A  t  .  O  A  O 


where  Cn  =  j£  j  -  M^^-Mp  75-  total  energy 

~tA  All 

Cg  =  i(5  5  -  5  5)  angular  momentum 


The  equation  for  the  variation  of  angular  momentum  may  be  obtained  by 
a  2 

eliminating  M  +  Mpb  from  Equations  (4.2-3). 


I  II  p 

Tin  i  o-n/m  11 N  cc 

no0  -r  cr\i  -  57  u 


p  — 
w2  " 


(4.5) 


When  damping  is  zero  (H  =  T  =  M  =0),  the  amplitude  variation  is  given  by 
Equation  (8.3-9)  which  is  repeated  here 


(62)'  =  -  J -  C2  +  4C162  +  4Mq6^  +  2^Bfc 


(4.6) 


The  actual  solutions  in  terms  of  the  sine  amplitude  function  are  given  in 
Section  8.3.  The  requirement  of  rotational  symmetry  once  again  leads  us  to’ 

O  PI 

the  restriction  that  the  aerodynamic  coefficients  are  functions  of  6  (5-) 

it,?  o»  .t.o  p  >  \  /  y 

and  1 1  j“.  But  (b1-)  and  |£  p  are  functions  of  C., ,  0o  and  b^  through 

-A-  I — 

Equation  (4.6)  and  the  definition  of  respectively.  Thus  Equations  (4.4-5) 
have  the  form 


C,  =  f,(C. ,  C„,  8*) 

J  J  1  27  - 


J  =  1.2 


(k. 

\  T  *  I  / 


For  a  conservative  system  f j  =  0,  and  Equation  (4.6)  applies.  We  assume 


that  the  damping  funetious,  f  ,,  represent  small  damping  and,  hence,  C,  and  C0 

J  X  c. 

vary  very  little  over  a  period  of  the  nutation.  If  this  is  the  case,  it  is 

2 

reasonable  to  use  the  no  damping  solution  for  5  in  f  and  to  assume  that  on! 

2*  J 

the  average  effect  of  the  b  is  important. 
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^  = 


u  J  a 


(4.8) 


«* 

i_  P 


o 

o 


■P  He  or*H 

"J  ““  — 


*  2 
P  is  a  period  of  S'-. 

c.2  _  -, 

Since  the  solutions  for  o  are  usually  given  in  terras  of  the  extreme  values 
and  5^,  this  average  function  will  contain  8^  and  5^  as  well  as  and  Cp . 
These  two  sets  of  variables  are  related  by  Equation  (4.6)  and  the  fact  that 

2  t  2 

(6  )  vanishes  for  extremum  values  of  8  . 


40,5^  +  4M  51!  +  2H.Sb  =  C?  j  =  1.2 

1  j  o  j  2  j  4 


(4.9) 


Thus,  either  set  could  be  eliminated  in  favor  of  the  other.  In  view  of  our 

previous  work  with  the  amplitude  plane,  a  third  set  is  much  more  convenient. 

These  variables  are  generalized  squared  modal  amplitudes  and  will  be  defined 

2 

so  that  the  sum  of  their  square  roots  is  the  maximum  amplitude,  5  ,  and  the 
difference  is  the  minimum  amplitude.  Instead  of  using  the  previous  symbols 
K*7,  as  was  done  in  Equations  (8.3*14-15),  ve  will  use  new  symbols,  x^,  so  that 
any  possible  confusion  with  the  epicylic  modal  amplitudes  can  be  eliminated. 


(n/  x±  -  -J  x2Y 


(4.10) 


(S 


(4.11) 


According  to  Equations  (4.10-11),  the  x,  axes  coi'i-espond  tu  constant  amplitude 

J? 

circular  motion  in  the  same  way  as  the  axes  did  in  the  quasi-linear  method. 


Equations  (,4.9-11/  may  now  be  combined  to  relate  the  Xj*s  and  CjSs. 

A  “  t  p  2 

2Mo(x1  +  x2)  +  -g-  (5x“  +  lOXjXg  +  5x~) 


C1  “  - 


(4.12) 


C2  =  2(Xl  -  x2)  J-  Mg  -  M2(x1  +  x2) 


(4.13) 
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According  to  Equation  (4.13),  the  angular  momentum  is  zero  when  x-j  =  x?  and, 
hence,  this  line  retains  the  locus  of  planar  motion*.  Equation  (4.13)  indicate 
a  second  possibility  for  planar  motion,  i . e . , 


ft 


Vxi 


+  3C2)  = 


u 


(4.14) 


An  inspection  of  the  various  solutions  of  Section  8.3  shows  that  this  can  occur 
for  the  type  (c)  moments  and  represents  oscillations  about  the  trim  equilibrium 
angle  of  J  *  /Xl  For  the  type  (c),  the  line  given  by  Equation  (4.14)  is  also 


V*? 

the  locus  of  planar  motions. 


If  Equations  (4.12-13)  are  differentiated  and  solved  for  x , ,  the  fine 

J 

equations  for  the  generalized  amplitude  plane  follow. 


V  — 
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2M 


‘  ‘32 
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2M 


1/  t  c,- 

‘'*2  v  ^2. 


+  3x0) 


J'  " o  -  *2  (X1  +  *2)C2 


Qa>2&2 


(4.15) 
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2Mq  +  ^2  (^xi  +  x2) 
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! — ~ 


o  +  ^2  (5xl  +  5x2^  J~  Mo  ‘  **2  ^X1  +  X2^C2 


XT 


-  8cs  to 


(4.16) 


2  a 

where  03  =  -  M  - 

o  2 


■-  w 

03  =  -  M  - 


3*1  +  2  +  5x2 


[o  ■  *2 


and 


C  *s  are  given  by  Equation  (4.8). 

u 


As  an  example  of  this  perturbation  method,  we  will  consider  the  special  case 
of  Section  9- 3c. 


This  motion  is  planar  in  the  roofed  coordinates  specified  by  Equation  (4.1). 
Only  when  spin  is  zero  will  the  motion  be  planar  in  the  eyes  of  a  nonspinning 
observer. 


<"'»  i  /T 

^XO 


(4.17) 


H  =  Ho  +  H25 


*  *  .  2. » 

M  =  ^(6  ) 


f  1 ,  T  '.J  \ 

,10  j 


P  =  0  . 

Equations  (4.8)  reduce  to 


(4.19) 


1  f 

/  *.  2 

A  2  4 

C.  =  -  i 

H  +  (H.  -  2M_ )5 

2C,  +  2M  5  +  M_5 

J-  L 

o  '12  12 ' 

1  o  *2 
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lvl2^2 


cn  =  - 

CL 


Ho  +  V. 
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f  I  on  \ 


Before  using  these  equations  in  Equations  (4.15-16),  we  must  compute  the  average 

value  of  S^11  for  n  =  1,  2,  3.  In  general,  this  is  quite  tedious  since  it  involve 

2n 

the  calculation  of  integrals  of  sn  cus .  This  can  be  done  numerically,  however, 

and  Equations  (4.15-16)  may  then  be  integrated  to  construct  a  generalized 

amplitude  plane  for  particular  values  of  the  coefficients.  To  derive  general 

conditions  on  the  coefficients  for  the  various  limit  motions,  it  is  necessary 

2n 

to  locate  singular  points  and  calculate  the  6q  *s  in  their  vicinities.  For 
two  special  cases  -  circular  singularities  and  planar  singularities  -  the 
necessary  algebra  is  not  excessive.  From  this  partial  information  about  the 
complete  generalized  amplitude  planes,  it  is  then  possible  to  infer  the  general 
form  of  these  planes  and  select  interesting  sets  of  values  of  the  coefficients 

for  numerical  construction  of  specific  generalized  amplitude  planes.  For 

/  .  9-9 

almost  circular  motion  (x^  >  !>  Xg),  it  can  be  shown  that 


e2n  .  n  n-1  (  mc  \ 

6a  *  *1  +  ”*1  x2  1  n  '  5  '+  3m 

\  c  / 


(4.22) 


Equations  (4.20-22)  may  now  be  placed  in  Equations  (4.15-16). 
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(4.25) 


n2  * 


Hg  +  +  pinj  Mg 


jr 


r-T 
CLX  ( 
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Equation  (4.25)  is  exactly  the  same  as  Equation  (3.33)  which  was  derived 
from  the  quasi-linear  method.  The  symbols  of  that  equation  may  be  easily 
related  to  those  of  Equation  (4.23)  by  a  consideration  of  the  coefficients 
defined  by  Equations  (3.31)  and  (4.17)  respectively 
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(4.26) 


Equation  t.4.24;  and  the  corresponding  quasi-linear  Equation  (3.3*0  do  differ 
markedly.  The  importance  of  this  difference  can  be  most  conveniently  evaluated 
by  comparing  the  two  sets  of  conditions  for  a  circular  singularity  which  is  a 
node.  The  quasi-linear  conditions  may  be  written  in  the  form 


i  -  -  hcA>  > 0 


(4.27) 
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(U.29) 


1  +  m 


Relations  (4.27-28)  axe  the  same  as  those  derived  from  the  perturbation  method. 

1  1  >'V\0  4“  T  /'>*'>  W1 A+ Vl  f  /“(f1  Tr.SrtJSfi  l  -*  -f---  f  ! !  fj  V. 


Th 
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turbation  method  form  of  Inequality  (U.29)  can  be  derived  by  setting 


X-L  =  0  and  x 2  <  0,  and  eliminating  Hq  between  Equations  (4.23)  and  (4.24) 


4  -  h  +  mc) 

\  <  2  +  >c 


(4.30) 


For  a  linear  static  moment  (m  =  0),  Inequalities  (4.29-30)  are  the  same. 
Their  asymptotic  values  -  (l  +  s[2)~±  and  (-l/3),  however,  differ  by  35$.  For 
a  type  (b)  moment,  the  quasi-linear  expression  limits  in,  to  the  interval 
(-  l/2,  0)  while  the  perturbation  expression  is  valid  for  the  complete  range 
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in  m^  for  which  almost  circular  motions  are  possible  (-  2/3  <  mc  <  0).  These 
differences  and  others  may  be  seen  in  Figure  9=4  where  both  upper  bounds  a re 
plotted  as  functions  of  m  . 

The  algebra  for  almost  planar  motion*  is  more  lengthy  and  will  only  be 
outlined  here.  If  the  common  value  of  x^  for  planar  motion  is  identified  by 
x  ,  the  maximum  amplitude  of  the  motion  is  2  */x  .Our  generalized  definition  of 


m  is 
P 


^XP 


(4.51) 


and  our  generalized  e.'s  of  Equations  (3.41)  are  defined  by** 

J 


v  =  -v  ( ~\  ±  2e  ) 

“j  'V"  '  J 


i  -  1  o 
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(\ j 

S'-  I 


For  these  variables,  the  average  values  of  &  can  be  determined 


S,  +  2B2n(£l  +  e2>]  (1*Xp,n 


(4 .33) 


where  A2n  and  B2n  are  combinations  of  complete  elliptic  integrals 
of  the  first  and  second  kind  which  in  turn  are  functions  of  m^.  These  quantities 
are  explicitly  derived  and  defined  in  Reference  9-10. 


The  amplitude  of  the  planar  singular  motion  2  \/x  is  fixed  by  the  condition 
t  P 

x  .  =  0  to  be 
J 
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(4.34) 


where 


A  =  A(m  )  = 
P' 


+  m  - 
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+  m  )A0  -  c 

p"  2 


b,  ■  B  A/ 
4  p  6 


Only  the  symmetric  planar  motion  for  which  =  x~  will  be  considered  here. 
Planar  motion  about  trim  is  studied  in  Reference  9-10. 


.  „2  __a 


The  2  appears  in  Equation  (4.32)  because  Xj  coi-x-es ponds  to  KTJ  and  1  +  2e j  is 
the  first  order  approximation  of  (l  +  ej)2  • 
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The  corresponding  equation  for  the  quasi-linear  method  is  Equation  ( 3 . 59 ) 
which  differs  from  the  above  by  the  absence  of  A.  One  interesting 
characteristic  of  Equation  (4.34)  is  that  for  certain  values  of 
and  H-,/m  and  a  type  (c)  moment,  it  has  two  roots.  The  m^'s  for  the  members 

4_  O  p 

of  such  a  pair  lie  in  the  intervals  (-  2.535,  ~  2.164),  and  (-  2.164,  -  2) 

respectively.  If  we  omit  consideration  of  the  second  member  of  these  pairs 

(m  <  -  2.164,  m  >  -  l),  the  conditions  for  a  stable  node  may  be  quickly 
P  P 

stated. 
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(4.36) 


The  lower  bound  in  Inequality  (4.36)  differs  from  that  of  Inequality  (3,6l), 

I  pVl  Q  4*1  »/>  1  rvT  mn  ^  -t  m  "0-?  rfi  O  K 

u. iiC.  U*iu  t-fcu.  c  v-wiu^jcij.  cu.  XAl  f  l^Uic 

Now  that  we  have  considered  the  implications  of  singularities  on  the 

coordinate  axes  and  the  line  —  Xg,  we  want  to  show  how  information  about 

the  remainder  of  the  generalized  amplitude  plane  can  be  inferred.  First,  it 

is  necessary  to  combine  the  various  boundaries  for  in  one  Plot'  To  do 

this,  we  must  compute  m  in  terms  of  m  . 

c  p 

From  Equations  (4.27,  31,  34), 
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By  means  of  Equation  (4.39),  Figure  9-6  can  be  replotted  in  terms  of  mc  and 
superimposed  on  Figure  9-4'  The  result  is  Figure  9-7» 
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This  figure  has  a  number  of  interesting  features.  For  a  Type  (b)  moment, 

*  / 

sen  that  as  M^/H 2  is  increased  a  circular  node  and  planar  saddle 
change  to  a  circular  saddle  and  a  planar  saddle  which  in  turn  change  to  a 
circular  saddle  and  a  planar  node!  This  disappearance  of  the  node  from  the 


_  QV1  . 


-L.  o  eu 


o 1  q+pv  o or'flnr'p  on  +  Vi 1  i  np  v 

ouiu.  wiiwo  wit  XXUv.  /v  ^ 


implis 


i  + 


the  generalized  amplitude  plane  and  the  possible  existence  of  a  limit  motion  v 

between  two  nonzero  modal  amplitudes.  This  conjecture  was  checked  by  the 

numerical  calculation  of  the  generalized  amplitude  plane  for  m  =  -  ,55  and 
*  c 
M^/tu  =  -  1.  (Figure  9-8).  This  figure  demonstrates  the  existence  of  the 

stable  node.  The  differential  equations  of  pitching  and  yawing  motion  for 

these  parametric  values  Were  then  calculated  on  an  analogue  computer  and  it  was 

found  that  the  motion  quickly  went  to  a  limit  motion  with  constant  maximum  and 

minimum  amplitudes.  (Figure  9- 9a)  These  maximum  and  minimum  amplitudes  are 

those  assigned  by  the  location  of  the  node  in  the  generalized  amplitude  plane. 


Finally,  we  note  that  the  dotted  region  of  Figure  9-7  appears  to  be 
quite  interesting  due  to  the  presence  of  two  stable  limit  motions.  A  sample 
generalized  amplitude  plane  is  shown  in  Figure  9-10  and  a  plot  of  the  angular 
motion  approaching  the  planar  motion  is  given  in  Figure  9-9h.  The  ability 
of  the  perturbation  method  to  predict  fairly  unexpected  motions  is  very 
impressive,, 


CHARLES  H.  MURPHY 
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prediction  of  missile  flight  performance.  Finally,  the  influence  of  strongly 
nonlinear-  static  moments  is  determined  by  a  perturbation  method  which  makes  use 
of  two  quasi -constants  of  the  motion  —  total  energy  and  angular  momentum. _ 

-AE; _ _  Accession  No. _  USCLAS3IFIED 

Ballistic  Research  laboratories,  AFG 

FREE  FLIGHT  MOTIOM  OF  SYMffiTRIC  MISSILES  Symmetric  missiles  - 

C.  H. Murphy  Motion 

I  Drag  and  rolling 

BRL  Report  No.  1216  July  1963  motion 

Linear  coefficients 

IRDT  &  E  Project  No.  1MD10301A005  Cubic  coefficients 

lUMCIAJ55IFIED  Report  Measurements  of 

coefficients 

The  linearized  angular  motion  of  a  symmetric  missile  is  developed  in  some 
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misaligned  control  surfaces  is  briefly  considered.  This  linear  theory  is,  then, 
japplled  to  the  analysis  of  ballistic  range  data. 
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Inonlinear  terms  arc  derived  and  their  use  in  the  construction  of  "amplitude" 
planes  indicated.  These  amplitude  planes  Iravo  proven  to  be  quite  useful  for  the 
prediction  of  missile  flight  performance.  Finally,  the  influence  of  strongly 
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The  linearized  angular  motion  of  a  symmetric  missile  :1s  developed  in  some 
detail  with  some  consideration  of  the  rolling  motion,  drag,  aerodynamic  Jump  and 
the  effect  of  varying  coefficients,  The  tricyclic  motion  of  a  missile  vith 
misaligned  control  surfaces  Is  briefly  considered.  This  linear  theory  is,  then, 
applied  to  the  analysis  of  ballistic  range  data. 
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nonlinear  terms  are  derived  and  their  use  in  the  construction  of  "amplitude" 
planes  indicated.  These  amplitude  planes  have  proven  to  be  quite  useful  for  the 
prediction  of  mlanlle  flight  performance.  Finally,  the  Influence  of  strongly 
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of  two  quasil-constants  of  the  motion  --  total,  energy  and  angular  momentum. 
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The  linearized  angular  motion  of  a  symnetric  missile  la  developed  in  some 
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j the  effect  of  varying  coefficients.  The  tricyclic  notion  of  a  missile  with 
misaligned  control  surfaces  is  briefly  considered.  This  linear  theory  is,  then, 
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nonlinear  terms  are  derived  and  their  use  in  the  construction  of  "amplitude" 
planes  indicated.  These  amplitude  planes  have  proven  to  be  quite  useful  for  the 
prediction  of  missile  flight  performance.  Finally,  the  influence  of  strongly 
I  nonlinear  static  moments  is  determined  by  a  perturbat ion  method  which  makes  use 
I  of  tvo^  quasi-constants  of  the  motion  --  total  energy  and  angular  momentum. 
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The  linearized  angular  motion  of  a  symmetric  missile  Is  developed  In  some 
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The  linearized  angular  motion  of  a  symmetric  missile  is  developed  in  some 
detail  with  some  consideration  of  the  rolling  notion,  drag,  aerodynamic  Jump  and 
the1  effect  of  varying  coefficients.  The  tricyclic:  motion  of  a  missile  with 
misaligned  control  surfaces  is  briefly  considered.  This  linear  theory  is,  then, 
japplled  to  the  analysis  of  ballistic  range  data. 
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